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ABSTRACT 


This  is  a  final  report  of  a  detailed  and  quantitative  study  of  the 
influence  of  changes  in  specific  defect  concentrations  on  the  microstructure 
and  the  electrical,  optical  and  mechanical  properties  of  Si-doped  GaAs . 
Preliminary  evidence  for  the  GaAs:Si  system  has  indicated  that  qualitative 
correlations  between  Si-related  defect  concentrations  and  a  diverse  set 
of  physical  properties  are  observed.  The  measured  changes  in  properties 
are  large  and  can  be  employed  to  construct  appropriate  defect  models. 

GaAs  is  a  material  of  interest  to  the  electronics  community  as  it  is  used 
in  a  number  of  devices.  The  role  of  defects  in  affecting  transport  properties 
is  of  fundamental  importance.  The  goal  of  this  research  is  to  obtain  a 
more  complete  understanding  of  the  defect  structure  of  Si-doped  GaAs 
through  a  quantitative  correlation  of  various  physical  properties  and 
concentrations  of  Si-related  defects  and  also  to  track  the  effect  of  defect 
species  such  as  solute  vacancy  pairs  on  several  different  physical  properties. 
The  final  report  summarizes  the  work  of  a  five  year  study.  It  also  includes 
a  description  of  the  work  of  Professor  J.  M.  Whelan  and  Mr.  John  Graham 
to  determine  total  Si  concentrations  which  is  necessary  to  determine  the 
absorption  cross  sections  of  all  of  the  Si-related  defects. 
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I.  INTRODUCTION 


The  compound  semiconductor  GaAs  Is  of  great  technological  Interest 
because  cf  Its  considerable  application  in  electronic,  optical,  and  opto¬ 
electronic  devices.  In  many  of  its  applications  the  defect  structure  and 
properties  are  of  major  importance  and  this  importance  has  stimulated  basic 
studies  of  the  defect  structure  of  doped  GaAs  and  its  influence  on  physical 
properties.  Our  interest  in  this  topic  led  to  the  original  proposal  which 
was  to  initiate  a  detailed  and  quantitative  study  of  the  influence  of 
changes  in  specific  defect  concentrations  on  a  diverse  set  of  physical 
properties  for  a  number  of  doped  crystal  systems.  One  system,  GaAs:Si, 
was  selected  for  the  major  effort  and  we  have  now  completed  this  study. 

As  will  be  discussed  most  of  the  original  objectives  have  been  achieved. 

One  major  task,  that  of  developing  an  inexpensive  and  accurate  method  for 
measurement  of  the  total  Si  concentration,  still  remains  incomplete  and  we 
anticipate  completion  by  the  end  of  this  summer  (1981). 

In  the  original  proposal  it  was  indicated  that  preliminary  evidence  for 
the  GaAs:Si  system  shows  that  qualitative  or  semi -quant! tative  correlations 
between  Si-related  defect  concentrations  and  a  number  of  physical  properties 
are  indeed  observed.  The  changes  in  properties  produced  by  annealing  can  be 
large  and  these  changes  are  useful  in  helping  to  construct  appropriate 
defect  models.  Electrical,  optical,  mechanical  and  microstructural 
characteristics  are  all  part  of  the  system  of  physical  properties  which  were 
measured.  To  our  knowledge,  this  study  is  the  first  one  in  which  the  changes 
In  quantities  such  as  carrier  density.  Infrared  absorption,  vacancy  loop 
size  and  density  and  yield  stress  caused  by  changes  in  thermodynamic 
conditions  are  correlated  and  all  explained  on  the  basis  of  a  common  defect 
model . 
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There  have  been  a  number  of  earlier  studies,  primarily  at  the  University 


of  Southern  California,  which  demonstrated  that  aging  produces  dramatic 
effects  on  the  properties  of  heavily  Si-doped  GaAs .  The  aging  effect 
is  observed  for  specimens  in  which  the  Si  concentration  is  less  than  the 
high  temperature  solid  solubility  limit.  The  effects  included  changes  in 
the  free  carrier  concentration,  defect-induced  localized  vibration  mode 
absorption,  photoluminescence,  TEM  measured  microstructure,  and  critical 
resolved  shear  stress.  The  results  of  these  earlier  investigations  suggested 
that  the  annealing-induced  changes  in  these  various  properties  are  inter¬ 
related.  Indeed,  in  the  present  study,  a  qualitative  correlation  has  been 
found  between  changes  in  these  properties  and  changes  in  the  concentration 
of  Si-vacancy  defects  such  as  (Si gaVQa )  pairs,  where  indicates  a  vacant 
M  site.  The  difficulty  with  detailed  comparisons  of  the  previous  measurements 
is  that  they  were  performed  using  material  from  different  ingots  or  different 
sections  of  an  ingot,  and  therefore,  the  Si  concentrations  were  different. 

The  primary  goal  of  this  research  was  to  determine  whether  a  quantitative 
correlation  exists  between  physical  properties  and  the  concentrations  of 
lattice  defects  such  as  (Slgavga)  pairs  in  heavily  Si-doped  GaAs.  This 
primary  goal  was  achieved  in  that  the  correlations  were  found  and  the 
details  are  given  in  the  body  of  this  report. 
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II.  RESEARCH  OBJECTIVES  AND  RESULTS 


A.  ABSTRACT 

The  purpose  of  this  study  is  to  measure  and  correlate  the  annealing- 
induced  changes  in  the  carrier  density,  the  microstructure  and  the  yield 
stress  for  different  Si  concentrations  In  GaAs.  The  measurement  which  permits 
the  correlation  of  these  properties  is  the  infrared  absorption  of  the 
localized  vibrational  modes  (LVM).  As  described  in  the  following  text 
the  LVM  measurements  allow  one  to  track  the  annealing-induced  changes  in 
specific  defect  concentrations  if  one  can  determine  the  absorption  cross 
sections  associated  with  the  different  defect  bands.  This  determination  was 
done  by  comparing  the  integrated  absorption  of  four  different  defect  bands 
with  the  measured  carrier  densities  for  a  large  number  of  samples  of 
different  Si  concentrations  (see  Section  II-B).  The  absorption  cross  sections 
were  then  used  to  measure  the  site  transfer  of  Si  from  Ga  to  As  sites  as  a 
result  of  Li  saturation  diffusions  at  elevated  temperatures  (see  Section  II-B). 
The  infrared,  electrical,  microstructural  and  mechanical  measurements  as 
well  as  the  analysis  have  been  completed.  In  brief  the  principal  results  are: 

(1)  Detailed  analysis  of  a  large  number  of  GaAs:Si  having  different 
[Si]  enabled  us  to  measure  the  absorption  cross  sections  of  the 
infrared  bands  attributable  to  4  different  Si  defects  Including 
Siga  donors  and  SiAs  acceptors. 

(2)  When  GaAs:Si  is  saturated  with  Li  by  diffusion  at  950°C,  Si  transfer 
from  Ga  to  As  sites  is  observed  and  the  degree  of  site  transfer 

Is  a  function  of  [SI]. 

(3)  Heavily-doped  GaAs:Si  have  Si-rich  faulted  loops  produced  by  Li 
diffusion  at  750°C  and  this  excess  Si  transfers  to  As  sites  If  Li 
is  diffused  at  950°C. 
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(4)  By  using  the  LVM  measurements  we  are  able  to  determine  the  changes 
in  [Siga].  [SiA$],  [Siga-SiAs]  an<*  tw0  unknown  but  Si-related 
defect  concentrations  as  a  result  of  annealing  at  400°C,  500°C, 
700°C,  800°C,  900°C  and  1000°C  ([  ]  =  concentration,  =  M  atom 
on  A  site).  Detailed  isothermal  anneal  curves  for  each  defect  were 
obtained  for  the  first  three  temperatures  for  three  different 
total  Si  concentrations. 

(5)  After  annealing  at  400°C  the  decrease  in  carrier  density  (nc) 
the  absence  of  observable  microstructure,  the  increase  In  yield 
stress  can  all  be  quantitatively  explained  on  the  basis  of  the 
decrease  in  observed  [Sigg]  from  the  LVM  measurements.  No  other  Si 
defect  concentration  changes  thus  making  a  Sigg  =  Vga  -*•  Sigg  -  V^a 
pair  formation  an  attractive  mechanism.  The  possible  SiGg  -  Vfts 
pair  defects  would  have  point  group  C3v  which  is  too  high  a  symmetry 
to  explain  the  yield  stress  results. 

(6)  For  samples  annealed  at  temperatures  above  400°C  the  changes  in 

ng  require  the  presence  of  a  new  and  as  yet  unidentified  acceptor. 
The  concentration  of  this  acceptor  increases  with  increasing  [Si] 
and  the  temperature  dependence  indicates  this  acceptor  concentration 
is  essentially  zero  at  400°C  and  1000°C  and  peaks  near  700°C. 

(7)  Microstructural  measurements  indicate  that  the  new  acceptor  is 
not  related  to  the  presence  of  dislocation  loops.  Extrinsic  loops 
are  observed  for  anneal  temperatures  near  700°C  and  are  essentially 
gone  at  1000°C.  The  loops  are  found  to  be  Si-rich. 


H-A-2 


(8)  At  the  higher  anneal  temperatures  the  changes  In  yield  stress 
still  qualitatively  correlates  with  the  changes  in  [SiGa]  which 

can  be  attributed  to  Siga-Vga  pair  formation  and  not  with  the  changes 
in  other  measured  Si  defect  concentrations.  This  indicates  the 
yield  stress  annealing  effect  is  still  controlled  by  the 
[SiGa-VGa].  To  our  knowledge  this  is  the  first  system  for  which 
this  type  of  data  has  been  available. 

(9)  With  very  long  anneal  times,  e.g.  ^1200h  or  more,  changes  in 
nearest  neighbor  pair  concentrations  [SiGa-SiAs]  are  observed  but 
only  for  a  very  limited  temperature  range  -  approximately  500-600°C. 
Above  and  below  these  temperatures  the  pair  concentration  is 
constant  in  time. 

In  this  next  section  we  have  given  detailed  accounts  of  the  research 
work.  Where  the  work  has  been  published  we  have  included  a  copy  of  the 
paper  as  the  account  of  that  portion  of  the  study.  In  this  way  the  final 
report  is  a  complete  account  of  all  of  the  research  done  under  this  grant. 
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B.  REVIEW  OF  RESEARCH  PROGRAM 

Infrared  absorption  bands  induced  by  Si-related  defects  in  GaAs:  Absorption 
cross  sections 

B.T.  Chen,  V.  Rana.*>andW.G.  Spitzer 

Departments  of  Physics  and  Materials  Science.  University  of  Southern  California,  Los  Angeles,  California 
90007 

(Received  29  June  1979;  accepted  for  publication  30  October  1979) 

Silicon  is  an  amphoteric  impurity  in  GaAs  and  is  present  in  different  defect  configurations.  This 
work  reports  measurements  of  the  infrared  absorption  cross  sections  for  some  of  the  localized 
vibrational  mode  bands  related  to  the  Si  defects.  Heavily  doped  melt-grown  samples  are  used  to 
measure  the  free-carrier  density  n,  and  the  integrated  absorption  of  all  bands  attributed  to 
electrically  active  defects.  The  resulting  data  are  analyzed  by  using  a  least-squares-fitting 
procedure  employing  two  different  models.  A  set  of  absorption  cross  sections  is  given  for  each 
model.  These  results  differ  from  others  reported  recently.  Using  the  two  sets  of  cross  sections  to 
interpret  the  absorption  spectra  of p-type  liquid-phase  epitaxial  samples  seems  to  reaffirm  a 
previously  reported  interpretation  in  which  an  additional  Si-related  acceptor  defect  of  unknown 
nature  is  present. 

PACS  numbers:  63.20.Pw,  81.40.Tv,  78.50.Ge 


I.  INTRODUCTION 

The  localized  vibrational  modes,  hereafter  called  LVM, 
of  Si-doped  GaAs  and  the  associated  infrared  absorption 
have  been  subjects  of  a  number  of  studies.1'13  Silicon  is 
known  to  be  an  amphoteric  impurity  in  GaAs,  and  materials 
grown  from  a  near-stoichiometric  melt  are  always  n  type. 
Absorption  bands  have  been  observed  and  attributed  to  sev¬ 
eral  defects  involving  Si,’  and  the  absorption  cross  sections 
for  two  of  the  more  numerous  and  electrically  important 
defect  species  have  been  reported.4  The  various  absorption 
bands  have  been  used  to  study  a  number  of  properties  of 
GaAs  :  Si  including  the  annealing-induced  changes  in  the 
carrier  density  of  heavily  Si-doped  material,*-*  the  Fermi 
level  effect  on  the  Si-site  distribution,9  the  behavior  of  Si 
when  it  is  introduced  by  ion  implantation,10  the  pairing  of  Si 
with  other  impurities,1 11 13  the  p-type  behavior  of  some  liq¬ 
uid-phase  epitaxially  grown  GaAs  doped  with  Si,4-11-13  lat¬ 
tice  strain  effects  resulting  from  fast  neutron  irradiation,3 
and  the  segregation  coefficient  of  Si  in  GaAs.3 14 

The  present  work  is  an  investigation  of  the  absorption 
cross  sections  for  the  absorption  bands  induced  by  Si-related 
defects  in  Si-doped  GaAs.  In  the  course  of  some  studies  of 
low-temperature  annealing  effects  it  became  obvious  that 
the  absorption  obtained  here  for  a  large  number  of  samples 
were  not  in  reasonable  agreement  with  predictions  based  on 
the  most  recent  and  presumably  most  accurate  values  re¬ 
ported4  for  the  cross  sections.  Because  of  this  difficulty  and 
because  the  cross  sections  are  needed  in  a  number  of  the 
studies  mentioned  earlier,  it  is  important  to  establish  reliable 
values. 

II.  BACKGROUND 

Because  of  the  amphoteric  nature  of  Si  impurities  in 
GaAs,  there  has  been  considerable  interest  in  the  use  of  mea¬ 
surements  of  the  LVM  frequencies  obtained  from  the  in- 
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frared  absorption  for  the  identification  of  defect  species.  Op¬ 
tical  measurements1-3  4  have  been  used  with  electrical  data 
and  expected  defect  symmetry  properties  to  identify  the  Si 
defects  responsible  for  the  observed  absorption  bands.  The 
“Si*;.  defect,  “Si  substitutional  on  a  Ga  site,  is  a  donor  and 
is  responsible  for  a  band  at  a  frequency  of  384  cm'*,  “Si*.  is 
an  acceptor  and  accounts  for  a  band  at  399  cm’1 ,  and  28SiG,  - 
“SiA,  are  nearest-neighbor  pairs  which  are  presumably  elec¬ 
trically  neutral,  account  for  bands  at  464, 393,  and  367  cm' 1 . 
All  frequencies  are  the  liquid-nitrogen-temperature  values. 
There  has  also  been  the  suggestion  that  a  band  near  369  cm'1 
seen  only  in  electron  irradiated  Si-doped  samples  is  related 
to  “Sio. -K0.,  where  V is  a  vacancy,  or  to  some  other  "Si- 
containing”  complex.2  All  of  the  bands  show  close  to  the 
expected  isotope  shift2  when  28Si  is  replaced  by  “Si.  The 
Sio.  -Si*,  bands  at  464  and  393  cm'1  show  mixed  band  fre¬ 
quencies  for  samples  containing  both  “Si  and  “Si,  i.e.,  bands 
*  '  J,Sio.  -“Su.  and  “S^,  -“Si*,  defects.  The  lack  of  such 
a  mixed  isotope  band  for  the  “Si  band  at  367  cm'1  has  led  to 
some  question  as  to  the  correctness  of  its  identification  as 
Si©. -Si a,  pair  band. 

In  order  to  make  the  infrared  measurements  which  gave 
the  above  frequency  assignments  it  was  necessary  first  to 
essentially  eliminate  the  free-carrier  absorption  which  effec¬ 
tively  obscures  the  LVM  absorption.  The  electrical  compen¬ 
sation  has  been  accomplished  by  saturation  diffusion  of  Li  oi 
Cu  at  an  elevated  temperature,  i.e.,  700  ST0S 1000  *C,  or 
by  penetrating  particle  irradiation,  i.e.,  1-2-MeV  electrons 
or  neutrons.  In  a  few  cases  partial  compensation  was 
achieved3 13  by  counterdoping  during  crystal  growth.  A  ma¬ 
jor  difficulty  with  all  of  the  above  procedures  is  that  the 
compensation  process  may  alter  the  distribution  of  Si  among 
the  various  defects  and  even  create  Si-related  defects  which 
were  not  present  in  significant  concentration  in  the  prooom- 
pensated  state.  Ample  evidence  exists  for  such  changes  with 
some  of  the  compensation  procedures.  Saturation  diffusion 
of  Li  or  Cu  results  in  the  formation1  of  Si^-Lio,  or  Sio,  ~ 
Cuq,  second-neighbor  pairs.  For  each  of  these  defects,  one 
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observes  a  new  set  of  LVM  bands.  Moreover,  it  is  known16 
that  Li  diffusion  at  TD  >900  "C  can  cause  site  transfer  of 
substitutional  Si  from  Ga  to  As  sublattice  sites.  Also  if  GaAs 
has  a  concentration  of  Si  in  excess  of  [Si]~2  X  10'*  cm  3, 
thermal  aging  at  temperatures  greater  than  7~400  *C  can 
cause  changes  in  the  free-carrier  concentration6  7  which 
have  been  shown  to  be  related  to  changes  in  Si-defect  concen¬ 
trations.  Thus  Li  and  Cu  diffusions  are  done  at  temperatures 
and  for  times  for  which  major  aging  effects  have  also  been 
observed.  The  compensation  by  electron  irradiation  is  fre¬ 
quently  done  with  the  sample  maintained  near  liquid-nitro¬ 
gen  temperature  during  the  irradiation  to  prevent  thermal 
annealing.  After  irradiation,  the  sample  is  allowed  to  warm 
to  room  temperature.  There  is  the  concern  that  mobile  na¬ 
tive  defects  introduced  by  the  irradiation  might  pair  with 
some  of  the  Si  defects  and  thus  alter  the  original  concentra¬ 
tions.  In  one  study4  the  irradiation  was  continued  to  fluences 
of  five  times  that  necessary  to  render  the  samples  transparent 
at  frequencies  where  there  is  little  lattice  or  LVM  absorp¬ 
tion.  No  significant  changes  were  observed  in  the  absorption 
at  the  LVM  frequencies.  However,  it  is  possible  that  the 
charge  states  of  the  native  defects  and  the  Si-defect  complex¬ 
es  and  hence  the  pairing  process  may  be  Fermi  level  depen¬ 
dent.  Nevertheless,  from  the  evidence  available,  the  e  -irra¬ 
diation  method  of  compensation  appears  to  be  the  least 
disturbing  to  the  precompensation  distribution  of  Si  defects, 
and  it  has  been  used  to  compensate  all  the  samples  in  this 
study  as  well  as  in  most  previous  studies. 

In  order  to  determine  the  concentration  of  each  of  the 
defects  which  induces  LVM  bands  one  must  know  the  ab¬ 
sorption  per  center  for  each  of  the  bands;  i.e., 

/b.n<i  oLvm  dv  /[defect],  wherea  is  the  LVM  absorption  coef¬ 
ficient,  v  is  the  frequency,  and  [  ]  indicates  concentration.  To 
obtain  the  absorption  per  center  a  comparison  is  made  be¬ 
tween  the  measured  free-carrier  density  n ,  prior  to  compen¬ 
sation  and  the  absorption  bands  observed  after  compensa¬ 
tion.  This  comparison  is  made  by  using 

",  =  a(a,4  )JM  -  b  (a„4  )j»*  (1) 

where  nt  is  the  carrier  density,  ap  the  peak  LVM  absorption 
less  background,  4  the  band  width  at  half-maximum,  and  a 
and  b  densities  of  centers  per  unit  absorption  with  J  oLVm  dv 
replaced  by  (apA  ).  Comparisons  of  a  large  number  of  sam¬ 
ples  of  different  n,  and  [Si]  should  lead  to  values  of  the  con¬ 
stants  a  and  b  satisfying  this  relation.  The  a  and  b  values 
obtained  depend  upon  the  validity  of  several  assumptions: 
(1)  the  absorption  of  a  LVM  band  is  proportional  to  the 
density  of  centers;  (2)  the  [SiG.  -SiA,  ]  does  not  contribute  to 
the  free-carrier  density;  (3)  the  compensation  technique  does 
not  alter  the  [SiG,  ]  and  [Si  *,  ]  values;  and  (4)  there  is  no  other 
defect  which  has  donor  or  acceptor  character  and  is  of  con¬ 
centration  significant  to  [Sio,  ]  or  [Si*,  ].  It  is  also  important 
that  the  values  of  [Sio,  ]  and  [Si*,]  in  the  samples  studied  be 
both  large  for  at  least  some  cases  in  order  that  values  of  both 
a  and  b  obtained  be  meaningful. 

The  most  recent  and  complete  measurements  of  a  and  b 
are  those  reported  by  Laithwaite  and  Newman,4  hereafter  L- 
N,  who  measured  14  samples  of  Bridgman -grown  GaAs  :  Si 
having  n,  values  ranging  from  3  X 1017  to  6. 1 X 10'*  cm"5  and 


10  LPE-grown  Si-doped  GaAs  layers.  The  samples  were  all 
compensated  by  e'  irradiation  and  generally  very  good 
agreement  was  obtained  between  the  measured  carrier  densi¬ 
ty  and  that  deduced  from  Eq.  ( 1 ).  The  values  for  a  and  b  used 
to  get  this  agreement  were  a  =  5.29  X  1016  cm'1  and 
b  =  6.82  X  10l6cm''.  However,  our  attempts  to  use  these 
values  for  the  measurements  to  be  discussed  in  Sec.  Ill  were 
unsuccessful.  Factors  as  large  as  between  3  and  5  were  noted 
between  the  elect  cally  measured  n,  values  and  those  ob¬ 
tained  from  Eq.  (1). 

Inspection  of  the  L-N  data4  reveals  that  for  all  of  the 
Bridgman-grown  samples  the  (apA  )3W/(apA  )3M  S0.15. 
The  n,  values  obtained  from  Eq.  ( 1)  by  using  the  L-N  a  and  b 
values  frequently  differ  from  the  electrically  measured  val¬ 
ues  by  as  much  as  10%.  Thus,  b  is  subject  to  considerable 
uncertainty.  The  value  of  b  comes  primarily  from  the  appli¬ 
cation  of  Eq.  (1)  to  the  LPE  samples.  However,  there  is  con¬ 
siderable  question 11  as  to  the  validity  of  assumption  (4)  for 
such  samples.  While  some  of  the  differences  between  the  L- 
N  results,  and  those  to  be  given  here  may  be  attributed  to 
differences  between  samples,  experimental  techniques,  and 
instrumental  resolution,  the  quantitative  differences  are 
such  that  it  appears  unlikely  that  these  are  the  primary 
sources  of  the  difficulty.  For  these  reasons  Sec.  Ill  contains 
some  detail  concerning  instrumental  tests  of  the  accuracy 
and  reproducibility  of  the  absorption  measurements  which 
would  not  usually  be  presented. 

III.  EXPERIMENTAL  METHOD 

The  samples  used  here  were  taken  from  a  series  of  hori¬ 
zontal  Bridgman-grown  ingots.  Three  additional  samples 
were  kindly  supplied  by  Newman.  In  some  cases  the  ingots 
were  polycrystalline  but  with  very  large  single-crystal  sec¬ 
tions  having  typical  dimensions  of  centimeters.  The  samples 
were  always  taken  from  these  single-crystal  regions.  A  total 
of  2 1  samples  were  used  for  the  a  and  b  determinations.  The 
total  [Si]  ranged  from  —  1 X 1018  to  ~  1 X  1020  cm-’  as  esti¬ 
mated  from  the  dopant  added  during  growth  and  the  known 
segregation  coefficient14  of  Si  in  GaAs  and  also  from  spec- 
trochemical  analyses.  In  a  few  cases  the  [Si]  was  also  mea¬ 
sured  by  using  an  electron  probe  microanalyzer.  Most  sam¬ 
ples  used  for  the  a  and  b  determinations  were  annealed  for  1 
h  at  1200’C  and  quenched  to  room  temperature  prior  to 
making  any  measurements  or  proceeding  with  electrical 
compensation  of  the  samples  by  e'  irradiation.  The  annealing 
treatment  establishes  a  common  thermodynamic  back¬ 
ground  for  all  samples  and  avoids  some  of  the  known  and 
previously  mentioned  lower-temperature  annealing  effects 
for  both  ",  and  the  Si-defect  concentrations.  These  anneal¬ 
ing  effects  are  generally  present  because  of  the  cool-down 
cycles  employed  in  the  growth  procedure.  Generally  the 
high-temperatue-quench  procedure  has  almost  no  effect  on 
the  ",  for  material  with  [Si]  Six  10"  cm'3  and  substantially 
increases  n,  for  samples  with  [Si]  >b  1 X 10”  cm'3.  The  result¬ 
ing  free -electron  densities  ranged  from  n,  ~  l  X 1018  to 
1 X 10”  cm'3.  After  the  high-temperature-quench  treatment 
sufficient  material,  ~0.1  mm,  was  removed  from  the  sur¬ 
faces  that  infrared  reflectivity  measurements  of  the  free-car- 
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tier  plasma  frequency  indicated  that  further  removal  of  ma¬ 
terial  did  not  change  the  carrier  density,  i.e.,  we  were 
measuring  the  bulk  carrier  density. 

Transmission  electron  microscopy  measurements  of 
the  high- Tannealed  +  quench  samples  showed  no  stacking 
faults,  no  precipitates,  and  almost  no  dislocation  loops.  The 
samples  appeared  free  of  microstructure  except  for  an  occa¬ 
sional  dislocation  loop.  Samples  which  had  not  been  given 
the  treatment,  i.e.,  in  the  as-grown  condition,  were  frequent¬ 
ly  found  to  have  extrinsic  stacking  faults,  particularly  those 
samples  with  the  larger  [Si].17 

Prior  to  electrical  compensation  the  free-carrier  density 
of  each  sample  was  measured  by  determining  the  frequency 
of  the  plasma  minimum  in  the  infrared  reflectivity.  The  val¬ 
ues  from  this  method  have  been  compared  with  Hall  mea¬ 
surement  results  and  very  good  agreement  was  obtained.' 
After  the  n,  measurements  all  samples  were  e  irradiated 
with  1  4-MeV  electrons  until  they  were  infrared  transparent. 
The  fluences  ranged  from  approximately  5  X  10' V/cm2  to 
2x  10* V/cm2  depending  upon  the  initial  n,  value.  In  all 
cases  dose  rates  of  approximately  5  X  10'V/cm2  sec  were 
used.  The  samples  were  maintained  at  a  temperature  below 
~  140  ‘K  during  the  irradiation  with  the  temperature  being 
monitored  by  a  thermocouple  mounted  along  side  the  sam¬ 
ples.  Transmission  electron  microscopy  measurements  of 
samples  after  e~  irradiation  showed  no  microstructural 
change  as  a  result  of  the  irradiation,  i.e.,  essentially  free  of 
any  structure. 

Because  of  the  differences  between  the  absorption  re¬ 
sults  to  be  given  here  and  those  reported  in  the  earlier  mea¬ 
surements  of  L-N,4  the  procedures  used  here  are  described  in 
some  detail.  After  electrical  compensation  the  infrared  ab¬ 
sorption  coefficient  at  liquid-nitrogen  temperature  was  mea¬ 
sured  for  all  samples.  The  measurements  were  made  by  using 
a  single-beam  grating  infrared  spectrometer  employing  in¬ 
terference  filters  to  isolate  the  first  order  spectrum.  The  sam¬ 
ple  Dewar  was  mounted  at  a  focal  point  in  the  exit  beam  after 
the  monochromator  section.  The  transmission  was  mea¬ 
sured  by  using  the  sample  in-sample  out  method  with  the 
sample  covering  one  of  a  pair  of  matched  aperatures.  A  radi¬ 
ation  thermocouple  with  a  commerical  13-Hz  phase-sensi¬ 
tive  amplifier  system  was  used  for  the  detection.  The  spectral 
range  covered  was  350  <  v<500  cm-1  which  is  known  to  cov¬ 
er  the  frequencies  of  all  the  absorption  bands  of  interest  in 
this  study.  The  samples  measured  were  generally  about  0.3- 
0.8  mm  thick  and  had  a  small  thickness  wedge  of  —0.2* 
eliminate  transmission  interference  fringes.  The  minimum 
transmission  which  is  at  the  absorption  peak  of  the  strongest 
band  ranged  from  —2  to  — 10%.  The  absorption  coefficient 
was  calculated  from  T  =  /Mmple  //0  =  ( 1  —  Rj1 
X  exp(  —  ax)/[  I  —  R 2  exp(  —  2 ax)],  where  R  is  the 
reflectivity. 

A  series  of  instrumental  tests  were  performed  in  order 
to  ensure  the  spectrometer  was  operating  properly.  The  tests 
were  to  ensure  that  a  low  scattered  light  level  existed 
throughout  the  measured  spectral  range,  that  the  second- 
and  higher-order  beams  were  properly  filtered  out,  and  that 
the  detection  system  was  indeed  linear  in  its  response  to  the 
input  signal.  Briefly,  some  of  the  tests  were  as  follows: 


(1)  Room-temperature  transmission  measurements  of 
—  1-cm-thick  crystals  of  NaCI  and  KC1  were  made.  Com¬ 
parisons  of  the  measured  absorption  with  values  reported  in 
the  literature"  for  v  5  600  cm'1  were  very  good.  The  trans¬ 
mission  of  both  crystals  for  30051  v  5  450  cm'1  was  less  than 
1%  as  expected. 

(2)  A  measurement  was  made  of  the  room-temperature 
transmission  of  a  high-purity  GaAs  sample  in  the  spectral 
range  of  the  two-phonon  absorption  (600  S  340  cm'1). 
The  measured  absorption  coefficient  was  within  <  20%  of 
the  literature14  values  at  all  frequencies  and  within  10%  for 
almost  all  frequencies.  Similar  measurements  for  Ge  gave 
similar  results.20 

(3)  Measurements  of  the  Sic,  absorption  band  were 
done  for  a  variety  of  spectral  slits  widths  (0.62  to  1.1  cm'1) 
with  qualitatively  the  expected  results  in  the  measured 
linewidth  being  obtained. 

(4)  Most  of  the  samples  were  remeasured  after  an  inter¬ 
val  of  several  months  with  good  reproducibility  observed  in 
(a,  A  )  for  all  bands,  generally  within  10%. 

(5)  A  sample  was  remeasured  with  the  /„  beam  (sample 
out)  aperature  covered  with  a  high-purity  GaAs  sample  of 
the  same  thickness.  Both  samples  were  at  liquid-nitrogen 
temperature.  The  calculated  absorption  coefficients  of  the  Si 
bands  of  the  Si-doped  sample  were  essentially  identical  to 
those  calculated  when  the  sample  out  position  was  an  open 
aperature.  This  test  eliminated  a  possible  error  due  to  the 
temperature  dependence  of  the  detector  sensitivity. 

(6)  A  different  blazed  grating  was  used  in  our  mono- 


F1G  I.  Localised  vibrational  mode  absorption  spectrum  of 
GaAs :  Si(  1 200  *C/ 1  h  anneal  compensated  by  electron  irradiation  at 
77  *K).  —  •  designates  WA-Si-2- J,  [Si]=70x  10”  cmJ;  —  +  designates 
KM-16-3,  [SiJ=r2.4x  10”cm ’. 
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chromator  and  also  room-temperature  transmission  mea¬ 
surements  were  made  of  a  few  samples  by  using  another 
commerical  infrared  system  with  the  results  in  both  cases 
being  again  very  satisfactory. 

(7)  Several  samples  were  measured  with  each  sample 
having  two  thicknesses  which  differed  by  a  factor  of  2.  The 
resulting  absorption  coefficient  spectra  for  a  given  material 
was  in  good  agreement  for  the  two  thicknesses,  i.e.,  within 
15%  at  all  frequencies  and  almost  always  within  10%. 

IV.  RESULTS  AND  DISCUSSION 

A.  Analysis  I:  Sia,  and  SiA,  only  electrically  active 

defects 

The  absorption  of  the  e'  irradiated  samples  is  illustrated 
in  Fig.  1  for  two  samples  of  different  [SiJ.  The  results  of 
similar  measurements  for  all  samples  are  summarized  in  Ta¬ 
bles  I  and  II.  In  a  few  cases  the  [Si]  is  known  from  electron 
microprobe  measurements  and  is  included  in  Table  II.  The 
peak  absorption  coefficient  less  background  and  half-width 
for  each  of  the  principal  bands  are  given  for  all  measured 
samples.  As  mentioned  previously  the  samples  were  mea¬ 
sured  twice  with  several  months  between  measurements, 
and  in  all  cases  the  apS  for  each  band  agreed  within  15%  for 
the  two  measurements.  The  mean  half-widths  are  1.5,  2.2, 
and  2.2  cm'1  for  the  SiG,  (384  cm-1),  Si^  (399  cm'1),  and 
SiG,-SiA,  (393  cm'1)  bands,  respectively.  The  instrumental 
resolution  during  all  of  these  measurements  was  ~0.75 


cm'1.  The  measured  half- widths  for  the  Sic,  and  SiAt  bands 
are  quite  close  to  the  L-N  values4  of  1 .4  and  2.0  cm' '.  Figure 
2  shows  a  plot  of  (ap4  ),**/«,  versus  (a,d  )JW/n, ,  where  n, 
is  the  measured  carrier  density.  All  of  the  samples  of  Tables  1 
and  II  except  the  last  three  samples  are  included  in  this  plot. 
These  three  samples  are  excluded  because  of  the  large  uncer¬ 
tainty  of  the  carrier  density  of  the  Litronic  LiD  sample,  i.e., 
n,  =  (l-2)x  10'*  cm'3,  and  the  significant  amounts  of  other 
acceptors,  particularly  Mg,  found  by  the  spectrochemical 
analyses  for  the  two  KM-26  samples.  The  straight  line  in 
Fig.  2  is  a  least-squares  fit  which  yields  values  for  a 1  and  b 1  of 

a1  =  10.1x10“  cm'1, 

6' =  9.3X10“  cm1, 

which  gives  n'  =  10.1  X  10“(a,d  )3M-9.3X  10“(apd ),»». 
In  calculating  n't  in  this  way  it  is  assumed  that  Sic,  and  SiA, 
are  the  only  electrically  active  defects  and  that  they  are  en¬ 
tirely  responsible  for  the  preirradiation  value  for  n, .  These 
calculated  n'  values  are  also  given  in  Table  II.  In  almost  all 
cases  the  optical  value  is  within  25%  of  the  measured  electri¬ 
cal  value  of  n, .  The  relative  chi,  x  =  [2"_ ,  [«'(/)  —  nr- 
(measured)]J/iV  —  n  —  1]I,J,  value  for  li  —  18  samples  anu 
n  =  2  parameters  is  1 .37  X 101*  cm'3.  The  values  of  a  and  b 
given  here  differ  significantly  from  those  given  by  L-N.4 

We  have  conducted  some  measurements  in  collabora¬ 
tion  with  Newman  in  order  to  help  establish  the  correct  ab¬ 
solute  values  for  a  and  b.  In  a  recent  communication,  New- 


TABLE  I.  Integrated  absorption  strength  of  the  LVM  bands  and  frte-carrier  density  for  the  Bridgeman  melt-grown  GaAs  :  Si  samples. 
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16.0 

2.70 

43  2 

570 

200 

114.0 

GA529 

1* 

2.50 

26.0 

1.40 

35.8 

2.7 

200 

5.4 

0 

0 

1.0 

200 

20 

GA503 

I" 

1.15 

11.8 

1.40 

15.7 

0.8 

2.00 

1.6 

0 

0 

0 

Litronic 

1" 

1-2 

10.2 

1.40 

14.3 

0.8 

2.00 

1.6 

0 

0 

0 

LiD 

KM-26 

1.07 

16.0 

1.30 

20.8 

1.4 

260 

3.6 

0 

0 

0 

2"' 

1.00 

15.0 

1.35 

20.3 

1.5 

220 

3.3 

0 

0 

0 

Average 

1.47 

218 

204 

249 

216 

mmmmm 

"Sample*  annealed  at  1200  "C/ 1  h  +  quenched  to  room  temperature.  "Samples  are  as-grown. 

“Sample*  annealed  at  1 100  ‘C/I  5  min  +  quenched  to  room  temperature.  'Sample*  excluded  from  leaxt-aquarewBtting  procedure  (aee  teat). 
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Table  II.  Comparison  of  electrically  measured  it,  value  and  calculated  optical  n,  values  for  tbe  samples  listed  in  Table  I. 


Ingot 

Sample 

it, (elect,  max.)  n\ 

(  X 10"  cm'*)  (X  10'*  cm5) 

(X 10"  car’) 

».lN 

(X 10"  cm'*) 

(SiKEPMA) 

(X 10"  cm'*) 

KM-16 

r 

5.42 

6.26 

5.79 

3.60 

2.7 

2* 

5.54 

5.97 

5.28 

3.38 

2.4  • 

3" 

5.42 

5.07 

5.37 

2.70 

2.4 

4‘ 

6.05 

5.82 

5.94 

3.07 

2.6 

5* 

7.20 

5.65 

6.22 

2.74 

3.1 

6* 

7.24 

5.33 

6.13 

2.47 

3.1 

WA-Si-2 

1* 

7.55 

6.64 

6.92 

3.30 

4.7 

2' 

7.S0 

7.77 

7.65 

3.86 

3* 

6.00 

649 

7.24 

2.98 

4* 

9.90 

8.08 

10.67 

2.85 

5.2 

5* 

9.90 

7.35 

11.28 

2.29 

KM-15 

1* 

S.65 

10.78 

9.66 

5.53 

WA-99 

1* 

"  S.90 

8.84 

8.20 

4.06 

2* 

9.40 

9.44 

8.14 

4.34 

KS 

1* 

8.20 

9.70 

8.94 

3.02 

2* 

5.10 

6.08 

5.33 

1.98 

GA529 

1* 

2.50 

3.10 

2.89 

.  1.99)  three  sample* 

GA503 

I* 

1.15 

1.43 

1.37 

0.94  >  *upp!ied  by 

Litronic  LiD 

!•■* 

1-2 

1.29 

1.23 

0.85  /  Dr.  Newman 

KM-26 

1** 

1.07 

1.76 

1.62 

1.12 

2*-* 

1.00 

1.74 

1.61 

l.ll 

•Samples  are  annealed  at  1200  "C/l  h  +  quenched  to  room  temperature.  ‘Samples  are  as-grown. 

•Samples  are  annealed  at  1 100  "C/15  min  +  quenched  to  room  ''Samples  excluded  from  least-squares-fitting  procedure  (tee  text), 

temperature. 


FIG.2.(a,d))M/n,  versus  (arA  for  all  the  samples  of  Tablet  I  and 

II.  The  line  is  a  least -squares  Bt  which  yields  valuta  of  a '  «•  10.1  X 10'*  cm'1 
and  h '  »  9.}  x  10'*  cm'1  for  it'  »s  '(a,  d  )m-k  ‘(a,  d  )?*  of  model  I  (see 
text). 
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man  indicated  that  a  spectrometer  inaccuracy  produced  a 
values  approximately  30%  too  large  for  aU  of  their  measure¬ 
ments.  This  error  should  not  influence  the  A  values  so  the  L- 
N4  values  for  a  and  b  should  be  raised  from  5.29  X 10'*  and 
6.82X  10“  cm'1  to  6.88 X 10“  and  8.87 X 10“  cm-1,  respec¬ 
tively.  The  values  of  nJ-'N  given  in  Table  II  are  calculated  by 
using  these  corrected  L-N  a  and  b  values  and  our  measured 
{apA  )  in  Eq.  (1).  Although  the  corrected  a  and  b  values  are 
closer  to  our  a 1  and  b 1  values,  the  differences  between  the 
values  of  n^~N  in  Table  II  and  the  n,  values  from  measure¬ 
ments  are  still  as  large  as  factor  of  2-4.  Clearly  the  30% 
correction  of  the  L-N  measurements  reduces  but  does  not 
resolve  the  discrepancies  discussed  above.  Itshould  be  point  - 
ed  out  that  the  present  a 1  and  b 1  values  were  derived  entirely 
from  crystals  grown  from  a  near-stoichiometric  melt  at  the 
GaAs  solidification  temperature.  No  samples  grown  at  low 
temperature  from  a  Ga  melt,  i.e.,  LPE  material,  were  used  in 
the  analysis.  This  latter  restriction  is  important  as  some 
studies4"*'12'13  have  indicated  that  the  defect  chemistry  of 
LPE  and  melt-grown  samples  may  be  quite  different. 

B.  Analysis  ||:  Two-donor  and  two  sccsptor  modal 

Although  the  principal  LVM  bands  found  in  GaAs :  Si 
samples  have  been  each  assigned  to  specified  Si-related  de¬ 
fects,  i.e.,  Sio.  (384  cm'1),  Si*,  (399  cm"'),  and  Si<j.-SiA, 
(464  cm*',  393  cm'1),  tbe  status  of  the  367-  and  369-cnr1 
bands  is  uncertain.  As  described  previously,  tbe  369-cm*1 
band  is  found  only  in  electron-irradiated  samples  of  Si- 
doped  GaAs.  It  has  been  suggested  that  this  band  is  related 
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FIG.  3.  |(a,4  W«,  +  |c"/o"|[(a,J  W»,]|  versus  (  |A  "/a"| 

X  [(0,4  ),«/»,  ]  +  |rf"/<i"|((a,4  )**/«,]!  for  all  the  samples  of  Tables  I 
and  II.  The  line  is  a  least-squares  fit  which  yields  values  of  a" 

=  +  10.0x10'*  cm'.  6  "=>=  — 12.7x10'*  cm',  c"  =  +25.3  x  10'*  cm', 
d"  =  -19.0X  10'*  cm’1  for  n"  =  a"(a,4  ),„  +  b  "(a,d  )w 
+  c"(a,4  )„,  +  d  "(a,4  )m,  of  model  II  (see  text). 


to  SiGl— Kg,  or  some  other  “Si-containing”  complex.  The 
367-cm  '  band  was  originally  identified  as  one  of  the  Sic  - 
Si  A,  pair  bands.  There  is  now  question  as  to  the  origin  of  this 
band  for  two  reasons.  The  first  was  mentioned  earlier  and  is 
the  lack  of  a  mixed  isotope  band  as  was  observed  for  each  of 
the  other  two  pair  bands.  The  second  is  that  it  has  been 


shown2’  that  the  (a, 4  )ibl  band  strength  can  decrease  upon 
prolonged  irradiations  but  the  other  pair  bands  do  not.  The 
ratio  of  (a,4  )MT  to  (a, 4  )393  for  the  samples  in  Tables  I  and 
II  shows  a  large  variation,  i.e.,  from  0.2  to  0.8.  Thus  it  ap¬ 
pears  that  the  367-cm' 1  band  is  not  one  of  the  Sic  -Si  a.  pair 
bands. 

The  (a,  4  ),67and(a,4  )J6,band  strengths  shown  in  Ta¬ 
ble  I  can  become  very  large  for  the  heavily  doped  samples. 
We  considered  the  effect  of  these  two  bands  with  the  pre¬ 
vious  384-cm'1  (Sic )  and  399-cm'1  (SiA.)  bands  by  using  all 
four  bands  in  the  least-squares-fitting  process.  One  can  as¬ 
sume:  n"  =  a"(a,4  )lt4  +  b  "(a,  4  )399  +  c"(cr,4 
+  d  "(t7,4  ),*,  in  a  four-dimensional  least-squares  fit  with¬ 
out  assuming  any  charge  status  for  these  four  bands.  The 
values  for  a",  b  ",  c",  and  d  "  can  be  determined  by  applyng 
the  four-dimensional  least-squares-fitting  procedure  and 
formula22  to  all  of  the  (a,  4  )  values  in  Table  I.  Figure  3 
shows  a  plot  of  [(a,  4  )iM/nt  +  |clVa"|[(a,4  )J67//f,]|  ver¬ 
sus  [  |6  "/a"|[(a,4  )399/n,  +  \d"/a"\[{arJ  W«,j j  for 
all  of  the  samples  of  Tables  I  and  II  when;  |  |  is  the  absolute 
value.  The  line  is  the  result  of  a  least-squares  fit  which  yields 
values  of 

a"=  +10.0X  lO^cm'1, 

6"  =  — 12.7X10"W, 
c"=  +25.3X  10,6cnr‘, 
d"=  — 19.0X  10l6cm'', 
which  gives 

n"  =  10.0 x  10l6(a,4  )3M  -12.7x  10,4( a, 4  )399 
+25.3 X  10,4(  a, 4  )*, -19.0X  10“(  a,4  )*». 

The  positive  and  negative  signs  for  the  values  of  a",  b  ",  r", 
and  d  "  represent  their  band  status  as  donor  and  acceptor, 

' respectively.  The  analysis  indicates  that  the  367-cm'1  band  is 
due  to  a  donor  defect  and  the  369-cm"1  band  to  an  acceptor. 
The  donor  and  acceptor  nature  of  the  defects  responsible  for 
the  384-  and  399-cm'1  bands,  respectively,  are  reaffirmed. 


TABLE  III.  LPE  layers  (GaAs  :  Si). 


Sample 

Growth  Temp. 

(°Q 

Sio. 

(a,  4  )>„ 
(cm1) 

Si*. 

(“,4)JW 

(cm*) 

367-cm"'  band' 
(a,  4  )„, 

(cm1) 

i. 

(elect,  mens.) 
(cm’) 

(cm*5) 

(cm’) 

|l 

750 

126.6 

80.7 

18.5 

4.8XI0'V 

5.3X  10”  n* 

7.1  X  10”  n 

/■  (2 

840 

338.4 

296.6 

32.6 

2.IX10'V 

6.6X  I0‘*  » 

9.4X  10”  « 

/♦ 

910 

165.4 

144.7 

24.4 

2.0X  10"; 

3.3X10'% 

4.3X  I0’% 

/  65 

Unknown 

288.5 

230.8 

Unknown1 

3x10'% 

7.7X10  '% 

Unknown' 

/  67 

900-807 

538.5 

442.3 

Unknown' 

3x  10'% 

13. 3x  10'*  n 

Unknown' 

1  68 

902-824 

961.5 

769.2 

Unknown ' 

3x  10'% 

25. 6x  10'*  n 

Unknown* 

\  79 

882-802 

111.5 

88  5 

Unknown' 

3x10  '% 

3  OX  10'*  n 

Unknown' 

]  80 

883-798 

96.2 

69.2 

Unknown' 

3xl0'% 

3.3x10 '% 

Unknown' 

//"  (  77  (No 

jSi) 

884-802 

0 

0 

Unknown' 

0 

0 

Unknown' 

1  89 

920-885 

69.2 

50.0 

Unknown' 

4.4X  10'% 

2.3x10'% 

Unknown' 

1  90 

919480 

215.4 

123.1 

Unknown' 

4.2X10'% 

10.3x10’% 

Unknown' 

1  128 

875 

153.9 

123.1 

Unknown' 

10”  P 

4.1x10'% 

Unknown' 

\  24 

900 

92.3 

61.5 

Unknown' 

10'*  n 

3.6X10'% 

Unknown' 

•Data  quoted  from  Kachan  n  al. '  *  'Data  quoted  from  L-N*  after  30*  a,  correction— *ee  text 

*Majority-camer  type.  *No.  369-cm*'  band  is  found  in  LPE  samples 

'Strengths  of  the  367<m'‘  bands  are  not  given  by  L-N.4 
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The  values  of  n "  calculated  in  this  way  are  also  given  in 
Table  II  with  n' .  In  almost  all  cases  the  optical  value  of  n"  is 
within  1 5%  of  the  measured  value  of  n, .  The  relative  chi 
value  for  N  =  18  samples  and  n  =  4  parameters  is 
0.88  X  10!8  cm'3  which  is  a  considerable  improvement  over 
the  first  model.  However,  this  result  may  be  simply  the  effect 
of  increasing  the  number  of  variables.  The  value  for  a"  is 
very  close  to  a  \  however,  b  "  is  about  30%  larger  than  b 
The  values  given  here  still  differ  significantly  from  the  modi¬ 
fied  values  of  L-N.  Although  this  four-band  least -squares- 
fitting  process  gives  a  better  fitting  result  than  the  previous 
model  with  only  two  bands,  its  validity  still  requires  further 
investigation,  especially  the  identification  of  the  367-  and 
369-cm' 1  bands.  It  should  also  be  pointed  out  that  our  calcu¬ 
lated  n\  and  n"  values  for  two  KM-26  samples  in  Tables  I 
and  II  are  much  higher  than  the  measured  nr  values.  This 
result  was  not  unexpected  as  significant  concentrations  of 
other  acceptors  were  found  in  these  samples  as  was  noted 
previously. 

C.  LPE  samples 

It  was  described  previously  that  heavily  Si-doped  GaAs 
crystals  grown  from  a  nearly  stoichiometric  melt  near  the 
melt  temperature  of  1238  *C  are  always  n  type.  However, 
GaAs  having  similar  Si  concentrations  but  grown  at  lower 
temperature,  frequently  between  750  *C  and  1000  ’C  and 
from  a  Ga-rich  solution  containing  dissolved  Si,  is  frequently 
p  type. 1313  The  conditions  which  determine  whether  a  solu¬ 
tion-grown  LPE  layer  will  bep  type  include  the  Si  concentra¬ 
tion,  growth  temperature,  and  crystallographic  orientation 
of  the  substrate  on  which  the  layer  is  grown.24  Some  major 
differences  in  the  characteristics  of  LPE  and  melt-grown 
GaAs  :  Si  samples  have  been  reported.1 1 31 7  Infrared  ab¬ 
sorption  results  of  these  two  types  of  materials  indicate  that 
only  bands  also  seen  in  melt-grown  samples  compensated  by 
e~  irradiation  are  observed  in  similarly  compensated  LPE 
samples.  However,  the  369-cm'1  band  which  is  present  in 
melt-grown  samples  is  absent  in  the  LPE  samples.  The  in¬ 
frared  absorption  results  fr„m  previous  studies  indicate  that 
(apA  )3M  >  (apA  ),„  even  though  the  LPE  layers  were  all  p 
type  before  irradiation.  If  these  defects  were  the  only  electri¬ 
cally  active  ones  and  the  cross  section  of  the  Si*.  band  is 
approximately  equal  to  that  of  SiG,  band,  the  samples 
should  have  been  n  type  prior  to  irradiation.  So  a  model,  in 
which  an  additional  acceptor,  which  is  not  observed  in  the 
infrared  absorption  spectra,  must  be  present  in  significant 
concentrations  was  proposed  in  some  of  the  previous  stud¬ 
ies.1213  Table  III  shows  integrated  absorption  data  for  the 
LPE  layers  quoted  from  the  Kachare  et  a/.13  and  the  L-N4 


studies  and  also  their  calculated  n\  and  /t"  by  using  the  two 
sets  of  cross  sections.  A  comparison  of  measured  n,  value 
and  calculated  n'  or  n"  value  for  each  p  type  or  even  n  type 
LPE  sample  seems  to  reafirtn  this  previously  reported  model 
in  which  an  additional  Si-related  acceptor  defect  of  un¬ 
known  nature  is  required. 
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ABSTRACT 

Infrared  absorption,  free  carrier  density,  and  transmission  electron  micros¬ 
copy  measurements  have  been  correlated  for  n-type  Si-doped  GaAs  samples 
with  1  x  1018  cm-8  <  [Si]  <5  x  101#  cm-3.  The  infrared  absorption  of  the 

localized  vibrational  modes  was  measured  after  the  samples  were  compen¬ 
sated  by  e~  irradiation  or  by  8Li  saturation  diffusion.  Comparisons  of  the  re¬ 
sults  of  the  two  compensation  methods  indicate  that:  (i)  the  Li-native  defect 
acceptor  complex  is  responsible  in  large  part  for  the  compensation  after  a 
750' C  Li  diffusion;  (ii)  the  degree  of  the  Si  site  transfer  from  Ga  to  As  lattice 
sites  as  a  result  of  the  Li  saturation  diffusion  at  950‘C  is  a  function  of  [Si] 
and  the  Li-native  defect  donor  is  probably  responsible  for  this  transfer; 
(iii)  samples  that  are  heavily  Si-doped  have  Si-rich  faulted  loops  produced 
by  a  heat-treatment  or  Li  diffusion  at  750‘C  and  this  excess  Si  transfers  to  As 
sites  if  Li  is  diffused  at  950’C;  (iv)  absorption  cross  sections  of  the  Li-native 
defect  donor  and  acceptor  bands  can  be  determined  which  account  for  the 
compensation  of  the  H-diffused  Si-doped  samples  and  also  fit  previously  pub¬ 
lished  data  for  Li-saturated  pure  GaAs.  There  is  also  the  suggestion  that,  at 
least  in  the  nondiffused,  nonannealed  more  heavily  Si-doped  material,  not 
all  of  the  Si  shows  up  in  the  localized  mode  measurements  f’en  though  there 
is  no  observable  reicrostructure. 


Results  are  presented  here  of  an  investigation  of  the 
behavior  of  Si  and  Li  when  they  are  simultaneously 
present  as  impurities  in  GaAs.  This  system  is  a  compli¬ 
cated  one  since  the  behavior  of  either  Si  or  Li  when 
individually  present  as  an  impurity  in  GaAs  is  com¬ 
plex.  However,  a  combination  of  optical,  electrical, 
and  microstructural  measurements  yields  considerable 
insight  into  the  nature  of  this  double-doped  system. 
In  this  study  one  observes  clearly  the  influence  of  Li 
saturation  diffusions  in  GaAs:  Si  at  various  tempera¬ 
tures  on  the  concentrations  of  both  Li  and  Si  defects 
for  a  variety  of  Si  concentrations.  The  presence  of  the 
Si  inhibits  the  formation  of  certain  Li-native  de¬ 
fect  complexes  while  the  presence  of  the  Li  can  cause 
the  transfer  of  Si  from  one  type  of  lattice  site  (Ga) 
to  another  (As).  The  measurements  of  compensated 
samples  permit  an  estimate  of  the  infrared  absorption 
cross  sections  of  the  prominent  localized  vibrational 
mode  (LVM)  bands  previously  attributed  to  the  two 
most  important  electrically  active  Li-native  defect 
complexes. 

The  Si  site  transfer  effect  in  Li-diffused  GaAsrSi 
has  been  observed  previously  but  was  not  studied  in 
any  detail  (1).  The  present  investigation  demonstrates 
that  the  site  redistribution  is  a  function  of  Si  concen- 
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tration.  It  also  shows  that  the  observed  LVM  absorp¬ 
tion  bands  produced  by  Si-related  defects  in  low  tem¬ 
perature  Li-diffused,  heavily  Si-doped  samples  do  not 
account  for  all  the  Si.  Transmission  electron  micro¬ 
scope  (TEM)  measurements  of  the  heavily  doped  sam¬ 
ples  after  Li  diffusion  reveal  some  microstructure 
which  could  account  for  the  unobserved  Si. 

Bockground 

Gallium  arsenide  doped  with  Si  has  been  of  partic¬ 
ular  interest  in  a  number  of  studies  because  of  the 
amphoteric  behavior  of  the  Si.  Silicon  substitution  on 
a  Ga  site,  Sice,  is  a  donor  and  on  an  As  site,  Si*,.  is 
an  acceptor.  The  distribution  of  Si  on  the  two  sub- 
I3ttices  is  known  to  vary  with  the  Si  concentration 
(2-4),  the  temperature  of  growth  (2-7),  the  tempera¬ 
ture  and  time  of  heat-treatment  (8-12),  the  nature 
and  concentration  of  other  impurities  (18-15),  and  the 
vacancy  concentrations  (1,8-11).  Electroluminescent. 
(18)  and  photoluminescent  (8-10)  characteristics, 
LVM  absorption  band  strengths  (1-7, 11-18),  and 
transport  properties  (17)  of  the  GaAs:  Si  are  dependent 
on  the  distribution  of  Si  on  the  two  lattice  sites.  Other 
Si  defects  such  as  Sic*-Si*.  nearest  neighbor  pain  and 
Sic,-native  defect  complexes  have  also  been  proposed 
and  studied. 

Lithium  is  an  electrically  active  impurity  in  semi- 
rnndurfnre  and  hue  been  the  subiert  nf  many  invettin- 


j 

I 


16b8 


JT.  Electrochem.  Soc.:  SOLID-STATE  SCIENCE  AND  TECHNOLOGY 


July  mo 


tiona  for  both  technological  and  scientific  reasons. 
While  lithium  behaves  as  an  interstitial  donor  in  the 
elemental  semiconductors  Si  and  Ge,  it  is  self-compen¬ 
sating  in  GaAs.  Lithium  can  form  either  donor  or  ac¬ 
ceptor  defects  in  GaAs  and  can  be  used  electrically 
to  compensate  either  n-  or  p-type  material.  The  rela¬ 
tive  concentrations  of  the  different  Li  defects  are  con¬ 
trolled  by  the  initial  concentrations  of  electrically  ac¬ 
tive  impurities  and  by  the  tendency  for  self-compen¬ 
sation.  There  has  been  a  considerable  volume  of  litera¬ 
ture  (2, 18-29)  dealing  with  the  detailed  behavior  of 
Li  in  pure  and  doped  GaAs,  and  as  a  result  of  these 
studies  a  number  of  specific  defect  species  involving 
Li  either  with  native  defects  or  paired  with  other  im¬ 
purities  have  been  identified. 

The  earlier  studies  of  the  GaAs: Si  +  Li  system  arose 
because  of  the  Li  behavior  as  an  electrical  compensator. 
In  order  to  make  infrared  measurements  of  the  LVM 
absorption  bands  in  GaAs:  Si  it  is  necessary  to  reduce 
the  free  carrier  absorption.  The  free  carriers  are  pres¬ 
ent  because  heavily  doped  GaAs: Si  is  strongly  n-type 
when  grown  from  a  nearly  stoichiometric  melt  One  of 
the  methods  used  for  carrier  removal  was  compensa¬ 
tion  by  Li  saturation  diffusion  at  an  elevated  tem¬ 
perature.  For  reference  Table  I  summarizes  all  LVM 
bands  observed  in  GaAs:  Si  +  Li  as  deduced  from  a 
number  of  studies  (1-7, 11-14, 25,  29-31).  Table  I  gives 
the  peak  absorption  frequency  for  each  band  measured 
at  liquid  nitrogen  temperature  and  the  detect  respon¬ 
sible  for  the  band,  where  Xa  refers  to  the  species  X 
on  the  A  sites.  According  to  the  theory  of  localized 
vibrational  modes  the  integrated  absorption  for  each 
band  is  proportional  to  the  volume  concentration  of 
the  defect  species  responsible  for  the  band.  We  refer  to 
the  constant  of  proportionality  as  the  absorption  cross 
section;  Jbud»  de/CXx],  where  a  is  the  absorption  co¬ 
efficient  and  [  ]  means  concentration.  The  *  =  384  and 
399  cm-1  bands  in  Table  I  are  due  to  Sic.  and  SiA* 
and  their  absorption  cross  sections  have  been  estab¬ 
lished  (3,4).  The  v  =  389  and  406  cm-1  bands  are  at¬ 
tributed  to  two  electrically  active  ‘Li-native  defect 
complexes,  and  the  absorption  cross  sections  for  these 
bands  are  given  in  a  later  section  of  this  work.  These 
latter  bands  do  not  involve  Si  as  they  occur  in  Li- 
diffused  undoped  GaAs  samples.  The  Si  concentration 
does,  however,  play  a  significant  role  in  determining 
the  concentrations  of  these  Li  defect  species.  In  Table  I 
the  band  at  „  =  367  cm-1  has  not  been  assigned  to  a 
specific  defect  It  was  previously  assigned  to  the  **Sic,- 
“Si..  nearest  neighbor  pair  defect,  however,  recent 
experimental  evidence  has  lead  to  the  conclusion  that 
this  assignment  is  not  correct  (•'.,  32). 

The  electrical  compensation  of  GaAs:  Si  which  is 
necessary  for  the  infrared  studies  has  been  accom¬ 
plished  by  saturation  diffusions  of  either  Li  or  Cu  at 
elevated  temperatures,  i.e.,  700*C  <  To  <  1000*C,  or 
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by  penetrating  particle  irradiation,  i.e,  1  to  2  MeV 
electrons  or  by  neutrons  (31).  In  a  few  cases  partial 
compensation  was  achieved  (6,  12)  by  counterdoping 
with  a  compensating  impurity  such  as  Zn  during  crys¬ 
tal  growth.  If  one  wants  to  study  the  distribution  of 
Si  among  the  various  defects  then  a  major  difficulty  of 
all  the  above  procedures  is  that  the  compensation 
process  may  alter  the  Si  distribution  and  even  create 
Si-related  defects  which  were  not  present  in  signifi¬ 
cant  concentration  prior  to  the  compensation  procedure. 
Ample  evidence  exists  for  such  changes.  Saturation 
diffusions  of  Li  or  Cu  result  in  the  formation  of  Sic- 
Lica  or  Sic«-Cuc>  second  neighbor  pairs,  and  Li  dif¬ 
fusion  at  Td  >  900‘C  can  cause  the  transfer  of  substi¬ 
tutional  Si  from  Ga  to  As  sites  (1).  Also  Li  and  Cu  dif¬ 
fusions  are  done  at  temperatures,  i.e.,  Td  >  400*C,  and 
times  for  which  major  aging  effects  have  been  ob¬ 
served  (11,  33),  and  these  effects  are  attributed  to 
changes  in  Si  defect  concentrations  as  well  as  the 
formation  of  new  defects.  The  compensation  by  elec¬ 
tron  irradiation  is  frequently  done  with  the  sample 
maintained  near  liquid  nitrogen  temperature  during 
the  irradiation  to  reduce  thermal  annealing.  Subse¬ 
quently  the  sample  is  allowed  to  warm  to  room  tem¬ 
perature.  There  is  the  concern  that  native  defects  intro¬ 
duced  by  the  irradiation  might  pair  with  some  of  the 
Si  defects  and  thus  alter  the  defect  concentrations 
(31).  In  one  study  (3)  the  irradiation  was  continued 
to  fiuences  approximately  five  times  that  necessary  to 
render  the  samples  transparent  at  frequencies  where 
there  is  little  lattice  or  LVM  absorption.  Essentially  no 
changes  were  observed  in  the  absorption  at  the  Sic.  or 
Sic.-SiAs  LVM  frequencies  as  a  result  of  this  extended 
irradiation  fluence.  Further  measurements  of  the  re¬ 
moval  rate  (31)  for  the  Six.  band  with  increases  in 
electron  fiuence  indicate  that  for  the  total  fiuences 
used  in  the  present  study  there  should  be  no  significant 
change  produced  in  this  band.  As  will  be  discussed  the 
present  measurements  are  in  accord  with  this  view. 
From  the  evidence  available  the  e"  -irradiation  method 
of  compensation  appears  to  be  the  least  disturbing  to 
the  precompensation  distribution  of  Si  defects  and  in 
the  present  work  it  is  assumed  that  the  Si  defect  con¬ 
centrations  are  the  same  as  the  precompensated  values. 

The  site  transfer  of  substitutional  Si  from  Ga  to  As 
sites  in  GaAs  has  been  inferred  in  several  studies.  Two 
reports  (8, 9)  describe  thermally  related  photolumi¬ 
nescence  changes  introduced  in  n-type  Si-doped  GaAs. 
In  one  case  the  change  was  by  heating  heavily  doped 
crystals  (n,  ~  4  x  10“  cm-3)  in  hydrogen  at  1050'C 
and  the  other  case  by  heating  lightly  doped  material 
(n.  ~  3  x  10”  cm-3)  in  evacuated  quartz  ampuls  at 
900*C.  These  two  studies  indicate  that  Si  atoms  may 
transfer  from  Ga  to  As  sites  via  an  As  vacancy  mecha¬ 
nism.  One  more  recent  photoluminescence  measurement 
(10)  of  lightly  doped  GaAs:  Si  (n.  ~  3.3  x  10”  cm-3) 
which  was  annealed  in  vacuum  at  a  temperature  of 
800*C  has  also  resulted  in  the  same  conclusion.  It  was 
also  reported  (13.  15)  that  for  heavily  double-doped 
crystals  of  GaAs:Si  +  Te  and  GaAs: Si  +  Se,  the  effect 
of  the  Te  or  Se  to  substantially  enhance  the  fraction 
of  the  Si  concentration,  [Si],  which  resides  on  As  sites 
and  thus  to  reduce  the  [Sic.]-  The  effect  of  redistribu¬ 
tion  of  Si  between  substitutional  sites  induced  by  the 
saturation  diffusion  of  Li  at  elevated  temperatures  has 
been  previously  observed  (1)  only  for  Si-doped  GaAs 
with  n,  ~  9.0  x  10”  cm-3  and  was  not  studied  in  any 
detail. 

In  the  present  study  a  comparison  is  made  of  the 
LVM  absorption  and  microstructure  for  identically 
doped  samples  compensated  either  by  *Li  diffusion 
or  by  irradiation  as  well  as  comparing  the  influence 
of  different  diffusion  temperature*  on  the  degree  of  Si 
site  transfer.  These  comparisons  are  made  by  using 
sets  of  identical  samples  for  the  different  treatments. 
Different  sets  of  samples  are  used  to  cover  a  rang*  of 
values  for  n*  and  [Si]. 
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Experimental  Merited 

The  temples  were  taken  from  several  horizontal 
Bridgman-grown  ingots  which  were  doped  with  dif¬ 
ferent  Si  concentrations.  Some  ingots  were  single  crys¬ 
tals  while  others  were  polycrystalline  but  with  very 
large  single  crystal  sections,  i.e.,  typical  dimensions 
of  centimeters.  The  samples  were  always  taken  from 
the  single  crystal  regions.  The  total  Si  concentration 
was  1  x  1018  cm-*  <  [Si]  <  5  x  101*  cm-3  as  esti¬ 
mated  from  the  dopant  added  during  crystal  growth 
by  using  the  known  segregation  coefficient  and  also 
from  spec trochemi cal  analyses.  In  a  few  cases  the  [Si] 
was  also  measured  by  using  an  electron  probe  micro¬ 
analyzer. 

All  samples  were  first  heat-treated  at  1200’C  for  1 
hr  and  then  quenched  to  room  temperature,  a  process 
hereafter  called  HT  +  Q.  The  HT  +  Q  treatment  es¬ 
tablishes  a  common  thermodynamic  background  for 
all  samples.  It  avoids  some  of  the  known  and  pre¬ 
viously  mentioned  lower  temperature  aging  effects 
for  both  the  free  carrier  density,  n,  and  the  Si  defect 
concentrations.  These  aging  effects  are  generally  pres¬ 
ent  in  the  as-grown  material  because  of  a  cool-down 
cycle  in  the  growth  process.  The  carrier  densities  after 
HT  +  Q  ranged  from  1.0  x  1018  cm-8  <  n,  <  8.0 

X  1018  cm-*.  The  HT  +  Q  procedure  has  essentially  no 
effect  on  the  n.  for  material  with  [Si]  <  1  x  1018 
cm-3  but  it  substantially  increases  n,  for  samples  with 
[Si]  >  1  x  101*  cm-3.  This  procedure  also  significantly 
improves  the  microstructural  quality  of  the  material, 
particularly  for  the  heavily  Si-doped  samples. 

After  the  HT  -f  Q  treatment  the  free  electron  den¬ 
sity  of  each  sample  was  measured  by  determining  the 
frequency  of  the  free  carrier  plasma  minimum  ob¬ 
served  in  the  infrared  reflectivity.  This  method  has 
been  compared  previously  (34)  with  Hall  measure¬ 
ment  results  and  very  good  agreement  was  obtained. 
Before  taking  the  measured  n,  as  the  bulk  value  suffi¬ 
cient  material  was  removed  from  the  surfaces,  —0.1 
mm,  that  the  infrared  reflectivity  measurements  of  the 
plasma  frequency  indicated  that  further  removal  of 
material  did  not  change  the  carrier  density. 

Several  samples  with  nearly  the  same  free  carrier 
density  were  cut  from  the  same  section  of  each  ingot 
and  were  divided  into  two  groups.  One  group  was 
compensated  by  e~  irradiation  while  the  other  group 
of  samples  was  saturation  diffused  with  8 Li  in  an  inert 
atmosphere,  some  at  Td  =  750*C  for  18  hr  and  others 
at  Td  =  950*  C  for  2  hr.  The  diffusions  were  done  with 
a  “sandwich"  arrangement  in  which  the  sample  of 
interest  was  placed  between  two  pure  GaAs  blocks 
which  had  previously  been  surface  alloyed  with  8Li 
on  the  outer  surfaces,  i.e.,  the  surfaces  not  in  contact 
with  the  sample  to  be  diffused.  This  procedure  sub¬ 
stantially  reduced  the  loss  of  material  which  occurs 
when  the  sample  of  interest  is  surface  alloyed  with  Li 
directly.  The  details  of  the  procedure  used  for  the  Li 
diffusions  have  been  reported  elsewhere  (22).  The 
samples  of  the  first  group  were  compensated  by  1.4 
MeV  e~  irradiation  while  at  temperatures  <  140’K 
and  with  electron  fluences  of  ~1  x  10“  e-/ cm*. 

The  infrared  absorption  coefficient  at  liquid  nitrogen 
temperature  of  the  Li-diffused  and  the  e~-irradiated 
samples  was  measured.  The  spectral  range  covered  was 
330  cm-1  <  p  <  900  cm-1  which  contains  the  frequen¬ 
cies  of  all  the  absorption  bands  of  interest  in  this 
study.  The  samples  were  generally  about  —0.3  mm 
thick  and  had  a  wedge  of  —0.2*  to  eliminate  multiple 
reflection  interference  fringes.  The  absorption  coeffi¬ 
cient  a  was  calculated  from  the  transmission  T  by  us¬ 
ing  the  usual  formula 

T  =  =  (1  —  R)*exp(— «r)/l  —  R*exp(—  2ox) 


where  R  is  the  reflectivity  and  x  is  the  sample  thick¬ 
ness. 

After  the  infrared  transmission  measurements  all 
Si-doped  samples,  as  well  as  some  pure  GaAs  samples 
which  had  been  saturation  diffused  with  *Li,  were 
used  to  prepare  specimens  for  transmission  electron 
microscope  measurements  (35).  The  TEM  specimens 
were  examined  by  a  Hitachi  HU-125C  microscope  oper¬ 
ated  at  125  kV.  Measurements  were  also  made  of  the 
HT  +  Q,  Si-doped  materials  before  Li  diffusion  as 
well  as  undoped  samples  and  they  were  free  of  micro- 
structure  except  for  an  occasional  dislocation  loop.  In 
particular  the  samples  showed  no  stacking  faults  or 
precipitates  while  the  Si-doped  samples  which  had  not 
been  given  the  HT  +  Q  treatment,  i.e.,  in  the  as-grown 
condition,  were  found  to  have  frequent  extrinsic  stack¬ 
ing  faults,  particularly  those  samples  having  the 
larger  [Si]  (36). 

Results  and  Discussion 

In  this  section  the  results  are  presented  in  the  fol¬ 
lowing  order: 

1.  The  changes  introduced  by  Li  diffusion  into  GaAs: 
Si  at  750‘C  are  examined  by  comparison  with  e” -ir¬ 
radiated  samples.  Some  of  the  irradiated  samples  will 
have  no  prior  annealing  while  others  were  thermally 
aged  at  700*C  prior  to  the  irradiation. 

2.  A  comparison  is  made  of  samples  which  have  been 
Li-diffused  at  two  different  temperatures  and  the  Li- 
induced  Si  transfer  from  Ga  to  As  sites  is  examined 
including  the  [Si]  dependence  of  the  transfer  process. 

3.  The  possible  role  of  the  known  Li-native  defect 
complexes  in  the  Si  site  transfer  process  is  explored. 

4.  The  LVM  absorption  measurements  of  the  Li  satu¬ 
ration  diffused  GaAs:  Si  samples  are  used  with  other 
data  to  estimate  the  absorption  cross  sections  of  the 
two  bands  arising  from  the  two  electrically  active  Li 
defect  complexes.  These  cross  sections  are  then  tested 
by  applying  them  to  previously  reported  LVM  mea¬ 
surements  of  Li-saturated  undoped  GaAs. 

5.  TEM  results  for  the  microstructure  of  the  Li- 
diffused  samples  are  discussed  and  correlated  to  the 
LVM  measurements  and  the  Si  site  transfer  observa¬ 
tions. 

Changes  in  defect  structure  and  concentration  re¬ 
sulting  from  *Li  saturation  diffusion  at  7 SO'C. — One  of 
the  primary  purposes  of  the  present  study  is  to  in¬ 
vestigate  the  site  transfer  of  Si  in  GaAs:  Si  which  re¬ 
sults  from  8Li  saturation  diffusions  at  an  elevated  tem¬ 
perature,  i.e.,  950*  C.  This  effect  is  demonstrated  by 
direct  comparisons  with  samples  diffused  at  a  lower 
temperature  of  750*C  where  the  transfer  process  is 
found  to  be  negligible.  However,  before  reviewing  the 
data  for  the  above  comparison  it  is  of  interest  to  in¬ 
quire  first  into  the  changes  introduced  by  Li  diffusion 
at  750*C  as  compared  to  the  prediffused  state.  The 
LVM  spectra  of  the  prediffused  state  is  taken  to  be 
that  of  the  e” -irradiated  samples.  As  mentioned  pre¬ 
viously,  the  heavily  doped  GaAs:  Si  undergoes  changes 
due  to  thermal  aging  at  temperatures  such  as  750'C 
and  therefore  the  effect  of  the  diffusion  on  Si  distribu¬ 
tion  could  be  the  result  of  the  750'C  diffusion  tem¬ 
perature  and/or  the  presence  ol  a  large  Li  concen¬ 
tration. 

To  demonstrate  the  changes  three  absorption  curves 
are  used,  one  for  a  sample  which  was  Li  saturated  at 
750*C,  a  second  e~  irradiated,  and  a  third  annealed 
near  730 'C  and  then  e~  irradiated.  These  comparisons 
allow  one  to  attribute  the  changer  as  due  primarily  to 
Li  saturation  or  due  to  thermal  annealing  at  the  dif¬ 
fusion  temperature. 

Figure  1  shows  the  LVM  absorption  bands  for  three 
samples  of  nearly  equivalent  [Si]  ~  5  x  10>(  cm”3 
and  n,vtx  101*  cm~3  after  the  HT  +  Q  treatment 
Sample  7-a,  see  curve  L,  was  compensated  by  *Li 
saturation  diffusion  at  790*C  and  sample  WA-2-3  was 
electron  irradiated,  see  curve  XL  The  data  ter  tbaee 


II-B2-3 
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384  cm" ' 


Fig.  1.  Infrared  absorption  spoctram  at  liquid  nitrogen  ttm- 
porataro  for  thrao  GaAs:Si  samples  haring  ntorly  cquiralont  [Si]  ~ 
5  X  10,#  cns_s  and  with  initial  n,  ~  8  X  10**  cm-3  oftsr 
1200*C/1  hr  +  qutnch  (sao  Table  II).  Curro  I:  sample  7-a,  750°C 
*U  diffused.  Curro  II:  tamplo  WA-2-3,  componsatod  by  c~  irradia¬ 
tion.  Curro  III:  sample  annealed  at  700rC  then  compensated  by  e~ 
irradiation. 

two  samples  including  the  integrated  absorption  (at,  A) 
of  the  LVM  are  among  those  of  Table  II,  where  op  is 
the  peak  LVM  absorption  coefficient  with  the  back¬ 
ground  removed  and  A  is  the  bandwidth  at  half-maxi¬ 
mum.  Curve  HI  of  Fig.  I  is  for  the  third  sample  which 
was  annealed  at  700 *C  before  being  compensated  by 
electron  irradiation.  The  (apA)ies  is  the  same  for  the 
three  samples  indicating  that  the  [Sia,]  is  independent 
of  the  compensation  method  or  thermal  aging  as  long 
as  the  temperatures  involved  are  ~750*C.  The  Sic- 
Sias  nearest  neighbor  pair  band  at  393  cm-1  is  the  same 
for  the  annealed  sample  and  the  diffused  sample  but 
they  both  are  approximately  twice  the  value  for  the 
nonannealed,  nondiffused  sample.  This  indicates  that 
there  is  a  substantial  increase  in  [Sics-Si*,]  for  the 
diffused  and  annealed  samples,  and  this  increase  is 
clearly  a  thermal  effect  rather  than  the  result  of  the 
presence  of  a  large  [Li].  The  (apA)ju  for  the  Sio  band 
of  the  annealed  sample  is  102.5  cm-*  which  is  below 


the  117.3  cm-3  of  the  nonannealed  sample,  but  sig¬ 
nificantly  larger  than  the  78.3  cm-3  of  the  diffused 
■ample.  The  decrease  of  [Sica]  in  the  diffused  sample 
therefore  appears  to  be  largely  a  combination  of  ther¬ 
mal  aging  and  the  formation  of  (Sic-Lice)  pairs 
which  give  rise  to  the  bands  at  374,  379,  and  405  cm*1, 
see  Table  I  and  Fig.  1-curve  L  Several  other  bands  are 
observed  in  the  absorption  spectra  of  Fig.  1.  The  weak 
bands  near  374  and  379  cm-1  in  the  nondiffused  sam- 
pies  are  due  to  the  LVM  of  the  **Sic  and  **SiCa  and 
they  superimpose  upon  two  of  the  MSic«-Lic,  bands, 
see  curve  L  The  367  and  369  cm-1  bands  are  poorly 
resolved  and  are  related  to  unknown  defects  involving 
Si;  see  Ref.  (4)  for  a  discussion  of  these  modes.  The 
389  cm-1  band  in  the  diffused  sample  is  due  to  a  Li- 
native  defect  complex  which  is  not  Si  related,  see 
Table  L 

These  results  necessitate  a  reinterpretation  of  some 
previously  published  data  in  which  the  samples  were 
not  given  the  HT  +  Q  treatment  prior  to  electron  ir¬ 
radiation  or  Tj  diffusion  at  700*C  (30).  This  prior  com¬ 
parison  between  Li  diffusion  and  e~  irradiation  was 
therefore  between  samples  which  had  been  thermally 
aged  in  an  unknown  temperature  profile  during  the 
growth  process.  Therefore,  the  effective  anneal  tem¬ 
perature  is  unknown  although  a  temperature  zone  of 
~800‘C  was  used  after  the  growth  zone. 

A  more  complete  comparison  of  the  nonannealed,  e~- 
irradiated,  and  750*C  *Li-diffused  cases  for  different 
[Si]  is  given  by  the  data  of  Table  IL  In  all  cases  the 
absorption  was  measured  at  liquid  nitrogen  tempera¬ 
ture.  The  data  of  Table  U  were  taken  from  samples 
such  as  those  in  Fig.  1.  The  data  for  the  irradiated 
samples  were  also  used  in  Ref.  (4)  where  they  are 
discussed  in  detail  From  Table  II  it  is  seen  that  the 
changes  observed  in  Fig.  1  and  the  conclusions  based 
on  those  changes  appear  to  apply  generally  with  larger 
changes  being  observed  at  higher  [Si].  In  particular,  it 
can  be  concluded  that  the  saturation  diffusion  of  Li  at 
750'C  produces  no  change  in  the  concentration  of  Six* 
acceptors,  significant  increases  in  the  concentration  of 
Sic»*SiAj  nearest  neighbor  pairs,  and  decreases  in  the 
concentration  of  Sic  donors  with  all  changes  being 
measured  relative  to  the  prediffused  state.  With  this 
information  we  may  now  examine  the  changes  which 
occur  when  the  saturation  temperature  is  increased. 

Si  site  transfer  induced  by  Li  saturation  diffusion. — 
Figure  2  and  curve  I  of  Fig.  3  show  the  absorption  co¬ 
efficient  of  three  samples  having  different  [Si]  which 
were  *Li  diffused  at  950*C.  Each  of  these  samples  is 
similar  to  one  of  those  diffused  with  »Li  at  750*C.  The 
results  of  750’  and  950°  C  saturation  diffusions  of  *Li 
are  summarized  quantitatively  in  Table  III  for  pairs 
of  initially  similar  samples.  Comparisons  indicate  that 
the  384  and  399  cm-1  bands  of  the  Sic  and  Six* 
respectively,  have  changed  their  relative  strengths 
while  the  other  Si-related  bands  remain  essentially 
unchanged.  The  384  cm*1  band  is  reduced  and  the  399 


Table  II.  Compariiani  of  LVM  retain  far  a~  irradiation  and  750*C  8U  diffusion 


Simple 

No. 

(cm-*) 

Compensation 

method 

(cm-*) 

(a».A)a« 

(cm-*) 

l>4l> 

(cm-*) 

(cm-*) 

(atS)m 

(cm-*) 

(cm-*) 

(tnl)jf 

(cm-*) 

KM-ZO-1 

t.er  x  io» 

e-  Irradiation 

too 

SO 

0 

0 

• 

1-a 

14  x  10“ 

•Li**  diffusion 

11.1 

0.7 

<•1 

— 

0 

-14 

-44 

SM-1S4 

M  x  10* 

70.1 

S14 

M 

14 

0.4 

to 

SJ  x  10“ 

■LI" 

00.1 

004 

104 

14 

114 

M 

KM-lt-4 

0.1  x  10“ 

01.0 

M.O 

0.0 

1.0 

0.4 

So 

0.7  x  10“ 

•LI" 

014 

00.0 

BJ 

04 

114 

so 

WA4M-S 

74  x  10“ 

r- 

U0.0 

074 

00.4 

14 

0.0 

to 

7.7  x  10“ 

BJ 

014 

40.0 

7.4 

174 

104 

WA-ffUU 

M  x  10“ 

#*■ 

1174 

07.4 

01.0 

7.0 

10.4 

7-0 

0.0  x  10“ 

u°° 

704 

004 

414 

04 

104 

114 

*  Itursrid 

••  «u  diffused  at  700°C  far  II  hr. 

t  Me  CEB’1  tend  only  tetrad  la  t-lmdUted  beevtly  dosed  (OB  pie. 
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Fig.  2.  Infrared  absorption  spectrum  at  liquid  nitrogen  tempera¬ 
ture  for  two  GaAs  ji  samples  haring  different  [Si]  and  initial  a(; 
both  *Li  diffused  at  950”C  for  2  hr.  Carre  I:  sample  6-b  (n,  ~ 
7.7  X  10»»  cm-*).  Curve  II:  sample  l-b  (n.  ~  1.0  X  1018  cm"8). 

cm-1  band  enhanced  for  the  higher  diffusion  tempera* 
ture  suggesting  site  transfer  of  Si  from  Ga  to  As  sites 
has  been  produced.  This  relative  decrease  of  the  384 
cm'1  band  after  the  higher  temperature  diffusion  de¬ 
creases  with  increasing  [Si]  while  the  relative  increase 
for  the  399  cm-1  band  increases.  The  SiG«-Six»  pair 
band  strengths  at  393  cm-*  are  almost  the  same 
for  two  diffusion  temperatures  for  each  set  of 
samples  indicating  that  the  large  annealing-induced 
increase  observed  at  750”  C  does  not  continue  with  in¬ 
creased  Li  diffusion  temperature.  There  is  no  consistent 
pattern  of  change  and  the  observed  changes  are  small 
between  the  two  diffusion  temperatures  for  all  other 
bands. 

One  can  examine  site  transfer  hypothesis  quantita¬ 
tively  if  the  absorption  cross  sections  of  the  Si  defect 
bands  are  known.  In  a  recent  study  (4)  the  cross  sec¬ 
tions  of  Sic*.  Si a>.  and  two  other  Si-related  defect 
bands  were  measured  for  a  number  of  melt-grown 
Bridgman,  HT  +  Q  samples  having  [Si]  and  n,  in  the 
range  of  those  used  here.  All  samples  were  compen¬ 
sated  by  e~  irradiation.  Assuming  the  *aSict  donor 
(384  cm-1  band).  **SiAt  acceptor  (399  cm-1  band),  and 
two  other  Si-related  defects  (367  and  369  cm-1  bands) 
to  be  electrically  active  defects  then  the  reciprocal  ab¬ 
sorption  cross  sections  were  determined  to  be  a  =  10.0 


406  cm'1 


Fig.  3.  lafrarud  absorption  sptetrum  at  liquid  nitrogen  tumpora- 
turo  for  four  GaAs:S<  samples  baring  an  identical  initial  »,  ~ 
5.3  X  1018  cm-*.  All  woro  8Li  diffused  and  in  ont  cut*  subst- 
qutnfljr  oufrfiffviod.  Curro  I:  sample  4-b,  950" C/2  hr  Hi  diffused. 
Curve  II:  sample  4-c,  900*C/3  br  Hi  diffused.  Curve  III:  sample 
4-a,  750”C/I8  hr  Hi  diffused.  Curve  IV:  sample  4-c’,  900”C/3  br 
•Li  diffused  +  700°  C/24  hr  8Li  outdiffused. 


X  10J«  crn-i,  b  =  12.7  x  10«  cm'1,  c  =  25.3  x  10" 
cm-1,  and  d  =  19.0  x  1018  cm-1,  where 

n,  =  a(opA)  jm  —  b(apA)sM  +  C(apA)se7  —  dlapAjsse 

=  [Sice]  —  [Sixs]  +  [Si]se7  —  [Si]j«» 

As  indicated  by  the  signs  the  charge  status  found  for 
367  and  369  cm-1  bands  are  donor  and  acceptor,  re¬ 
spectively.  For  all  of  the  Li-diffused  samples  shown  in 
Table  III  the  369  cm'1  band  is  not  observed  and  the 
367  cm-1  band  is  always  small  and  has  almost  no 
change  in  value  at  different  diffusion  temperatures. 
If  one  assumes  that  site  transfer  of  Sic*  to  Six*  is  the 
only  mechanism  responsible  for  the  changes  in  the  Si 


Tabts  III.  Comparisons  of  LVM  melts  of  750*  ond  950*C  Hi  diffusion  mtosurtmants 


Sumpla 
set  No. 


Sample 
No./Td  of 
•U  C C) 

(cm-») 

( aP  A ) 
(cm-*) 

(cm-*) 

(a^llxa 

(a*A)»r 
( cm-*/ 

(cm-*/ 

(cm-*) 

(cm-*) 

1 -4/730 

14  x  10“ 

in 

4.7 

—3 

0 

**14 

44 

0 

1-0/030 

LB  x  10“ 

4.4 

174 

—3 

0 

44 

414 

137.0 

S-u/730 

1J  x  1(P* 

344 

104 

4.3 

0 

4.0 

14.1 

4 

1-0/030 

14  x  10** 

M4 

30.0 

-44 

4 

3.0 

144 

■44 

S-x/730 

1.1  x  10** 

43.1 

111 

-4.0 

4 

4.4 

11.0 

4 

3-0/400 

1.1  x  10“ 

344 

414 

-4.0 

4 

44 

444 

33L7 

4-U/730 

44  x  10“ 

M.l 

304 

134 

14 

11.4 

33.0 

4 

44/M0 

>4  x  10“ 

444 

44.0 

13.1 

4 

13.4 

444 

3444 

lx/ 730 

4.7  x  10“ 

414 

334 

M4 

34 

114 

374 

0 

M/0H 

4.7  X  10“ 

47.0 

71.4 

33.4 

4 

114 

444 

4414 

4-0/730 

7.7  x  10“ 

414 

334 

434 

7.4 

17.1 

344 

• 

44/030 

7.7  x  10“ 

73.4 

110.7 

sac 

34 

114 

344 

817.1 

T-a/730 

44  x  10“ 

734 

M4 

414 

44 

14.3 

344 

4 

74/M0 

44  x  10“ 

M4 

43.0 

41.0 

44 

14.4 

344 

1444 

Siuusd  vahsaa 


II-BP-5 
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Tab!*  IV.  The  ratoa*  *#  0  calculated  from  Tobl*  III  os  a  fuactiaa 
of  initial  *, 


Sample 

■at  No. 

*•* 

(cm-*) 

l(o»4)<*4** 

(cm-*) 

(cm-*) 

— «(o*5)«4 

4(a*5)«i 

1 

1.0  x  10** 

—11.7 

417.1 

4  047 

7 

14  x  10** 

-»-5 

454 

4040 

S 

J.l  x  10** 

-145 

4749 

4041 

4 

54  X  10** 

-144 

4714 

4  0.45 

s 

5.7  X  10>* 

-144 

4154 

40.41 

s 

7.7  x  10“ 

-5.1 

4  574 

4  0.15 

7 

10  X  IP* 

+  9.0 

4744 

-047 

*  Miuund  nlusi. 

•*5(cwl )  m  < a»4.>  BtOSO'C  -  (0*4)  at  750'C. 


band  strengths  of  the  isolated  Si  species,  then  one  has 
6  [Sic*]  =  —  6  [Six*]  and  thus 

0  =  ((aph)sM/&(apA)3M  —  4*  b/o  =  ■+■  1-27 

Table  IV  gives  the  values  for  0  calculated  from  Table 
HI  and  Fig.  4  gives  0  as  a  function  of  n*  At  the  lower 
values  of  n.  most  of  the  change  is  indeed  accounted 
for  by  Si  site  transfer.  However,  as  tv,  and  [Si]  in¬ 
crease  the  value  of  0  decreases  and  the  change  becomes 
progressively  less  due  to  Sic.  -*  Six*.  Since  the  other 
Si  defect  concentrations  do  not  appear  to  change  sig¬ 
nificantly,  the  large  increases  in  (apA)SM  must  be  due 
to  transfer  of  Si  to  Six*  from  some  state  or  states  not 
detected  by  the  LVM  measurements  made  after  750*  C 
Li  diffusion.  For  the  samples  of  largest  [Si]  this  latter 
process  dominates  the  change.  From  the  TEM  mea¬ 
surements  mentioned  previously  it  is  known  that  there 
are  no  precipitates,  inclusions,  or  extrinsic  stacking 
faults  in  the  prediffused  samples.  Later,  TEM  mea¬ 
surements  will  be  presented  of  both  ".JO*  and  950*C 
Li-diffused  samples  which  show  that  Si-rich  loops 
form  at  750*C  and  disappear  at  950*C  for  the  more 
heavily  Si-doped  material.  Also  GaAs:Si  samples 
which  were  annealed  at  900 •  and  1000*C  and  then  «” 
irradiated  showed  no  evidence  for  the  Si  site  transfer 
process.  The  large  Li  donor  concentration  present 
under  the  conditions  of  the  high  temperature  satura¬ 
tion  diffusion  apparently  enhances  the  formation  of 
Six*  acceptors.  One  may  thus  conclude  that  site  trans¬ 
fer  of  Si  from  Ga  to  As  sites  is  indeed  observed  for  the 
high  temperature  Li  diffusion.  The  site  transfer  is  not 
simply  a  thermal  process,  i.«.,  the  presence  of  the  Li 
is  required.  Comparisons  between  low  and  high  tern- 


B*.  4.  A  plot  of  0  =r  —  8(apA)it*/&(*pA)}M  n.  loitlel  a*  far 
oM  of  the  iomqI**  of  TeMa  IV.  for  site  treat**  oalp  0  =1  JO. 


perature  Li-saturated  samples  indicate  that  at  the  low 
[Si]  almost  all  of  the  [Six*]  gain  for  the  high  tem¬ 
perature  case  is  the  result  of  Sic*  site  transfer  while  at 
the  highest  [Si]  almost  all  the  [Six*]  gain  results  from 
some  other  source  of  Si. 

Possible  rote  of  Li-native  defect  complexes  in  Si  site 
transfer. — From  the  preceding  section  we  concluded 
that  the  presence  of  the  Li  was  required  to  produce 
the  site  transfer  and  increase  in  [Six*].  It  is  of  interest 
to  inquire  whether  one  can  detail  in  any  way  how  the 
Li  is  instrumental  in  producing  this  effect  and  why 
the  higher  diffusion  temperature  is  required. 

The  strengths  of  the  Li-lattice  defect  complex  bands 
induced  by  *Li  diffusion  of  pure  GaAs  are  very  small 
for  a  low  diffusion  temperature  (21),  Le.,  at  Td  < 
750’C,  and  they  increase  monotonically  with  To-  The 
389  cm-1  (acceptor)  and  406  cm-1  (donor)  bands  are 
known  to  be  the  only  two  bands  induced  by  Li-lattice 
defect  complexes  which  are  electrically  active.  For 
pure  GaAs  the  ratio  (ap)« e/(«p)sa*  is  approximately  2 
at  higher  diffusion  temperatures,  To  >  800*C  (21,  22). 

Examination  of  all  absorption  curves  for  'Li-diffused, 
Si-doped  GaAs  samples  with  different  initial  n*  and 
[Si]  which  were  diffused  at  lower  temperatures,  i.e., 
TD  <  750’C,  indicates  that  the  389  cm-1  acceptor  band 
is  the  only  Li-lattice  defect  band  observed  with  sig¬ 
nificant  strength  and  it  increases  in  samples  with  higher 
n(  and  Si  values  as  shown  in  Table  V.  This  indicates 
that  the  growth  of  the  389  cm-1  Li-native  defect  ac¬ 
ceptor  band  may  be  a  significant  factor  in  the  com¬ 
pensation  of  Si-doped  GaAs  for  the  lower  temperature 
diffusions.  As  a  comparison  this  result  is  quite  differ¬ 
ent  from  that  of  Li-diffused  p-type  GaAs  for  a  similar 
diffusion  temperature.  In  this  latter  case,  no  Li-lattice 
defect  bands  are  observed.  In  particular  the  389  cm-1 
band  is  absent  The  compensation  U  completely  ac¬ 
complished  (28)  by  forming  LirMgc*  (or  -Znc„  or 
-Cdc*.  or  -Mnc«)  neutral  pairs. 

For  the  high  temperature  diffusion,  TD  ~  950*C, 
the  Li-lattice  defect  complex  bands  in  the  present 
samples  have  increased  substantially.  To  illustrate  this 
Fig.  3  shows  the  absorption  curves  for  four  GaAs:  Si 
samples  having  a  common  origin  and  identical  tv*  val¬ 
ues  but  *Li  diffused  or  outdiffused  at  different  tem¬ 
peratures.  Comparisons  between  these  curves  indicate 
that  the  saturation  concentration  of  Li  increases  due 
to  a  higher  diffusion  temperature,  the  389  and  406  cm'1 
Li-lattice  defect  bands  both  increase,  however,  the 
growth  rate  of  406  cm*1  Li-lattice  donor  band  is  much 
larger  than  that  of  the  389  cm~>  acceptor  band.  It  is 
this  donor  formation  which  likely  provides  the  mecha¬ 
nism  producing  Si  acceptor  formation  by  the  transfer 
of  Si  to  As  sites  from  other  defects.  When  Li  is  outdif¬ 
fused  at  a  lower  temperature  the  absorption  curve 
tends  to  return  to  the  curve  it  would  have  if  indif- 
fused  at  that  temperature;  see  curve  IV  of  Fig.  3.  In 
particular,  the  site  transfer  produced  in  a  900*C  dif¬ 
fusion  has  reversed  and  the  406  cm'1  band  has  disap¬ 
peared.  Furthermore,  a  950*C  *Li  rediffusion  after  a 
750*C  diffusion  produces  Si  site  transfer  with  similar 
LVM  curves  as  a  sample  diffused  initially  at  950*  C. 


Table  V.  Coa^arttaa  of  the  (o>)ut  value*  ter  all  750*C 
*U-4i  (fated  sample* 


Santa  No./ 
TaeTnj  CC) 

>u*  (em-*) 

(•>■»  (em-*) 

Pm  OaAj/750 

—44 

I-e/750 

14  x  10“ 

u 

t-e/790 

14  x  10“ 

54 

»*(1H 

5.1  x  10“ 

9.0 

40/750 

14  x  10“ 

104 

Se/710 

6.7  x  10“ 

114 

Se/750 

7.7  x  10“ 

14.4 

7-0/750 

(4  x  10“ 

144 

•  Maamred  vtkiaa. 
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ToUt  VI.  CaaforiiM  of  the  («p}«o*/(<>p)sM  ratio*  far  all  WC 
•U-diffvwMi  laxpila* 


Sample  No./ 

T»  e<  *LI  (*C) 

a.  * 

(cm-*) 

Max 

tern-1) 

<«»>«. 

(cm-1) 

>*»)«* 

<•>— 

Fui*  GaAa/SSTC 

0 

S3 

~04 

-1 

1-b/MO 

1.0  x  10“ 

ISA 

11.0 

S 

S-b/SSS 

l.a  x  10“ 

U 

170 

S 

S-b/MO 

1.1  x  10“ 

17J 

0S.0 

♦•b/sae 

M  x  10“ 

82.0 

105.1 

-e 

s-b/sse 

0.7  x  10“ 

21.0 

104.0 

S-b/SSS 

7.7  x  10“ 

14A 

SS.l 

-7 

7-b/OOS 

M  x  10“ 

14.0 

11.4 

<—4 

•  Maamiad  valuta. 


Table  VI  shows  that  the  (ap)*o«/(«b)sM  ratio  in  gen- 
feral  increases  with  increases  of  n.  for  samples  which 
ore  “Li  diffused  at  950  “C  for  2  hr.  This  observation  is 
consistent  with  the  result  that  in  more  heavily  doped 
samples  more  Si  transfers  to  As  sites  producing  more 
acceptors. 

The  conclusion  here  is  that  the  primary  Li  defect 
responsible  for  the  enhancement  of  the  Six*  concen¬ 
tration  is  the  Li  donor  responsible  for  the  LVM  band  at 
406  cm-1  (for  *Li).  The  increase  in  the  concentration 
of  this  Li  defect  with  temperature  as  measured  by  the 
LVM  absorption  is  the  reason  why  the  site  transfer 
requires  the  high  diffusion  temperature.  The  concen¬ 
tration  of  this  Li  defect  and  the  Sis*  both  increase 
with  the  total  [Si]. 

Absorption  cross  sections  of  Li-native  defect  com¬ 
plex  bimds. — From  the  discussion  in  the  previous  sec¬ 
tion  we  see  that  the  electrically  active  defects  in  the 
Li-saturated  samples  are  the  SiCm  donor,  the  Six* 
acceptor,  the  Si-related  donor  responsible  for  the  367 
cm-1  band,  the  Li  donor  responsible  for  the  406  cm-1 
band,  and  finally  the  Li  acceptor  giving  the  389  cm-1 
band.  One  can  now  use  this  information  to  obtain  the 
absorption  cross  sections  for  the  Li  bands.  For  all  of 
the  750"  and  950“C  “Li-diffused  samples  of  Table  III 
one  can  use  the  fact  that  the  samples  are  highly  com¬ 
pensated  after  diffusion  to  obtain 

7*«  =  0  =  —  b(apA)jM  +  c(apA)s97 

+  EiapAlus  +  F(apA)«oe 

where  a,  b,  and  c  have  been  previously  given.  If  the 
two  Li  defects  are  singly  charged  then  E  and  F  are 
the  reciprocal  cross  sections  for  the  Li  defect  bands. 
If  either  Li  defect  is  a  doubly  charged  donor  or  ac¬ 
ceptor  then  the  corresponding  coefficient  contains  a 
factor  of  2.  For  the  369  and  406  cm-1  bands  our  best 
estimated  values  of  A  are  2.4  ±  0.2  and  3.4  ±  0.2  cm-1, 
respectively.  Both  of  these  bands  overlap  other  bands 
which  frequently  made  determinations  of  A  difficult 
Figure  5  shows  a  plot  of  vs.  (opAljit/n-  Here 

n  =  —  [a(apA)sM  —  b(«pA)3M  +  c(«pA)3*7] 

=  EioyAljj*  +  F(  apA)4M 

where  all  samples  of  Table  HI  except  sample  2-b  are 
included.  Sample  2-b  was  omitted  because  e  was  orders 
of  magnitude  smaller  than  for  the  other  samples.  Also 
note  that  dUpA)***  is  omitted  as  that  band  does  not 
occur  in  Li -diffused  samples.  The  line  in  Fig.  5  is  a 
least  squares  fit  which  gives  values  for  E  and  F  of 

I  =  -  9.7  x  10“  cm~* 

F=  +  3.8X  10“  cm"* 

where  the  signs  indicate  the  389  cm~*  band  is  an  ac¬ 
ceptor  and  the  one  at  406  cm-1  is  a  donor  as  expected 
from  previous  assignments.  Application  of  these  values 
for  X  and  F  and  the  n,  equation  given  above  to  the 


10'“  c»"»  end  t  =  +18  x  10“  cm"*  (w*  text). 

samples  of  Table  III  yields  n,  values  for  all  cases  save 
one  which  are  <6%  of  the  sum  of  all  the  defect  con¬ 
centrations  used  in  the  n,  equation.  The  single  excep¬ 
tion  is  sample  2-a  where  the  figure  is  11%. 

When  7Li  is  used  as  the  diffusant  in  place  of  “Li 
then  the  Li  bands  have  an  isotope  shift  (22)  with  the 
389  cm-1  -*  364  cm'*  and  the  406  cm-*  -*  379  cm-*. 
For  the  same  Li  concentrations  there  is  also  a  change 
in  bp  with  (ap)ua/(ap)3*4  =  (ap)«oe/(ap)s7»  =  1.30. 
Since  the  A  shows  almost  no  change  with  Li  isotope  one 
has 

E(*Li)  =  1.30E  (“Li)  =  -12.6  x  10“  car* 

Fit Li)  -  i.30F(«Li)  =  +3.6  x  10“  cm"* 

As  a  test  of  the  validity  of  these  values  they  can  be 
applied  to  the  LVM  results  previously  published  (22) 
for  HJ  saturation  diffusion  in  pure  GaAs  for  700*C 
<  To  <  950*C.  Near  700*C  the  £  and  F  values  with  the 

previous  data  indicate  almost  exact  compensation  in 
agreement  with  experiment  However,  at  900*C  the 
estimated  n,  is  approximately  10%  of  the  total  Li 
donor  and  acceptor  charge  concentration  if  one  still 
assumes  that  the  389  (364)  and  408  (379)  cm"*  defects 
are  the  only  electrically  active  ones.  Experimentally  tho 
samples  are  highly  compensated  and  high  resistivity. 

One  can  also  use  the  above  values  of  £  and  F  to 
estimate  the  total  [Li]  in  the  saturated  pure  material. 
From  previous  work  (22)  it  is  known  that  the  389 
(364)  cm"'  acceptor  complex  has  2  Li  per  defect  and 
the  406  (  379)  cm"*  donor  defect  has  3  Li  Therefore 
the  total  [Li]  due  to  these  two  defect  centers  only  is 

N  =  2|E  (*Li)  |  («*A)  s*Jo  +  3|F(rLi)|(apA)t7«/m 

where  n  and  m  are  the  acceptor  and  donor  defect 
charge  states,  respectively.  Assuming  n  =  m  =  1  the 
calculated  values  for  N  are  the  points  in  Fig.  6  and 
they  can  be  compared  with  the  total  [Li]  which  is 
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l/kTptev'1) 

Fig.  i.  A.plet  of  InN  n.  1/kTo  for  the  7 Li  diffuse  pur*  GaA> 
Hinpl«*  in  Rtf.  (22).  Her*  N  -  2  f(7U)  |  («PA)SM  + 
3  |  f(7Li)  |  (opA)*?*  wfckh  is  Hi*  total  [Li]  du*  to  that  two  defect 
«*l*n  (auumiwf  •  =  m  =  1,  m*  text).  A  straight  lit*  ii  Hi* 
fatal  [Li]  measured  by  Fuller  and  Wolfsttin  [m  Ref.  (IS)]. 


given  by  the  straight  line  determined  by  the  flame 
spectrophotometric  analysis  measurements  of  Fuller 
and  Wolfstim  (18).  While  there  is  reasonable  agree¬ 
ment  between  the  calculated  points  and  the  line  there 
are  several  difficulties.  First,  there  is  no  independent 
evidence  that  n  =  m  =  1.  Second,  the  LVM  measure¬ 
ments  show  additional  bands  due  to  the  other  Li  defects 
and  they  are  not  included  in  the  calculated  N  values. 
Third,  the  determination  of  E  and  F  depends  on  the 
prior  knowledge  of  a,  b,  and  c  which  means  that  the 
value  of  N  involves  a  number  of  parameters  known 
with  only  limited  accuracy.  In  view  of  the  uncertain¬ 
ties  the  qualitative  agreement  in  Fig.  6  is  considered 
to  be  very  satisfactory.  Therefore  the  reciprocal  cross 
sections  tor  the  Li-donor  and  Li-acceptor  defects  are 
consistent  with  the  data  for  the  compensation  of  Li- 
saturated  GaAs:Si,  the  compensation  of  Li-saturated 
pure  GaAs  and  the  total  [Li]  in  the  Li-saturated  pure 
GaAs. 

Microstructures  of  ®Li-diJJuied  samples _ The  ex¬ 

tent  and  nature  of  the  microstructure  of  'Li -diffused 
samples  depend  both  on  the  starting  material,  particu¬ 
larly  the  [Si],  and  on  the  temperature  of  the  diffusion. 
It  will  help  to  clarify  the  results  if  they  are  discussed 
by  groups  of  samples  as  follows: 

(i)  As  previously  mentioned  nondiffused,  nonan¬ 
neal  ed,  HT  +  Q  samples  had  no  significant  observable 
microstructure  even  for  the  highest  [Si]  used  in  this 
study. 

(ii)  Pure  GaAs  samples  *Li  diffused  at  Td  750’, 
900*,  and  9S0‘C  were  measured  and  showed  Li  pre¬ 
cipitates  and  some  dislocations  which  tend  to  be 
pinned  by  the  precipitates.  Some  extrinsic  stacking 
faults  on  {111}  planes  which  may  be  related  to  the  Li 
metastable  phase  were  also  observed  in  the  750*C 
diffusion.  The  size  and  density  of  Li  precipitates  de¬ 
crease  with  decreases  in  diffusion  temperature.  In  the 
case  of  the  950*C  diffused  samples,  the  average  size  of 
Li  precipitates  is  approximately  1000A  and  their  den¬ 
sity  is  approximately  9.5  x  10“  particles/cm*.  These 
observations  are  in  agreement  with  those  reported 
previously  (24). 

(iii)  * Li-diffuied  samples  at  To  s  750*C  with  low 
[Si]  and  n,  <  3.0  x  10“  cm-*,  as  in  sets  No.  1  to  3  of 
Table  III,  have  a  microstructure  which  is  similar  for 
the  730'C  *Li -diffused  pure  GaAs  material  described  in 
case  (ii).  There  were  however,  no  extrinsic  stacking 
faults  observed.  Thermal  annealing  with  To  >  400*C  of 


nondiffused  samples  with  [Si]  ~  1.5  x  10“  cm-*,  i.e., 
set  No.  1,  resulted  in  no  observable  microstructure  as 
in  the  HT  +■  Q  state. 

(it?)  *Li  diffusion  at  To  —  750*C  for  the  samples 
with  high  [Si]  >  2  x  10“  cm-*  and  n,  >  5.0  x  10“ 
cm-*,  as  in  sets  No.  4  to  7  of  Table  m,  induces  a  high 
concentration  of  defects.  Samples  4-a  and  5-a  with 
[Si]  ~  2  x  10“  cm-*  show  a  microstructure  with  pri¬ 
marily  triangular-shaped  faulted  loops,  as  shown  in 
Fig.  7a.  The  loops  are  on  four  {111}  planes  of  a  tetra¬ 
hedron  and  have  the  displacement  vector  Rr  =  a/3 
[111].  To  the  best  of  our  knowledge  this  kind  of  defect 
has  not  been  reported  previously.  Similar  defects  are 
also  observed  in  the  material  with  similar  [Si]  ~  2  x 
10“  cm-*  when  annealed  at  750*C  without  *Li  diffu¬ 
sion,  and  their  density  is  similar  to  that  of  750*C  'Li- 
diffused  samples,  but  they  are  about  three  time  smaller 
in  linear  dimension.  By  using  two  different  and  inde¬ 
pendent  TEM  analyses  (37,  38)  we  have  established 
that  these  triangular  loops  are  extrinsic  faults  and  mi¬ 
croanalysis  with  a  scanning  electron  microscope  (SEM) 
shows  that  they  have  a  large  Si  content.  We  therefore 
conclude  that  these  triangular  loops  are  due  to  Si 
precipitation  from  GaAs: Si  solution,  however,  the 
mechanism  of  their  formation  is  unknown  as  is  the 
reason  for  the  observed  morphology.  The  microstruc¬ 
ture  of  samples  6-a  and  7 -a  which  have  the  highest 
[Si]  ~  5  x  10“  cm-*  show  primarily  circular  loops 
which  are  also  on  four  {111}  planes  of  a  tetrahedron 
as  shown  in  Fig.  7b.  The  TEM  analysis  (38)  has  estab¬ 
lished  that  these  circular  loops  are  also  of  the  extrinsic 
type.  As  seen  in  Fig.  7b  the  loops  here  are  much 
smaller  than  those  in  Fig.  7a  making  the  SEM  mea¬ 
surement  difficult.  However  it  seems  reasonable  to  as¬ 
sume  that  the  circular  loops  are  due  to  Si  precipitation 
from  GaAs:  Si  solution  as  in  the  case  of  the  triangular 
loops.  Typical  dimension  of  these  loops  is  —2500A  for 
the  sample  with  [Si]  ~2x  10“  cm-*  (Fig.  7a),  de¬ 
creasing  to  —450 A  for  the  sample  with  the  largest  [Si] 
~  5  x  10“  cm-3  (Fig.  7b).  The  loop  density  in  the 
former  case  is  —8.0  x  10“  cm-*  and  is  roughly  thirty- 
five  times  higher  or  —2.9  x  1014  cm-*  in  the  latter 
case.  This  indicates  that  the  amount  of  Si  precipitation 
through  the  formation  of  these  triangular  or  circular 
loops  increases  with  increasing  [Si]  in  the  750*C  'Li- 
diffused  samples.  Assuming  that  these  extrinsic  tri¬ 
angular  or  circular  loops  are  forming  completely  by  Si 
with  double  interstitial  layers  on  {111}  planes,  one 
can  estimate  [Si]  —  3.2  x  10“  cm-*  for  the  triangular 
loops  observed  in  sample  4-a  (Fig.  7a)  and  [Si]  ~  8.8 
X  10“  cm-*  for  the  circular  loops  observed  in  sample 
7-a  (Fig.  7b). 

(t?)  'Li-diffused  samples  at  To  a  950*C,  all  seta  In 
Table  QI,  all  had  similiar  microstructure,  i.e.,  Li  pre¬ 
cipitates  and  some  dislocations  which  tend  to  be 
pinned  by  the  precipitates.  This  is  similar  to  the  950*C 
•Li-diffused  pure  GaAs  material.  However,  the  size 
~700A  and  density  <  3.3  x  10“  particles/cm*  of  pre¬ 
cipitates  is  small  compared  to  that  of  the  pure  GaAs. 
One  previous  study  (18)  concluded  that  the  solid 
solubility  of  Li  in  doped  GaAs  is  frequently  larger  than 
that  in  pure  GaAs  which  can  account  for  the  reduction 
in  precipitates  in  the  Si-doped  material.  A  few  large 
dislocation  loops  are  also  observed  in  samples  with 
high  [Si]  >  2  x  10“  cm-*  as  in  sets  No.  4-7  In  Table 
III  The  character  of  these  large  loops  is  not  known, 
however,  they  would  account  for  only  a  very  small 
amount  of  [Si]  even  if  one  assumes  their  origin  is  due 
to  Si  precipitation.  As  previously  mentioned  *Li  re¬ 
diffusion  at  950*  C  of  sample  7-a  after  it  was  'Ll  dif¬ 
fused  at  750’C,  also  induces  site  transfer  and  shows  a 
similar  LVM  curve  as  that  of  sample  7-b  with  only 
950‘C  diffusion.  After  this  double  diffusion  the  micro- 
structure  does  not  show  the  high  density  of  circular 
loops  seen  after  the  first  lower  temperature  diffusion. 
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Fig.  7.  Bright  fi.ld  microgrophi  showing  tht  preienc.  of  a  high  concentration  of  loops  which  an  extrinsic  and  on  four  {111}  pianos 
of  a  tetrahedron  and  observed  after  750‘C  eLi  diffusion,  (a.  left!  Sample  4-0,  n,  ~  5.3  X  10“  cm-3  shows  triangular-shaped  faulted  loops 

-¥ 

with  density  '—8.0  X  1012  cm-3  and  the  displacement  vector  Kf  =  a/3  till].  Half  of  theta  triangular  faulted  loops  (y,  6)  do  not  show 
fringe  contrast  under  the  deflection  condition  being  used  here.  Foil  orientation  ~<I14>.  (b,  right)  Sample  7-a,  initial  n,  ~  8.0  X  1018 
cm-3  shows  circular  loops  with  density  ~Z9  X  I014  cm-3.  Foil  orientation  — <1 12>. 


The  resulting  final  microstructure  of  sample  7-a  is 
very  similar  to  that  of  sample  7-b. 

The  microstructural  results  described  above  can  be 
correlated  with  the  LVM  measurements  and  some  of 
the  results  discussed  in  the  four  preceding  sections.  For 
the  samples  of  lower  [Si]  no  Si-rich  loops  are  observed 
at  any  point  in  the  study  and  the  gain  in  (opA)  for  the 
399  cm-1  band  for  Six*  is  largely  the  result  of  Si  site 
transfer  from  Sic.  defects,  i.e.,  p  values  of  +0.97, 
+0.89,  and  +0.81  for  sets  No.  1-3  of  Table  IV  as  com¬ 
pared  to  the  predicted  value  of  +1.27.  On  the  other 
hand,  the  higher  [Si]  samples  show  the  Si-rich  loop., 
after  the  75 Q°C  diffusion,  but  essentially  none  after  a 
950°  C  diffusion.  These  same  materials  have  very  low  p 
values  indicating  that  the  increases  in  [SUi]  come 
largely  from  a  source  other  than  the  Sic  defect  con¬ 
centration  observed  at  750°C.  Since  no  other  Si-related 
LVM  bands  show  significant  changes  between  the  two 
temperatures  it  is  reasonable  to  inquire  into  the  possi¬ 
bility  that  the  Si-rich  loops  could  be  the  source  of  the 
high  temperature  increases  in  [Su,]-  By  using  the 
reciprocal  cross  sections  a  and  b  given  earlier,  one  can 
estimate  the  net  changes  in  [Si]  in  the  399  and  384 
cm-1  band  jointly  from  the  data  of  Table  IV.  For 
sample  4-b  the  increase  in  [Si]  is  2.7  x  10ls  cm-3  and 
5.7  x  1013  cm-3  for  sample  7-b.  These  values  are  quite 
close  to  the  estimated  [Si]  in  the  loops,  3.2  x  10<* 
cm-3  and  8,8  x  1013  cm-3,  respectively,  obtained  in 
the  previous  section.  Therefore  we  conclude  that  the 
Si  in  the  loops  formed  at  750*C  is  indeed  the  source  of 
the  additional  Si  observed  after  the  950*C  diffusion. 

There  remains  the  fundamental  question  of  the 
original  source  of  the  Si  which  appears  as  the  Si-ricb 
loops  after  a  750°  C  Li  diffusion  and  as  enhanced  [Si a.] 
after  950”C  Li  diffusion.  We  know  from  the  data  of 
Table  II  that  the  750°C  Li  diffusion  causes  a  reduction 
of  [Sic.],  •  reduction  of  both  the  367  and  369  cm-1 
Si-related  detects  and  a  substantial  increase  in  the 


[Sic'SiA.).  Also  the  Sic.-Lic.  bands  are  created  by 
the  diffusion.  The  [Si*,]  is  unaffected.  Since  the  cross 
sections  of  the  Sic.-SiA.  band  and  the  Sic.-Lic.  bands 
are  not  known  one  cannot  quantitatively  determine 
whether  the  total  [Si]  represented  by  the  bands  re¬ 
mains  constant.  It  is  conceivable  that  the  Si  in  the 
loops  and  the  enhanced  [Si**]  comes  from  net  de¬ 
creases  of  Si  in  the  Si  bands  observed  in  the  predif¬ 
fused  samples.  If  this  is  the  case  then  the  source  must 
be  either  the  Sic.  or  the  Si-related  defects  causing  the 
367  and  369  cm-1  bands  and  cannot  be  the  Sic.-SiA, 
defect.  Also  the  Si  in  the  loops  cannot  come  from  the 
Si  a.  defects. 

An  alternate  explanation  for  the  original  source  of 
the  Si  in  the  loops  or  enhanced  [Six,]  for  the  more 
heavily  doped  samples  is  the  possibility  that  some  of 
the  Si  dopant  does  not  give  rise  to  LVM  spectra  in  the 
prediffused  (e--irradiated)  samples.  An  indication 
of  this  possibility  comes  from  a  comparison  of  the  mea¬ 
sured  total  [Si]  with  that  calculated  from  LVM  spectra 
by  using  the  known  absorption  cross  sections.  It  is 
reasonable  to  estimate  Si  concentration  by  assuming 
[SiJc'.cui.irt  =;  [Sic.]  +  [Su,]  +  [Sic.-Su,]  +  [Si]3«r 

+  [Sljsso  ~  a(ap^)3S4  +  b(apA)  399  +  2a(aj>A)393  + 
c ( atp5) 387  +  d(apA)38»  for  e~-irradiated  samples  in 
Table  II  The  results  for  KM-16-3,4  and  WA-Si-2-2,3 
show  that  values  for  [Si]c.ieuu«<!  account  for  only 
60%-70%  ot  the  total  [Si]  values  measured  by  elec¬ 
tron  microprobe  analysis  which  were  2.3  X  10“  cm-3 
and  —5  X  10“  cm-3,  respectively.  In  addition,  New¬ 
man  (37)  and  co-workers  have  recently  studied  irradi¬ 
ated  samples  simultaneously  doped  with  Si  and  a 
group  VI  donor.  The  (apA)u.  band  of  Sic.  Is  about  the 
same  as  (0,4) jm  for  Si*,.  Prolonged  irradiation  (31) 
can  reduce  the  (erpA)s»#.  They  have  found  that  after 
such  an  irradiation  and  upon  annealing,  [S1a»]  recov¬ 
ers  and  (a^isM  can  overshoot  its  original  value  by  as 
much  as  a  factor  of  five  even  though  the  [Sic]  does 
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not  change.  Furthermore  in  a  most  recent  communi¬ 
cation,  they  have  also  observed  that  annealing  Si- 
doped  GaAs  samples  ([Si]~2  x  10u  cm"')  after  being 
irradiated  with  fast  neutrons,  the  strength  of  the  Six, 
band  also  increases  and  may  overshoot  its  original 
value  by  a  factor  of  three  but  without  any  measurable 
change  in  the  strength  of  other  LVM  bands  involving 
silicon.  They  have  also  concluded  that  there  is  some 
silicon  in  as-grown  crystals  which  does  not  give  rise 
to  LVM  absorption. 

Conclusions 

The  use  of  measurements  of  the  carrier  density,  the 
LVM  absorption,  and  the  microstructure  of  GaAs:  Si 
samples  has  provided  information  which  allows  one 
to  describe  in  considerable  detail  the  behavior  of  a 
number  of  defects  in  Li-saturated  material  Specifi¬ 
cally  the  conclusions  which  have  been  made  from  the 
present  study  are: 

L  The  effect  of  a  saturation  diffusion  of  Li  at  a  tem¬ 
perature  (750*0)  which  is  below  that  at  which  Si 
transfers  from  Ga  to  As  sites  results  in  significant 
changes  in  the  Si  defect  concentrations  as  compared  to 
the  prediffused  state.  These  changes  are:  (i)  a  de¬ 
crease  in  [Sic*];  (it)  an  increase  in  [Sic-Six,] ;  and 
(iii)  the  creation  of  Sic*-Lic.  defects.  There  is  no 
change  in  [Six*]-  The  changes  increase  with  increases 
in  [Si]. 

2.  For  samples  of  lower  [Si],  i.e„  [Si]  <  5  X  10'* 
cm-*  and  n,  <  3  x  10l*  cm-*,  saturation  diffusion  of 
Li  at  a  temperature  of  950*C  does  result  in  Sic,  -»  Si*, 
and  most  of  the  increase  in  [Six,]  is  attributable  to 
the  decrease  in  [Sic*]-  The  Li  defect  responsible  for 
inducing  the  Si  site  transfer  is  the  Li  donor  which 
has  an  LVM  band  at  406  cm-1  (for  *Li).  The  micro- 
structure  of  these  lower  [Si]  samples  includes  some  Li 
precipitates  and  some  dislocations  which  tend  to  be 
pinned  by  the  precipitates.  There  is  no  evidence  of  Si 
precipitation  or  extrinsic  faults. 

3.  For  samples  of  higher  [Si],  te.,  [Si]  >  2  x  101* 
cm-*  and  n.  >  5  x  101*  cm-*,  the  comparisons  of  LI 
saturation  diffusion  at  750*  and  950*C  are  more  com¬ 
plicated.  The  [Sic*]  generally  decreases  and  [Six*] 
increases  for  the  higher  temperature  case,  however. 
Sic*  -*  Six*  site  transfer  is  no  longer  the  dominant  proc¬ 
ess.  After  the  750*C  saturation,  TEM  measurements 
show  faulted  {111}  extrinsic  loops  which  are  Si  rich. 
These  loops  are  not  present  after  a  950*  C  saturation 
diffusion.  Semiquantitative  comparisons  indicate  that 
the  large  increase  in  [Six*]  at  the  higher  temperature 
could  come  from  the  Si  content  of  these  loops.  The 
question  of  the  origin  of  the  Si  in  the  loops  could  not 
be  answered  unambiguously  although  two  possible 
sources  were  identified  and  discussed. 

4.  By  using  the  fact  that  the  Li-saturated  sabiples 
are  electrically  compensated  and  using  the  known 
values  of  the  absorption  cross  sections  of  the  electri¬ 
cally  active  Si  defects,  the  absorption  cross  sections  of 
the  bands  for  the  principal  Li-donor  and  -acceptor 
defects  were  determined.  These  cross  sections  were 
used  with  the  knowledge  that  there  are  3  Li  per  donor 
defect  and  2  Li  per  acceptor  and  some  previously 
published  data  for  Li-saturated  undoped  GaAs  to  esti¬ 
mate  the  total  [Li]  as  a  function  of  saturation  tem¬ 
perature.  The  results  are  in  reasonable  agreement 
with  earlier  quantitative  analysis  data. 
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3.  SI-DEFECT  CONCENTRATIONS  IN  HEAVILY  Si-DOPED  GaAs: 

ANNEALING-INDUCED  CHANGES-II 
by 

R.T.  Chena^  and  W.G.  Spitzer 
Departments  and  Materials  Science  and  Physics 
University  of  Southern  California 

Los  Angeles,  California  90007 

ABSTRACT 

Isothermal  annealing  produces  changes  in  the  free  carrier  density, 
defect-induced  localized  vibrational  mode  (LVM)  infrared  absorption, 
microstructure  as  measured  by  transmission  electron  microscopy  (TEM), 
and  critical  resolved  shear  stress  of  heavily  Si-doped  GaAs.  The 
changes  have  been  measured  and  correlated  for  three  different  Si 
concentrations  for  several  annealing  temperatures.  The  measurements 
reveal  temperature  dependent  annealing- induced  changes  in  several 
specific  defect  concentrations.  The  observations  indicate  the  following 
behavior  for  two  Ingots  with  [Si]>2xl0  cm  :  (1)  When  the  anneal 
temperature,  T^=  400°C,  the  concentration  of  Siga  donors,  as  determined 
from  LVM  spectra  decreases  probably  due  to  the  generation  of  Vga  defects 
followed  by  the  formation  of  Sigg"^  pairs.  This  change  is  responsible  for 
observed  decreases  in  carrier  density  and  the  large  Increase  In  yield 
stress.  The  yield  stress  shows  a  dependence  of  the  form  o-o«  [S1ga-Vga]^4. 
(2)  When  TA=500°C,  the  LVM  spectra  Indicate  that  all  of  the  observed  Si 
defect  concentrations  change  .  The  decrease  in  [Si'Ga]  alone  cannot 
explain  the  decrease  in  carrier  density,  and  a  previous  suggestion  that  a 
new  acceptor  is  required  is  confirmed.  Both  the  LVM  measurements  and  the 
shear  stress  Indicate  that  only  a  small  fraction  of  the  [Sl^]  reduction 
Is  by  the  formation  of  $1Ga-VGa  pairs.  (3)  When  T^  *  700°C, 
a  new  acceptor  is  still  required  and  the  other  experimental 
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observations  at  =  500°C  are  also  still  seen  here.  There  is  a  large 

decrease  in  [Sigfl]  and  [Si^s]  observed  for  short  anneal  times  which 

coincides  with  the  formation  of  Si-rich  extrinsic  loops  and  the  loop 

area/vol  increaseswith  [Si].  (4)  When  >  700°C,  all  of  the  changes 

18  -3 

become  smaller  as  increases.  For  lower  [Si]~  1.5x10  cm  ,  no 
significant  annealing-induced  changes  are  observed  for  any  of  the 
given  above. 

a ^Present  address:  Electronic  Research  Center,  Rockwell 
International,  Thousand  Oaks,  CA  91360. 
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I.  Introduction 

There  have  been  a  number  of  studies^1 ”7^  of  annealing-induced 

effects  on  the  properties  of  heavily  Si-doped  GaAs.  Annealing  produces 

(3  41 

changes  In  several  properties  Including  the  free  carrier  density,'  ’  1  the 
(1-41 

photoluminescence,  '  the  defect-induced  infrared  localized  vibrational 
mode  (LVM)  absorption, the  mlcrostructural  characteristics  as 
determined  by  transmission  electron  microscopy^  (TEM),  and  the  critical 
resolved  shear  stress. The  studies  have  suggested  that  the  observed 
changes  in  these  various  properties  are  related.  The  present  work  is  an 
extension  of  some  of  these  prior  Investigations,  particularly  the  earlier 
work  of  Kung  and  Spitzer  in  Ref. 4,  hereafter  referred  to  as  Part  I.  The 
present  work  Is  a  detailed  and  systematic  study  of  the  annealing-induced 
changes  In  specific  defect  concentrations  and  the  influence  of  these 
changes  on  several  diverse  physical  properties  of  heavily  Si-doped  GaAs. 
Measurements  are  reported  here  of  the  carrier  density,  the  infrared  LVM 
absorption,  the  dislocation  loop  size  and  density,  and  the  yield  stress 
as  functions  of  annealing  for  samples  from  the  same  material.  The 
concentration  of  SI  Is  varied  by  about  two  orders  of  magnitude,  and  the 
annealing  temperatures  T^  cover  the  range  from  400°C  tn  1000°C.  The 
results  of  all  measurements  for  the  same  material  are  correlated  and  are 
Interpreted  on  the  basis  of  common  defect  models.  Specifically,  the 
present  work  extends  the  previous  studies  Including  those  of  Part  I  in 
the  following  ways: 

(1)  A  detailed  study  is  presented  which  covers  a  variety  of  annealing 
temperatures  and  times  as  well  as  three  different  Si  concentrations 
(hereafter  concentration  Is  Indicated  by  [  ]). 

(2)  Four  different  physical  measurements  are  made  for  essentially  all 
temperatures  and  times  by  using  samples  with  the  same  origin  and 

Identical  [SI].  The  common  origin  facilitates  detailed  comparisons 
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of  the  different  data.  The  previous  studies  did  not  use  material 
of  a  common  origin  and  properties  were  measured  for  only  a  very 
limited  set  of  annealing  conditions. 

(3)  The  study  is  a  semi -quantitative  one  in  that  specific  defect  and 
microstructure  concentrations  are  calculated.  The  previous 
studies  are  largely  qualitative. 

(4)  Annealing-induced  changes  in  the  carrier  densities,  LVM,  TEM,  and 
yield  stress  are  discussed  quantitatively  in  terms  of  the  changes 
of  the  specific  defect  concentrations.  It  is  found  that  the 
changes  in  defect  concentrations  indicated  by  the  LVM  spectra  can 
largely  explain  the  changes  in  the  other  measured  properties. 

The  subject  of  annealing-induced  changes  in  physical  properties 

of  GaAs  has  a  long  history  and  the  results  of  these  different  studies 

are  often  seemingly  contradictory.  Most  of  these  studies  have  Involved 

relatively  lightly  doped,  low  carrier  density  material  or  doped 

layers  on  the  surface  of  the  samples.  For  this  reason  one  might 

draw  the  conclusion  that  measurements  such  as  those  reported  here 

are  Indications  of  the  behavior  of  the  particular  samples  studied 

but  would  not  apply  to  randomly  chosen  samples.  It  should  be  realized 

that  the  three  principal  studies  of  the  annealing-induced  bulk 

changes  in  heavily  Si-doped  GaAs  have  been  quite  consistent,  see 

references  3  and  4  as  well  as  the  present  work.  These  studies  involve 

over  200  samples  taken  from  7  different  ingots  grown  by  two  different 

methods  in  two  different  laboratories.  Of  even  greater  significance 
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is  that  all  samples  having  a  high  Si  doping  level,  [S1]*2xl0  cm  , 
and  grown  from  a  nearly-stochiometric  melt  show  the  annealing  effects 
to  be  described.  As  will  be  discussed,  attempts  at  quantification 
and  modeling  these  changes  does  require  that  all  samples  start  from 

the  same  thermodynamic  equilibrium  state. 

II-B3-4 

/ 


The  mechanical  measurements  for  the  material  used  here  have  already 
been  reported,^  and  therefore  much  of  the  experimental  detail  concerning 
tnese  measurements  is  omitted  in  the  present  discussion. 


II.  Background 

The  amphoteric  nature  of  Si  in  GaAs  is  well  know^and  it  has  made 
this  impurity  of  particular  interest  in  a  number  of  studies.  It  is  known 


that  Si  substitutes  on  a  Ga  site,  Siga,  producing  a  shallow  donor 


state  and  on  an  As  site,  Si^s,  producing  an  acceptor  state. 


If  GaAs: Si 


is  grown  from  a  nearly  stochiometric  GaAs  melt  then  the  resulting  material 


is  always  n-type,  and  the  behavior  of  the  Si  is  reasonably  well  understood. 
At  low  Si  concentrations,  [Si]  <  1x10^®  cm"'*,  the  Si  preferentially 
occupies  the  Ga  sites, and  the  room  temperature  free  electron 


density ,n  .approximately  equals  the  [Si].  If  the  [Si]  increases  to 
larger  values  then  ng  does  not  increase  proportionally. )  The  Si 
forms  significant  concentrations  of  other  defects  such  as  SiAs  acceptors, 
^Ga'^As  nearest  neighbor  pairs, and  si-native  defect 
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complexes.^  J3-16)  Qeneral ly ,  for  as-grown  samples  with 

[Si]  ;t  10  cm"  ,  the  value  of  nfi  becomes  almost  independent  of  [Si] 

18  -3  (8) 

and  nfi  ^  5x10  cm"  .  '  The  exact  saturation  value  for  is 

dependent  upon  growth  method  and  thermal  history.  The  high  temperature 

solid  solubility  limit  for  Si  in  GaAs  is  -^0.5  at.  %  or  ^2x10^°  atoms  cm~3.^) 

Heat  treatments  at  elevated  temperatures  are  frequently  a  part  of 

the  procedure  of  processing  GaAs  materials  and  devices.  Examples  where 

the  heat  treatment  is  important  include  expitaxial  growth.  Impurity 

diffusion,  and  annealing  of  implanted  layers.  When  GaAs  crystals  are 

annealed  at  high  temperature,  i.e.  T^  >  700°C,  large  concentrations  of 

lattice  defects  are  frequently  introduced  near  the  surfaces.  There  have  been 

several  studies  which  reported  that  annealing  n-type  GaAs  increases  the 

acceptor  concentration  and  changes  the  photoluminescent  characteristics 

near  the  sample  surface.  These  changes  depend  on  the  As  pressure  during 

annealing.  Different  explanations  have  been  proposed  for  this  behavior 

including  contamination  by  Cu^18,19^  or  si^20’^,  the  creation  of  Ga  or 

As  vacancies the  transfer  of  Si  between  types  of  lattice  ; 

HI  20  30-331 

sites/  ’  ’  '  and  the  formation  of  complex  defects  involving  both 

impurity  atoms  and  lattice  defects.^ 30-34 ^  In  some  of  these  studies  the 
diffusion  coefficients  for  Ga  and  As  vacancies  have  been  estimated. 

Munoz  et  al  .^)  estimated  the  diffusion  coefficients  at  850°C  to  be 

-12  2  -1  -10  2  -1 
10  cm  sec  for  Ga  vacancies  and  10  cm  sec  for  As  vacancies. 

(271 

In  a  more  detailed  study  Chiang  and  Pearson'  ’  determined  the  temperature 
dependent  relations  for  the  diffusion  coefficients  which  give  values  at 
850°C  of  8x10"^ 3  and  1x10”^ 4  cm2  sec"^  for  Ga  and  As  vacancies,  respectively. 
Because  the  present  work  involved  thermal  anneals  at  temperatures  as  high 
as  1000°C  and  our  primary  interest  is  in  bulk  properties,  a  sufficient 
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amount  of  material  will  be  removed  from  the  surfaces  to  insure  that  the 
effects  of  the  changes  in  the  surface  layer  properties  on  the  measurements 
are  eliminated. 

As  previously  mentioned,  the  present  work  is  an  extension  of  some 
earlier  studies, and  for  the  sake  of  clarity  a  brief  review  of  these 
prior  studies  will  be  given  here.  In  particular  some  previous  LVM  studies  of 
GaAs:Si  are  considered  which  involve  annealing  effects  on  the  carrier 
density,  the  microstructure,  and  the  yield  stress.  For  further  details, 
including  more  extensive  reviews  of  the  early  studies,  the  reader  is 
referred  to  the  literature/3-*^ 

A.  Localized  Vibrational  Mode  (LVM)  Absorption 
Because  of  the  amphoteric  behavior  of  Si  in  GaAs  a  number  of  studies 
have  used  LVM  frequencies  obtained  from  infrared  absorption  measurements 
for  the  identification  of  defect  species  and  estimates  of  their  concentra¬ 
tions.  Infrared  measurements^1  *  ,13,35^  have  been  used  with  electrical 
data  and  expected  defect  symmetry  properties  to  identify  the  Si  defects 
responsible  for  the  observed  absorption  bands.  Several  of  the  observed 
LVM  bands  have  been  assigned  to  specific  point  defects.  The  bands  having 

liquid  nitrogen  temperature  peak  absorption  frequencies  at  384  and  399  cm  1 
28  28 

are  due  to  Siga  and  Si^$  defects,  respectively.  The 

-1  28  28 

bands  at  464  and  393  cm  are  assigned  to  Siga  -  Si^$,  where  the 

Ga  and  As  sites  are  nearest  neighbors.  These  pairs  are  presumed  to  be 
electrically  neutral.  Two  other  bands  located  at  367  and  369  cm  1 
have  been  attributed  to  two  other  unidentified  but  Si-related  defects. 

It  has  been  suggested^)  that  the  369  cm-1  band  which  is  seen  only  in 
electron-irradiated.  Si -doped  samples  and  the  367  cm-1  band  which  was 
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previously  incorrectly  assigned  to  the  Siga  -  SiA$  nearest 
neighbor  pair  defect,  are  both  caused  by  defects  involving  Si  paired 
with  native  lattice  imperfections.  The  bands  and  their  assignments  are 
given  in  Table  I. 

All  of  the  above  LVM  results  have  been  obtained  from  electrically 
compensated  material.  The  compensation  is  necessary  to  reduce  the 
density  and  hence  the  absorption  due  to  free  carriers  which  would  other¬ 
wise  effectively  obscure  the  LVM  absorption.  For  GaAs:Si  grown  from 
a  nearly  stochiometric  melt,  which  is  always  n-type,  saturation 
diffusion  of  Li  or  Cu  at  an  elevated  temperature,  i.e.  700°C  <  Tn  <  1000°C, 
and  penetrating  particle  irradiation,  i.e.  1  to  2  MeV  electrons,  have 
been  used  as  compensation  techniques.  If  however  one  wishes  to  study 
the  pre-compensation  distribution  of  Si  among  the  various  possible  defects 
then  one  difficulty  of  these  procedures  is  that  the  compensation  process 
may  alter  the  distribution  of  Si  defects  and  even  create  Si  related  defects  which 
were  not  present  in  significant  concentration  prior  to  the  compensation 
procedure.  Ample  evidence  exists  for  such  changesv  *  ’  '  where  the 

compensation  is  produced  by  saturation  diffusions  of  Li  or  Cu^and 
the  effect  of  this  compensation  process  on  the  Si  distribution  has  been 
previously  discussed. From  the  evidence  available^®’^  »37 ,38) 
e’-irradiation  method  of  compensation  appears  to  be  the  least  disturbing 
to  the  pre-compensation  distribution  of  Si  defects,  and  it  has  been  used 
to  compensate  all  the  samples  in  this  study  as  well  as  in  many  of  the 
previous  studies.  This  compensation  method  has  the  further  advantage 
of  reducing  the  number  of  defect  bands  compared  to  those  in  the  Li  or 
Cu  compensation  cases  and  thus  reducing  the  complexity  of  the  interpretation. 
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According  to  the  theory  of  localized  vibrational  modes  the 

integration  of  the  linear  absorption  coefficient  over  frequency  for 

each  band  is  proportional  to  the  volume  concentration  of  the  defect 

species  responsible  for  the  band.  We  refer  to  the  constant  of 

proportionality  as  the  absorption  cross  section;  /  adv./[Defect],  where 

band  ' 

a  is  the  LVM  absorption  coefficient  and  v  is  the  frequency.  In  the 

present  study  the  integrated  absorption,  /  adv,  is  replaced  by 

band 

(apA)band»  where  ap  is  the  peak  LVM  absorption  less  background  and 

a  is  the  bandwidth  in  cm'1  at  half-maximum. 

The  LVM  absorption  cross  section  for  each  of  the  Si  defect  bands 

must  be  known  if  the  bands  are  to  be  used  to  determine  the  changes  in 

defect  concentrations.  In  a  recent  study^11^  the  cross  sections  of 

Siga*  Si As,  and  the  two  Si-related  defect  bands  of  Table  I  were  measured 

for  a  number  of  samples  from  Bridgman,  melt-grown  ingots.  The  samples 

had  all  been  subjected  to  the  same  initial  treatment  of  annealing  for 

lh  at  1200°C  followed  by  quenching  to  room  temperature.  The  purpose  of 

this  treatment  was  to  remove  any  thermal  aging  effects  produced  by 

temperatures  in  the  400°C  to  900°C  range  during  the  growth  period. 

The  importance  of  this  procedure  will  be  discussed  later  in  the  present 

18  -3  20  -3 

work.  The  samples  had  concentrations  of  1x10  cm  [Si]  £1x10  cm 

18  -2  19-3  _ 

and  1x10  cm  ng£  1x10  cm  and  were  compensated  by  e  -irradiation. 

28  -1  28  -1 
The  Siga  donor  (384  cm  band),  Si^$  acceptor  (399  cm  band)  and  two 

other  Si  related  defects  (367  and  369  cm-1  bands)  were  assumed  to  be  the 

only  electrical ly-active  defects  prior  to  the  irradiation  and  thus,  they 

determine  the  carrier  density  after  the  high  temperature  plus  quench 

treatment.  With  this  assumption  the  reciprocal  absorption  cross  sections 

were  found  to  be  a^lO.OxlO1  ^cm’1 ,  b=12.7xl01^cm"1 ,  csZS.SxlO^cm  1 , 


and  d  =  19.0x10^  cm"^ ,  where 


ne  =  a^apA^384  "  b^apA^399  +  c^apA^367  "  d^apA^369 
-  [Si.  ]  +  [Si]1<;7  -  [Si],ca- 


The  electrical  behavior  of  the  defects  responsible  for  the  367  and  369  cm 
bands  are  donor  and  acceptor,  respectively.  These  values  of  the  constants 
will  be  used  in  the  quantitative  analysis  of  the  LVM  results  in  the 
present  work. 

B.  Carrier  Concentration 

Kung  and  Spitzer  (KS)  were  the  first  to  report^  the  annealing 

effects  in  heavily  Si -doped  GaAs.  The  effects  are  qualitatively  similar 

to  those  reported  earlier^”^3 Vor  GaAs  doped  with  Te  and  Se  impurities.  It  was 

observed  by  KS  that  the  n  could  be  decreased  by  as  much  as  a  factor  of  6 

e 

I  Q  O 

by  annealing  GaAs:Si  with  [Si]~  4x10  cm  at  temperatures  between 

400  and  900°C.  In  a  later  study  called  Part  I^KS  attempted  to  determine  the 

source  of  the  annealing-induced  decrease  in  ng  by  measuring  the  LVM 

infrared  absorption  and  low  temperature  photoluminescence  for  samples 

similar  to  those  that  were  employed  in  the  ng  investigation.  It  was  concluded  that 

the  decrease  in  ng  was  related  to  the  decrease  in  the  [SIq  ]  donors. 

After  annealing  at  400°C  the  loss  in  [Siga]  could  quantitatively  account 
for  the  decrease  in  ng  while  after  annealing  at  600  or  750°C,  the 
decrease  in  f Siga ]  was  insufficient  to  account  for  the  ng  change  and 
thus,  an  additional  unknown  acceptor  species  was  postulated.  After 
the  600  and  750°C  anneals,  the  growth  of  a  new  photoluminescence  band 

r- 

near  1.0  ev  was  attributed  to  this  new  acceptor  defect.  It  was  also 
observe  at  all  annealing  temperatures  that  the  concentrations  of  S1^$ 
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and  Siga  -  Si^s  pairs  as  determined  by  strengths  of  the  LVM  absorption 
bands,  remained  essentially  constant.  This  observation  indicated  that 
neither  site  transfer  nor  Si-Si  association  reactions  involving  these 
Si  defects  were  related  to  the  decrease  in  strength  of  the  Siga  band. 
Based  on  these  results,  KS  tentatively  suggested  that  the  observed 


decrease  in  the  strength  of  the  Siga  band  might  be  due  to  the  formation 
of  Siga  -  Vga  pairs,  where  represents  a  vacancy  on  the  Ga  site  which 


-is  a  second  neighbor  of  the  Siga 


This  suggestion  of  the  Siga  -  VGa 


defect  was  based  largely  on  earlier  studies  of  photoluminescence  spectra. 


It  is  clear  that  if  Vfia  is  an  acceptor  species  as  Is  generally 

assumed^then  pairing  Sifia  donors  with  pre-existing  acceptors 

will  not  result  in  a  decrease  In  n  .  The  decrease  in  n  must  result 

e  e 

from  the  generation  of  Vga  defects  from  internal  sources  and 

these  vacancies  can  then  pair  with  the  Si „  and  are  observed  by 

the  decrease  in  [$i_  ].  The  generation  must  involve  internal 
bd 

vacancy  sources  since  for  the  anneal  temperature  and  the  sample 
thicknesses  used,  the  vacancy  diffusion  rate  is  far  too  small 
for  the  interior  of  the  samples  to  be  In  equilibrium  with  the  surface. 
The  bulk  generation  of  acceptor  vacancies  during  annealing  Is  favored 
by  the  fact  that  the  samples  are  strongly  extrinsic  and  n-type 
at  the  anneal  temperature  employed.  Although  this  explanation 
is  attractive  it  must  be  emphasized  that  there  is  no  direct 
experimental  LVM  evidence  that  Sigfl-Vga  formation  is  responsible 
for  the  decrease  In  [S1Ga]  and  therefore, the  mechanism  Is  regarded 
as  speculative. 
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C.  Microstructural  Studies 

The  effects  of  annealing  on  the  defect  structure  of  heavily  Si-doped 
GaAs  have  been  examined  by  Narayanan  and  Kachare^  who  used  thin  foil 
transmission  electron  microscopy  (TEM)  for  samples  taken  from  different 
regions  of  the  same  GaAs  ingot  as  was  used  by  KS.^  In  brief,  this  study 
revealed  that  the  annealing  treatments  are  accompanied  by  pronounced 
changes  in  the  microstructure  of  the  material.  Annealing  at  1100°C 
followed  by  quenching  to  room  temperature  produced  a  homogeneous  solid 
solution  which  was  free  of  observed  microstructure;  however,  subsequent 
annealing  of  these  samples  at  a  temperature  in  the  range  of  600  to  850°C 
induced  prismatic  dislocation  loops  of  the  vacancy  type.  It  was  found 
that  the  density,  size,  and  distribution  of  the  loops  exhibited  an 
anomalous  temperature  dependence.  Samples  annealed  at  750°C  had  a 
high  density  of  small  loops  while  those  annealed  at  600°C  exhibited  a 
smaller  density  of  larger  loops.  Samples  annealed  at  400°C  contained 
no  dislocation  loops.  In  order  to  explain  this  behavior  it  was 
postulated  that  an  appreciable  fraction  of  the  quenched-in  vacancies 
become  trapped  by  the  Sigfl  defects  to  form  Siga  -  V^a  at  the  lower 
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annealing  temperature  so  that  the  vacancies  are  not  available  for  the 
nucleation  of  prismatic  loops. 

Qualitatively  similar  microstructural  observations  have  been 
reported  for  heavily  je-/44"49^  Se-/50^  and  S-doped^49^  GaAs.  In  a 
detailed  study^4^  of  heavily  Te-doped  material,  dislocation  loops  were 
observed  after  annealing  at  a  temperature  as  low  as  400°C.  The  loop  area 
per  unit  volume  was  found  to  be  constant  for  all  annealing  temperatures 
between  400  and  1000°C.  The  conclusions  of  the  various  investigations 
regarding  the  nature  of  the  prismatic  loops  in  heavily  Te-doped  samples 
are  not  mutually  consistent.  Interstitial-type^46’49^  and  vacancy-type^44’45^ 
loops  have  been  reported  by  different  investigators.  Thermal  annealing 
does  not  produce  dislocation  loops  in  heavily  Zn-doped,^49'5^in  undoped, 
or  in  lightly  Si-^  or  Te-doped  GaAs  samples. 

D.  Yield  Stress 

Swaminathan  and  Copley  (SC)^54^  investigated  the  mechanical  behavior 
of  GaAs  single  crystals  doped  with  various  impurities.  The  temperature 
range  used  was  420  to  680°K.  As  observed  in  earlier  studies ^2,53) 
doping  with  a  donor  impurity  increased  the  yield  stress,  while  an  acceptor 
decreased  the  yield  stress  relative  to  the  undoped  state.  To  explain 
the  donor  behavior,  it  had  been  previously  proposed  that  dislocations 
were  pinned  by  the  impurity  atmosphere.  Because  of  the  absence  of 
either  abrupt  or  serrated  yield  behavior  in  their  stress-strain 
experiments,  SC^54^  claimed  this  explanation  was  unacceptable,  and  they 
proposed  that  the  observed  hardening  was  due  to  the  interaction  of  moving 
dislocations  with  static  defect  complexes.  To  explain  the  magnitude  of 
hardening,  they  proposed  that  the  defect  complexes  were  solute-vacancy 
pairs  similar  to  those  known  to  strengthen  alkali  halide  crystals. ^55^ 
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If  the  solute-vacancy  pair  consists  of  a  Si  atom  and  a  vacancy  located 
on  a  second  nearest  neighbor  site,  i.e.  Siga  -  VGa,  then  the  elastic 
distortion  resulting  from  the  pair  has  deviatoric  as  well  as  hydrostatic 
components,  and  there  will  be  an  elastic  interaction  between  the  pair 
and  both  edge  and  screw  components  of  the  glide  dislocations.  A  SiGfl  -  V^s 
pair  would  not  have  the  type  of  elastic  distortion  necessary  for  the 
interaction.  Also,  in  GaAs,  a  Si  atom  occupying  a  Ga  site  is  a  single 
donor, ^  and  there  is  evidence  suggesting  that  a  gallium  vacancy  may  be  a 
single  acceptor. Thus,  there  may  be  a  short  range  electrostatic 
interaction  between  pairs  and  dislocations  associated  with  the  displacement 
of  the  oppositely  charged  defects  relative  to  each  other,  while  they  are 
at  the  dislocation  core,  produced  by  the  shearing  action  of  the  gliding 
dislocations.  An  interaction  of  this  type  has  been  proposed  to  explain 
impurity  hardening  in  ionic  crystals.^57^ 

The  effect  of  various  thermal  treatments  on  the  yield  stress  and 
dislocation  substructure  of  heavily  Si-doped  GaAs  have  also  been 
investigated  by  SC.^  The  material  they  investigated  was  taken  from 
a  source  ingot  which  had  a  somewhat  higher  [Si]  than  the  material 
investigated  by  KS^’^  and  by  Narayanan  and  Kachare.^  SC  measured 
samples  which  were  first  annealed  at  1100°C  and  then  quenched  to  room 
temperature.  The  samples  were  subsequently  aged  at  different  temperatures. 
They  found  that  the  yield  stress  measured  at  400°C  increased  with  aging 
time  and  approached  a  saturation  value.  The  saturation  value  decreased 
with  increasing  aging  temperature.  This  behavior  is  consistent  with 
the  proposed  interaction  between  moving  dislocations  and  static  Siga  -  VQa 
pairs.  The  yield  stress  would  be  expected  to  increase  during  isothermal 
aging  because  of  the  formation  of  the  pairs  and  to  approach  a  saturation 
value  corresponding  to  the  establishment  of  the  steady-state  concentration 
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of  pairs  at  the  aging  temperature.  This  concentration  of  pairs 
is  expected  to  decrease  with  increasing  aqing  temperature 
because  the  thermal  energy  increases  relative  to  the  energy  bonding 
the  Si  atoms  to  the  vacancies.  Thus  the  saturation  value  of  the  yield 
stress  also  would  be  expected  to  decrease.  This  interpretation  was  qualit¬ 
atively  consistent  with  the  LVM  and  microstructural  behavior  just 
discussed. 

III.  Experimental  Methods 

Figure  1  is  a  flow  diagram  which  indicates  the  measurements  made  in 
the  present  study  and  their  sequence.  The  [Si]  range  indicated  at  the 
top  of  this  figure  includes  the  values  for  three  different  horizontal 
Bridgman-grown  ingots  which  were  doped  with  different  [Si]  as  indicated 
by  the  second  column  of  Table  II.  The  [Si]  in  Table  II  are  those  in 
the  first  portion  of  the  ingot  solidified.  Since  the  equilibrium 
segregation  coefficient^5®^  for  Si  in  GaAs  is  0.14,  samples  with  different 
[Si]  can  be  obtained  by  cutting  them  from  different  parts  of  the  ingot. 

In  the  present  experiments  a  large  number  of  samples  of  each  of  three 
different  [Si]  were  desired.  Therefore  ingots  of  several  hundred  grams 
were  grown  in  long  boats  and  approximately  the  first  one-third  of  each 
ingot  was  used.  In  this  portion  the  [Si]  should  change  slowly  with  the 
fraction  solidified  if  nearly  ideal  growth  conditions  are  approximated. 
Thus  samples  taken  from  the  same  region  in  the  front  portion  of  the  same 
ingot  should  have  nearly  identical  values  of  [Si]. 

The  total  [Si]  was  estimated  in  three  ways.  The  first  estimate 
was  based  on  the  dopant  added  to  the  melt  and  the  known  segregation 
coefficient.  The  second  method  was  by  using  spectrochemical  analyses. 

The  third  and  most  reliable  method  for  the  more  heavily  doped  samples 
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was  a  measurement  made  by  using  an  electron  probe  microanalyzer. 

The  two  ingots  having  the  lower  [Si]  were  single  crystals  and  were 
grown  by  the  authors.  The  ingot  of  highest  [Si]  was  obtained  from 
Crystal  Specialties,  Inc.^59^  and  was  polycrystalline  but  with  very  large 
single  crystal  sections,  i.e.  typical  crystallite  dimensions  were 
centimeters.  In  this  latter  case  groups  of  samples  having  nearly  the 
same  [Si]  were  always  taken  from  the  same  region  of  the  same  crystallite. 

The  total  number  of  samples  required  from  all  three  ingots  for  all 
of  the  measurements  indicated  in  Fig.  1  is  quite  large,  approximately 
140  samples  not  including  those  used  for  mechanical  measurements  by 
Ganesan  and  Copley^.  The  latter  will  not  be  discussed  in 
this  section  as  they  have  already  been  described. 

To  obtain  the  samples  for  the  ng,  LVM,  and  TEM  measurements, 
wafers  approximately  1,5  mm  thick  were  cut  from  the  ingots  with  the 
large  area  surfaces  perpendicular  to  the  ingot  axis.  Each  wafer  was 
then  divided  into  a  number  of  individual  samples  having  typical  large  area 
surface  dimensions  of  3x10  mm  .  For  the  KM-16  ingot  the  measurements  of 
[Sij  indicated  that  the  values  within  a  given  wafer  were  essentially 
identical.  The  most  heavily  doped  ingot,  WA-Si,  had  a  factor  of~2 
change  between  the  [Si]  value  near  the  bottom  of  a  wafer,  i.e.  next  to  the 
boat,  and  the  top.  This  meant  that  for  the  WA-Si  ingot,  samples  of 
the  same  [Si]  had  to  be  taken  from  the  same  regions  of  adjacent  wafers. 
Fortunately  this  ingot  was  the  largest  one  and  frequently  more  than  15 
samples  could  be  cut  from  a  single  wafer. 

As  indicated  in  Fig.  1  a  few  samples  were  measured  in  the  as-grown 
state.  Except  for  these  cases  all  samples  were  heat  treated  in  small 
evacuated  quartz  ampoules  at  1200°C  for  lh  and  then  quenched  to  room 
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temperature  {a  process  hereafter  called  HT+Q)  prior  to  making  any 

measurement  or  proceeding  with  the  isothermal  anneals.  Both  the  quartz 

and  samples  were  subjected  to  special  washing  and  cleaning  procedures 

to  reduce  the  possibility  of  Cu  contamination  during  HT+Q.  As  mentioned 

earlier  the  HT+Q  treatment  establishes  a  conmon  thermodynamic  background 

for  all  samples.  Also,  as  will  be  discussed,  this  process  significantly 

improves  the  microstructural  quality  of  the  material,  particularly  for 

the  most  heavily  Si -doped  samples.  The  carrier  densities  after  HT+Q 
18  -3  19  -3 

ranged  from  1.0x10  cm”  to  1.0x10  cm  .  The  HT+Q  procedure  has 

18  -3 

essentially  no  effect  on  ng  for  material  with  [Si]  i  2x10  cm  and 
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ng  <1.5x10  cm  ,  but  it  substantially  increases  ng  for  as-grown  samples  with 

19  -3 

[Si]  *1x10  cm”  .  In  some  cases  the  as-grown  nfi  almost  doubled  as  a 
result  of  the  HT+Q  treatment. 

After  the  HT+Q  treatment  the  ng  for  each  sample  was  measured  by 
determining  the  frequency  of  the  free  carrier  plasma  minimum  observed  in 
the  infrared  reflectivity.  This  method  has  been  compared  in  detail 
previously^3, ^  with  Hall  measurement  results  and  very  good  agreement 
was  obtained.  In  making  these  measurements  it  was  necessary  to  remove 
surface  material  from  each  sample  after  the  HT+0  treatment  to  insure 
that  the  measured  ng  was  the  bulk  value.  It  was  found  that  removal 
of  <v0.1  wn  by  lapping  and  chemical  polishing  gave  ng  values  which  did 
not  change  with  further  removal  of  material.  A  similar  removal  of  at 
least  0.1  mm  sample  surface  layer  is  also  performed  after  annealing. 

Since  a  large  number  of  samples  of  essentially  identical  [Si] 
and  ng  were  needed  from  each  ingot,  the  samples  from  the  same  wafer 
or  neighboring  wafers  were  formed  into  groups  Identified  by  their  common 
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origin,  the  post  HT+Q  value  of  ng,  and  the  [Si].  Again  following 
Fig.  1,  a  few  samples  from  each  group  were  used  for  TEM  and  liquid 
nitrogen  temperature  infrared  LVM  measurements.  These  and  all  other 
samples  used  for  the  infrared  measurements  were  e"-irradiated  with 
1.4  MeV  electrons  prior  to  the  measurements.  The  sample  temperature 

was  kept  below  140°K  during  the  irradiation  and  the  e"-fluence,£  IxlO1^ 

2 

e  /cm  .  The  LVM  measurements  of  these  HT+0  samples,  which  were  not 
subsequently  isothermally  annealed,  have  been  used^11^  to  obtain  the 
absorption  cross  sections  for  the  four  different  Si  defects  as  mentioned 
previously.  The  TEM  measurements  were  made  for  both  e~-1rradiated 
and  non-irradiated  samples,  and  no  changes  in  microstructure  attributable 
to  the  irradiation  were  observed. 

The  remaining  samples  were  sealed  in  quartz  ampoules  as  described 
previously  and  taken  through  the  isothermal  anneals  indicated  in  Fig.  1. 
The  groups  were  arranged  so  that  several  samples  of  each  group  could 
be  annealed  for  different  times  at  each  temperature.  The  detailed 
schedule  used  is  given  in  Table  II.  Unless  otherwise  indicated  each 
anneal  time  and  temperature  combination  required  two  samples,  one  for 
the  post-annealing  measurement,  of  nfi  and  then  to  be  electron  irradiated 
and  measured  for  the  LVM  absorption,  the  other  for  TEM  analysis.  Again 
additional  samples  were  required  for  mechanical  measurements, and  they 
have  not  been  included  in  this  description.  In  a  few  cases  the  annealed 
samples  used  for  ng  and  infrared  measurements  were  subsequently  TEM 
analyzed  to  determine  if  the  e”-irradiation  produced  changes  in  the 
annealed  sample  microstructure.  As  in  the  non-annealed  samples  the 
comparisons  made  here  and  those  made  previously^^  for  the  annealed  and 
Irradiated  samples  with  those  annealed  but  non-irradiated  showed  that  no 
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TEM  sensitive  microstructural  differences  occurred. 

The  infrared  absorption  measurements  were  made  with  a  single¬ 
beam  grating  spectrometer  equipped  with  a  set  of  transmission  filters 
to  isolate  the  first  order  spectrum.  The  transmission  was  measured 
in  a  low  temperature  optical  dewar  by  using  the  sample  in-sample  out 
method  with  the  sample  mounted  over  one  of  a  pair  of  matched  aperatures. 

The  spectral  range  covered  was  350  cm"1  £  v  £  500  cm"1  which  contains 
all  of  the  absorption  bands  of  interest. 

In  a  few  cases,  larger  frequency  ranges  were  covered  to  ensure  that  no 
new  bands  were  being  introduced.  The  samples  used  for  the  infrared 
measurements  were  generally  about  0.3  mm  thick  and  had  a  wedge  of  -xO. 2° 
to  eliminate  multiple  reflection  interference  fringes.  The  absorption 
coefficient  a  was  calculated  from  the  transmission  T  by  using  the 
standard  expression: 

T  =  ^ample^o  =  0-R)2exp(-ax)/l-R2  exp(-2ax), 

where  R  is  the  reflectivity  and  x  is  the  sample  thickness. 

All  of  the  samples  for  TEM  studies  ~e  lapped  and  polished  on 

both  surfaces  to  a  thickness  of  x0.2  mm.  Specimens  of  %3x3  mm2  were 

cut  and  were  thinned  by  using  a  chemical  jet  polishing 
(62l 

technique'  '  to  obtain  samples  suitable  for  the  TEM  measurements.  The 
specimens  were  examined  by  a  Hitachi  HU-125C  transmission  electron 
microscope  operated  at  125  KV.  TEM  measurements  of  the  HT+Q,  Si -doped 
samples  as  well  as  undoped  samples  showed  them  to  be  free  of  microstructure 
except  for  an  occasional  dislocation  loop.  In  particular  the  samples 
showed  no  stacking  faults  or  precipitates  while  Si -doped  samples  which 


had  not  been  given  the  HT+Q  treatment,  i.e.  in  the  as-grown  condition, 
were  found  to  have  frequent  extrinsic  stacking  faults,  particularly 
those  samples  having  the  largest  [Si],  I.e.  WA-Si. 

IV.  Results 

The  body  of  experimental  data  to  be  presented  here  is  organized  in 
sections  as  follows: 

A.  Annealing-induced  changes  in  ng 

B.  Annealing-induced  changes  in  LVM  absorption 

C.  Annealing-induced  effects  in  the  TEM  measured  microstructure. 
In  each  case  the  data  will  include  the  results  for  the  three  [Si] 

and  for  all  anneal  temperatures.  The  results  of  mechanical  measurements 
will  be  discussed  later. 

A.  Annealing-induced  changes  in  n  : 

Ingot  KM-26,  [Si]  ^  1.5x1QT8  cm"3,  nc  ^  l.OxlO18  cm3 

Measurements  of  samples  cut  from  this  ingot  and  given  the  HT+Q 
treatment  showed  no  change  in  ng  after  annealing  at  400,  500,and  700°C. 
Ingot  KM-16,  [Si]  1  2.5x10^  cm~^,  nc  1  5-7x10^  cm~^  and 
Ingot  WA-Si,  [Si]  13-6x10^  cm"3,  nc  17-9xi018  cm3. 

Samples  from  these  ingots  show  dramatic  annealing-induced  changes 
in  ng,  and  the  results  are  shown  in  Figs.  2,  3, and  4.  The  T^  >  600°C 
anneal  curves  appear  to  have  reached  equilibrium  by  approximately  150h 
while  the  400  and  500°C  curves  have  not  reached  equilibrium  for  the 
longest  anneal  times  used.  Decreases  in  ng  approaching  an  order  of 
magnitude  are  observed.  Even  though  the  equilibrium  values  are  not 
reached  at  the  lower  temperatures,  it  appears  from  the  data  that  the 
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equilibrium  value  for  ng  increases  with  increasing  T^.  From 

(41 

previously  reported  measurements'  '  it  is  known  that  the  annealing- 
induced  changes  are  reversible,  i.e.  a  HT+Q  treatment  restores  ng 
to  the  pre-anneal ed  value. 

B.  Annealing-Induced  Changes  in  LVM  Absorption: 

Ingot  KM -26 

Because  of  the  relatively  low  [Si]  in  this  ingot,  only  the  Siga 
and  Si^s  bands  at  384  and  399  cm’1  were  observed.  These  absorption 
bands  showed  no  measureable  change  in  integrated  strength,  0pA,  as  a 
result  of  annealing  at  any  of  the  temperatures.  The  ng  values  calculated 
from  the  free  carrier  plasma  minimum  in  the  reflectivity  before  e’- 
irradiation  are  known  to  be  close  to  those  obtained  from  Hall  measurements 
and  are  referred  to  as  "measured"  values.  By  using 
floptical)  =  a(apA)384  -  b(apA)3gg»  where  a  =  lO.OxlO16  cm"1  and 
b  =  12.7xl016  cm-1,  values  of  na  (optical)  can  be  calculated  from  the 
LVM  data  and  they  are  in  good  agreement  with  the  measured  values.  The 
comparison  will  be  discussed  later.  Hereafter  ne(optical)  refers  to  the 
ng  calculated  from  the  LVM  data.  It  may  also  be  noted  that  [$iga]=a(apA)3g4 
and  [Si^]  =  Map^ggj  and  thus,  the  LVM  results  can  be  presented  in 
terms  of  the  concentrations  of  the  related  defects. 


Ingot  KM-16  and  Ingot  WA-Si 


The  LVM  measurements  were  made  with  the  same  samples  used  for 
the  ng  ■  measurements  described  above.  As  an  illustrative  example,  the 

infrared  absorption  specturm  of  a  500°C/1200h  annealed  sample  from  KM-16 
is  shown  in  Fig.  5  along  with  the  spectrum  for  an  identical  but  non- 
annealed  sample.  The  band  frequencies  observed  in  Fig. 5  are  the  same  as  those 
observed  previously  and  given  in  Table  I.  In  this  illustrative  case  it 
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is  seen  that  the  anneal  causes  changes  in  all  of  the  LVM  bands  and 
hence  in  the  concentrations  of  the  related  defects.  One 
can  calculate  ne(optical)  from 

ne(optical)  =  a(y )3g4  -  b(apA)399  +  c(apA)367  -  d(«pA)369 

*  Piga]  "  +  ^367  cm_1{1efect]  "  [369  cm'^defect]^ 

where  a  and  b  are  the  same  as  above,  c  =  25.3x10^  cm  \  and  d  =  19.0x10^ 
cm"^.  Again, each  term  on  the  right  is  the  concentration  of  the  related 
defect.  The  LVM  strengths  for  all  annealing  conditions  have  been  measured 
and  some  of  the  results  are  given  In  terms  of  defect  concentrations  in  Figs. 6 
through  9.  Since  the  constant  of  proportionality  between  [Siga-$iAs]  ar)d 
(apA)393  is  not  known,  the  [Sigg-Si^]  in  these  figures  is  not  the 
actual  concentration  but  proportional  to  it.  Figures  6  through  8  give 
the  results  for  ingot  KM-16  and  Fig. 9  is  for  WA-Si  after  700°C  annealing. 

The  relationship  between  ng  (optical)  and  ng  (measured)  for  all 
non-annealed  samples  is  shown  in  Fig.  lO.The  straight  line  of  slope 
one  indicates  ng  (optical)  =  ng  (measured),  and  as  it  is  seen  in  this 
figure  all  non-annealed  samples  lie  close  to  this  line.  In  one  sense 
this  agreement  is  not  surprising  since  the  LVM  measurements  for  these 
samples  are  those  which  were  used^1^  to  determine  a,  b,  c,  and  d  which 
yielded  ng  (optical).  However,  the  close  fit  of  all  18  points  is 
significant  when  one  realizes  that  there  were  only  4  adjustable  parameters. 

The  changes  in  ng  (optical)  with  annealing  time  are  illustrated  for  2 
different  values  of  T^  in  Fig.  Hand  compared  with  the  annealing 
curves  for  ng  (measured)  for  the  same  samples  and  conditions.  The  two 
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400°C  anneal  curves  are  very  close  to  one  another  while  those  for 
500°C  strongly  diverge.  The  results  of  similar  measurements  for  all 
ingots  and  annealing  conditions  are  summarized  in  Figs. 12  and  13, 
where  n  (optical)  is  plotted  against  n  (measured).  In  all  cases  where 
the  points  do  not  lie  close  to  the  lines  of  slope  one  they  are  high 
because  the  decrease  in  ng  (optical)  is  too  small  to  account  for  the 
observed  decrease  in  ng  (measured). 

C.  TEM  Measured  Microstructure: 

As  mentioned  previously,  the  HT+Q  treated  samples  which  were  not 
annealed  all  had  no  significant  microstructure  for  all  three  ingots. 

Only  a  few  dislocations  were  observed.  In  particular  there  were  no 
stacking  faults,  no  dislocation  loops  and  no  precipitates.  After 
annealing,  some  significant  changes  were  seen  for  the  more  heavily  doped 
samples. 

Ingot  KM-26 

Similar  to  the  ng  and  LVM  results  for  this  ingot  of  lowest  [Si], 
there  were  no  microstructural  changes  observed  for  any  of  the  annealing 
temperatures  or  times.  The  samples  remained  essentially  featureless. 
Ingot  KM-16  and  Ingot  WA-Si 

Figure  14  shows  a  set  of  micrographs  for  samples  from  Ingot  KM-16 
which  were  given  the  HT+Q  treatment  and  then  annealed  at  different 
temperatures.  No  dislocation  loops  were  observed  for  annealing  temper¬ 
atures  of  400  and  500°C.  As  T^Z  600°C  a  high  density  of  loops  is 
observed.  In  general  the  loop  size  increases  and  the  loop  density 
decreases  for  increasing  T^.  Ingot  WA-Si  has  a  similar  behavior  except 
that  the  total  loop  area  per  unit  volume  is  larger.  Tables  III  and  IV 
summarize  the  TEM  results  for  the  annealed  samples  from  the  two  ingots. 
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The  loop  densities  and  areas  are  rough  estimates  from  micrographs. 

For  600°C  £  £  800°C  anneals,  increasing  the  anneal  time  increases 

the  loop  size  and  decreases  the  density.  For  each  ingot  the  total 
loop  area  shows  some  increase  for  short  anneal  times  and  then  tends  to 
remain  essentially  constant  for  longer  annealing  times  at  600,  yOO^and 
800°C.  The  total  loop  area  appears  to  start  to  decrease  at  900°C  and  by 
1000°C  is  significantly  decreased  for  both  ingots. 

The  effect  of  the  anneal  temperature  and  time  on  the  habit  planes 
of  the  loops  is  also  indicated  in  Tables  III  and  IV.  A  short  term 
anneal  produces  small  loops  lying  on  {111}  planes  with  b  =  a/3[l 11]. 
Increasing  the  anneal  time  increased  the  loop  size  and  the  loops  were 
increasing  on  {110}  planes  with  b  *  a/2[110].  For  annealing  temperatures 
T^  £  800°C  and  times  >  150h  the  {111}  loops  had  disappeared  and  the  {110} 
loops  increased  in  size  and  decreased  in  density. 

All  of  the  above  TEM  results  are  similar  to  those  previously  reported 
for  annealed,  heavily  Te-doped  GaAs.^47^  The  primary  differences  are 
that  the  loops  in  the  Te-doped  material  are  formed  at  T^  as  low  as  380°C 
and  the  loop  area  is  constant  for  all  T^  up  to  1000°C  and  for  all  times. 

An  analysis^)  of  the  TEM  results  has  been  made  for  different  T^ 
by  using  12  of  the  annealed  samples  from  both  ingots  In  order  to  determine 
the  nature  of  the  loops.  In  all  cases  the  loops  were  found  to  be  of  the 
extrinsic  type.  This  result  differs  from  that  previously  given  by  Narayanan 
and  Kachare'  '  for  annealing-induced  loops  in  similar  material.  They  found 
the  loops  were  of  the  intrinsic  type,  i.e.  vacancy  loops.  In  order  to 
check  the  present  results,  a  different  and  independent  TEM  analysis^^ 
was  doubly  conducted  on  some  of  the  above  12  annealed  samples,  which  show 
faulted  loops  on' {111}.  Moreover,  two  additional  annealed  samples 
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were  sent  to  the  Boeing  Aerospace  Corp.  (Seattle,  Wash.)  for  a 
second  analysis  with  a  different  microscope  (EM30Q)  and  operator. 

The  results  all  confirmed  the  conclusion  that  the  loops  are  of  the 
extrinsic  type. 

In  an  attempt  to  resolve  the  differences  between  the  present 
results  and  the  previous  work  samples  from  KM-16  and  samples  from 
the  ingot  used  by  the  previous  authors  were  given  the  same  HT+Q 
treatment  employed  in  the  previous  study  (1100°C/15  min)  and  then 
annealed  at  exactly  the  same  condition,  i.e.  750°C/320h  which  was 
an  annealing  condition  of  the  previous  study.  The  TEM  analysis^3^ 
were  carefully  conducted  and  the  results  confirmed  that  the  loops 
are  indeed  of  the  extrinsic  type.  We  have  therefore  been  unable  to 
reproduce  the  results  of  the  previous  study  even  when  the  experimental 
conditions  are  the  same  as  those  specified  by  the  previous  authors. 

In  all  of  the  present  work  for  all  heavily-SI  doped  ingots  the  thermal 
annealing-induced  loops  are  always  extrinsic- 


In  a  previous  study^^of  Li-saturated  GaAs:Si,  where  the  Li 

saturation  was  done  at  temperatures  in  the  range  of  the  values 

used  here  and  where  the  same  material  used  here  was  also  used  for  that 

study,  extrinsic  loops  were  also  observed  and  scanning  electron 

microscope  measurements  employing  an  energy  resolution  system  established 

that  the  loops  were  Si-rich.  By  assuming  the  loops  in  the  present  study 

contain  double  layers  of  Si  interstitial  atoms,  one  can  roughly  estimate 

the  Si  content  in  the  loops,  i.e.  CSi ] ^ 00p »  ^or  eac*1  annealed  sample. 

These  estimates  are  also  given  in  Tables  III  and  IV.  Further  justification 

of  this  assumption  was  supplied  by  one  KM-16  sample  which  was  double 

annealed  at  950°/2h+700°C/65h.  Large  extrinsic  and  triangular  loops 

o 

having  a  size  of  approximate  7000  A  were  observed.  The  loop  size  is 
similar  but  even  larger  than  that  observed  for  the  900°C  annealed  case. 
Microanalysis  of  these  large  loops  with  a  scanning  electron  microscope 
showed  that  they  have  a  large  Si  content  similar  to  the  loops  in  the  Li 
saturated  GaAs:Si.  We  therefore  assume  that  all  the  observed  extrinsic 
loops  are  due  to  Si  precipitation  from  the  GaAsrSi  solution. 

V.  Discussion 

The  first  part  of  this  section  examines  the  conclusions  which  can 
be  drawn  from  the  results  of  the  previous  section  and  the  relationships 
between  the  different  measurements.  As  will  be  seen  the  rather  extensive 
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body  of  experimental  data  permits  one  to  examine  the  behavior 
of  a  relatively  complicated  system.  The  discussion  is  unusual  in  that 
the  measurements  allow  conclusions  to  be  formed  on  the  basis  of  the 
knowledge  of  the  changes  of  individual  defect  concentrations  rather 
than  considerations  only  in  terms  of  quantities  such  as  ng  which  depend 
upon  the  net  behavior  of  two  or  more  defect  concentrations.  The  second 
part  of  this  section  will  utilize  the  results  in  an  attempt  to  interpret 
the  annealing-induced  changes  in  mechanical  properties. 

A.  Carrier  Density,  LVI^and  Microstructure: 

Ingot  KM-26 

in  this  relatively  lightly  Si -doped  ingot  there  appears  to  be  no 
detectable  effect  of  annealing.  The  values  of  ng,  [Sigg],  [Si^]  and 
the  microstructure  are  all  unchanged  by  the  longest  annealing  times  at 
all  temperatures.  This  result  is  in  agreement  with  the  previous  less 
detailed  results  for  ng  by  KS^  for  GaAstSi  with  a  similar  [Si]. 

Ingot  KM-16  and  MA-Si 

In  the  more  heavily  doped  ingots  there  are  a  large  number  of 
annealing-induced  changes.  In  general,  the  pattern  of  the  chanqes  are 
similar  for  both  ingots.  The  results  will  be  examined  for  each 
annealing  temperature. 

(4Q0°C)  After  annealing  at  this  temperature  there  are  essentially 
no  loops  or  other  microstructural  features  although  both  [Sigfl]  and  ng 
showed  large  decreases.  The  other  Si  defect  concentrations  including 
[Si^],  tSiQa-Si^sJ,  [367  cm-1  defect],  and  [369  cm'1  defect]  all 
remained  constant  or  showed  relatively  small  changes  even  for  the  samples 
with  the  highest  [Si].  It  is  also  noted  in  Figs.  12  and  13  that  all  of 
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the  400°C  points  are  close  to  the  straight  lines  of  slope  one  indicating 
that  ng(optical)  is  essentially  equal  to  ng(measured)  for  all  anneal 
times  and  for  both  ingots.  This  conclusion  can  also  be  reached  by 
considering  Fig. 11.  Therefore  the  decrease  In  ng(measured)  with  anneal 
time  is  attributable  entirely  to  the  decrease  in  [Siga]  donors.  It  Is 
also  observed  that  the  decrease  In  $igg  defect  concentration  cannot  be 
allocated  to  increases  in  any  of  the  other  observed  defect  concentrations 
In  Fig. 6  and  similarly  for  WA-S1  for  the  same  anneal  conditions  unless 
one  assumes  that  each  defect  Involves  several  Si  atoms.  The  latter 
possibly  can  be  rejected  since  earlier  Si  Isotopic  studies^)  showed 
that  the  367  and  369  cm"1  bands  Involve  one  Si  atom  per  defect.  Since 
no  precipitates  or  extrinsic  loops  are  observed  in  the  400°C  annealed 
material  we  attribute  the  change  in  [Sigfi]  to  the  formation  of  a  new 
neutral  point  defect.  It  has  been  proposed  previously^4-6^  that  this  new 
point  defect  results  from  the  generation  of  Vgg  acceptors  followed  by 
a  pairing  reaction  to  forms  s^ga_vGa-  It  is  clear  that  this  new  defect 
cannot  be  responsible  for  either  the  367  or  369  cm'1  bands,  and  it  is 
therefore  disturbing  that  no  new  bands  attributable  to  the  LVM  of  this 
defect  are  observed.  It  is  possible  that  the  presence  of  the  second 
neighbor  vacancy  lowers  the  effective  force  constants  at  the  Si  site 
and  the  LVM  frequencies  shift  down  Into  the  region  of  high  lattice 
absorption  of  the  GaAs  and  hence  are  not  observable.  For  the  purposes 
of  the  present  study  the  formation  of  Si^-V^  defects  Is  regarded  as 
the  mechanism  responsible  for  the  decrease  in  [Siga].  There  is  however  no 
evidence  here  that  the  Siga-Vga  formation  is  indeed  responsible  for  the 
changes  In  [S1ga]. 

(500°C)  As  after  annealing  at  400°C  the  samples  have  no  observable 
microstructural  features  and  in  particular,  no  loops  or  faults.  The  defect 
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behavior  is  much  more  complicated  than  that  at  400°  C  in  that  all  of 
the  defects  change  concentration  in  significant  amounts.  The  [Sigg], 
[SiA$],  and  ng  (measured)  decrease  while  all  of  the  other  defects  show 
increases,  especially  for  the  longest  anneal  times.  Of  particular 
interest  is  the  increase  in  CSiQa-Si^s3  after  1200h  which  indicates  that 
there  is  detectable  Si  diffusion  at  this  temperature.  Of  the  3  temper¬ 
atures  studied  in  detail,  400,  500,  and  700°C  the  only  temperature  at 
which  changes  in  Siga-Si^s  pair  concentrations  are  observed  is  at  500°C. 
KS^  have  previously  reported  this  effect  for  600°C  annealing.  Assum¬ 
ing  the  [Si]  to  be  uniform  and  the  diffusion  length  to  be  equal  to  the 
distance  between  two  nearest  Si  atoms,  a  diffusion  constant  of  about 
3xl0‘20  cm2  sec-1  at  500°C  is  obtained.  This  value  is  in  reasonable 
agreement  with  that  estimated  by  KS^  from  earlier  data. 

As  indicated  in  Figs.  11,  12^and  13,  n£  (optical)  does  not  equal 
ng  (measured)  after  annealing  at  500°C.  In  particular  the  ng  (measured) 
decreases  substantially  while  nfi  (optical)  has  very  little  decrease.  The 
n&  (optical)  decreases  about  10%  for  KM-16  and  about  30%  for  WA-Si  while 
ng  (measured)  decreases  by  a  factor  of  six  to  eight  for  the  longest 
anneal  times  for  both  materials.  Therefore  one  must  postulate  the 
formation  of  a  new  "unknown  acceptor"  which  is  not  observed  in  the  LVM 
spectra  and  is  not  the  SiQ0-VGa  center  proposed  to  explain  the  400°C 
change  in  [SiGa3 .  The  concentration  of  the  new  acceptor  is  given  by 
n  (optical)  -  n,  (measured)  for  the  longest  anneal  times  and  this 
quantity  is  plotted  in  Fig.  15  as  a  function  of  the  anneal  temperature. 
The  curves  on  the  figure  have  no  significance  other  than  to  smoothly 
connect  the  experimental  points.  If  the  loss  in  [Siga]  plus  the  loss 
in  [Sifls]  is  compared  with  the  gains  in  the  [367  cnf^  defect]  and 


[369  cm”^  defect]  and  if  one  assumes  one  Si/defect  for  the  latter 
two  cases  as  indicated  by  the  prior  isotropic  studies, then 
there  could  be  a  sufficient  net  [Si]  loss  to  account  for  both  the 
increase  in  CSiga~Si^s]  and  the  [unknown  acceptor],  again  assuming 
the  latter  has  one  Si /defect.  Therefore  after  annealing  at  500°C 
one  does  not  have  to  postulate  the  presence  of  any  Si  in  the  system 
other  than  that  present  or  proposed  for  the  point  defects  discussed. 

A  major  uncertainty  here  is  the  lack  of  information  on  changes  in 
[SiGg]  due  to  Siga-Vga  formation.  Moreover,  as  will  be  seen  a  different 
conclusion  seems  to  be  required  from  the  700°C  annealing  data. 

(700°C)  At  this  annealing  temperature  microstructure  is  intro¬ 
duced  in  the  form  of  loops.  The  loops  density  decrease  and  the  average 
loop  size  increases  with  time,  see  Tables  III  and  IV.  The  loops  are 

of  the  extrinsic  type.  There  is  no  listing^65^  of  any  Ga  Si  or 

x  y 

Si  As  compounds,  and  therefore  equal  numbers  of  Ga  and  As  inter- 
x  y 

stitials  or  of  Si  interstitials  are  the  only  two  likely  possibilities 
for  the  extrinsic  loops.  As  previously  reported  in  the  study  of  the 
annealing  effect  for  GaAs:Te  material,  the  loops  consisting  of  Ga  and 
As  interstitials  started  to  form  even  for  an  annealing  temperature  as 
low  as  380°C  while  in  the  present  study  the  loops  formed  only  for 
T^  >  600°C.  Because  of  the  lower  diffusion  coefficient  for  Si  in  GaAs 
for  the  temperature  range  of  concern  ,it  is  not  unreasonable  to  find  a 

(471 

larger  T^  is  required  for  extrinsic  loops  formed  of  Si  interstitials.  ' 
Silicon  has  a  diamond  structure  with  a  lattice  parameter  of  5.43  H  which 
is  close  to  that  of  GaAs  and  hence  layers  of  Si  could  be  easily 
accomodated  in  the  GaAs  lattice. 


As  observed  in  Pigs.  12,  13,and  15,the  ng  (optical)  and 

nfi  (measured)  behave  in  a  manner  similar  to  that  observed  at  500°C. 

Again  a  large  [unknown  acceptor]  is  necessary  to  account  for  the  fact 

that  ng  (optical )»  ng  (measured)  after  prolonged  annealing.  However, 

the  individual  defect  concentrations  showed  a  significantly  different 

behavior  at  700°C  than  at  500°C.  The  [Sig  -S1^  ]  remained  constant, 

the  [Siga3  and  [Si^]  both  decreased  in  the  first  0.5  h  and  were 

essentially  constant  thereafter  and  the  367  and  369  cm~^  defects 

increased  somewhat .particularly  in  WA-Si.  Again,  if  it  is  assumed  that 

there  is  one  Si/defect,  then  the  total  loss  of  Si  from  Siga  and  Si^s 

defects  and  the  total  gain  in  the  367  and  369  cm"^  defects  are  both 

approximately  10x10  cm"  .  The  2  to  3x10  cm  Si  atoms  in  the  loops 

is  relatively  small  compared  to  these  changes  (see  Table  IV  for 

WA-Si)  and  the  loops  could  come  from  the  net  changes  in  these  defect 

concentrations.  This  is  an  attractive  explanation  since  the  large 

decreases  in  the  SiGa  and  Sifts  concentrations  and  the  loop  formation 

both  take  place  in  the  first  0.5h  of  annealing.  However,  it  is  very 

unlikely  that  the  changes  in  these  same  point  defect  concentrations 

18  -3 

could  also  supply  the  5x10  cm  Si  atoms  in  the  "unknown  acceptor". 
Therefore  if  the  "unknown  acceptor"  does  involve  Si  it  would  appear  to 
require  a  Si  source  beyond  those  considered  here.  Other  studies^’^ 
have  suggested  that  in  heavily  Si -doped  GaAs,  not  all  the  Si  is  in  the 
defects  responsible  for  the  LVM  spectra. 

(800,900,1000qC)  In  this  temperature  range  all  of  the  annealing- 
induced  changes  reduce  as  one  increases  T^.  The  ng  (measured)  comes 
closer  to  the  pre-annealed  HT+Q  value,  ng  (optical)  becomes  nearly 
equal  to  ng  (measured)  so  that  the  [unknown  acceptor]  becomes  small. 
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the  individual  defect  concentration  changes  are  reduced,  and  the  loop 
area  and  hence  Si  in  the  loops  is  also  reduced. 

B.  Relationship  of  Yield  Stress  to  ng,  LVM.and  TEH  Results. 

Recently,  Ganesan  and  Copley ^ 7 )  reported  yield  stress  measurements 
at  400°C  on  two  sets  of  samples  with  different  [Si]  which  were  taken 
from  two  of  the  ingots  (KM-16  and  WA-Si)  used  for  the  present  study.  The 
samples  for  the  stress  measurements  were  prepared  from  regions  of  each 
ingot  adjacent  to  those  from  which  samples  were  taken  for  the  measurements 
reported  here.  Figures  16  and  17  show  the  yield  stress  as  a  function  of 
annealing  time  for  these  samples  for  T^  =  4C0,  500,  and  700°C.  The 
results  given  in  these  two  figures  are  similar  to  those  of  previous 
measurements.^  An  unusual  feature  seen  in  Fig.  17  is  the  Tft  =  700°C 
curve  for  samples  with  the  largest  [Si],  They  seem  to  have  reached 
saturation  even  before  the  annealing.  Although  the  crystal  orientation 
of  the  samples  from  two  ingots  were  shown  to  be  different,  the  Schmidt 
factors^66)  turn  out  to  be  identical.  Hence  it  can  be  inferred  that 
the  difference  in  yield  stress  values  arises  from  the  [Si]  difference 
between  the  ingots  rather  than  an  orientation  effect.  In  general,  the 
yield  stress  values  before  annealing  increase  with  increasing  [Si]. 

One  primary  conclusion  of  the  previous  section  is  that  the 
decrease  of  [Siga]  after  T^  =  400°C  annealing  is  produced  by  the  generation 
of  Vga  and  the  formation  of  neutral  S1ga“Vga  pair  defects.  This  process  is  also 
responsible  for  the  decrease  in  ng.  As  previously  discussed  it  has 
been  concluded^)  that  these  Siga-Vga  pair  defects  are  responsible  for 
the  dramatic  increase  In  yield  stress  after  annealing.  If  one  were  to 
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assume  for  all  the  decrease  in  [Siga]  =  &[Siga]  arises  completely 
from  the  formation  of  Siga"Vga  pairs,  i.e.  6[Sigg]  =”6^SlGa"VGa^ 
one  can  attempt  to  correlate  the  LVM  results  for  6[Siga]  to  yield 
stress  changes  shown  in  Figs.  16  and  17.  A  customary  relation  is 
of  the  forn/37^  o-o0  =  A( Nq+5 tSiga] )m ,  where  o  is  the  yield  stress, 
aQ  is  the  yield  stress  for  undoped  GaAs,  NQ  is  the  pre-anneal, HT+Q 
[Si6a-VGa],  s[Siga]  is  the  assumed  change  in  [Sl'ga“VGa^and  A  and 
m  are  constants.  Fig.  18  shows  that  a  good  correlation  for  both  ingots 
is  obtained  for  =  400°C  anneals  and  requires  that  m  =  1/4.  The 
value  of  A  and  the  two  NQ  values,  one  for  each  ingot,  are  determined 
by  the  data  of  three  samples:  one  non-annealed  sample  from  each  ingot 
and  a  400°C/1200h  annealed  sample  from  KM-16.  The  calculated  Nq  value 
is  lxlO15  cm"3  for  KM-16  and  3.8xl0,?  cm'3  for  WA-Si .  The  plot  of  the 
resulting  equation  is  shown  by  the  straight  line  in  Fig.  18.  For 
crystals  of  NaCl ,  NaBr,  KCl,and  KBr,  Chin  et  al/55^  experimentally 
showed  that  the  yield  stress  measurements  give  m  =  1/2.  The  defect 
density  was  assumed  to  be  equal  to  the  doping  concentrations  of  divalent 
impurity  ions  such  as  Ca2+,  Sr2+,  Ba2+  added  during  crystal  growth. 

If  however,  the  defect  responsible  for  altering  the  yield  stress 

is  a  vacancy-impurity  pair  similar  to  that  proposed  here  then  a 

2+ 

likely  process  for  pair  formation  involving  the  Ca  is 

A  I  U 

CaNa  +  VNa  ^CaNa"VNa^  ’  where  •>  f »  and  x  mean  singly  positive, 
singly  negative,  and  uncharged  species  with  respect  to  the  charge 

normally  present  at  the  lattice  site.  If  the  concentration  is 

2+ 

determined  by  the  Ca  concentration  then  mass  action  gives 

^CaNa_VNa^  *  kp^CaNa^  ^VNa^  *  VCaNa^'  If  most  of  the  Ca  imPurity 
Is  present  as  CajJa  then  [CajJa]  »  [Ca2+],  and  o-oQ  •«  [Ca2+]’U(CaNa-VNa)x],< 

Thus  the  result  becomes  Identical  to  that  determined  here  where 

the  Si ga”Vga  pair  concentration  is  measured  in  a  much  more  direct  manner. 
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Indeed,  the  results  reported  here  show  the  first  experimental  data 
for  the  dependence  of  the  yield  stress  on  the  measured  concentration 
of  the  specific  defect  responsible  for  strengthening. 

As  indicated  In  Fig. 18,  the  data  of  the  long  term  annealed  samples 
at  500  and  700°C  for  both  ingots  are  not  in  agreement  with  the  results 
of  the  Ta  =  400°C  anneals.  This  indicates  that  the  assumption  that 
decrease  in  [SiGa]  for*  these  samples  after  500°C  or  700°C  annealing 
comes  completely  from  the  formation  of  the  Si Qa~VGa  defect  is  not  valid. 
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From  the  data  in  Fig.  18  it  is  clear  that  the  yield  stress  indicates  that 


only  a  small  fraction  of  the  decrease  of  [Siga]  can  be  resulting  from 


Siga~Vga  ^ormation  for  these  cases. 


This  conclusion  is  consistent 


with  the  LVM  results  for  both  ingots  which  show  that  Siga~Si^s  pair 
defect  and  other  Si  related  defects  grow  after  =  500°C  anneal,  and  the 


Si-rich  loops  as  well  as  the  other  Si  related  defects  grow  after  the 


=  700°C  anneal.  Indeed,  we  have  already  concluded  for  the  same 
values  that  almost  all  the  decrease  in  [Sigfl]  as  well  as  that  for  the 
[SiAsl  was  required  to  explain  the  observed  increases  in  the  other  Si 
point  defects  and  microstructure.  The  most  striking  example  of  the 
correlation  of  the  different  results  is  for  the  T^  =  700°C  annealing  of 
WA-Si.  The  large  decrease  of  [Si6a]  in  the  first  0.5h  was  required  to 
explain  the  extrinsic  loop  formation.  If  the  decrease  in  [Sigfl]  were 
all  by  the  formation  of  [SiGa-VGa]  then  the  final  a-oQ  should  have 
been  a  factor  of  ^2  higher  than  the  observed  value.  The  constancy 


of  [Siga]  after  this  initial  decrease  is  also  consistent  with  the 
constancy  of  the  yield  stress  for  different  annealing  times. 


VI.  Summary 

Measurements  of  the  isothermal  annealing-induced  changes  in  the 
carrier  density,  LVM  spectra,  TEM  microstructure  and  critical  resolved 
shear  stress  of  GaAs:Si  have  been  compared  for  three  different  [Si] 
ingots  at  several  annealing  temperatures.  The  principal  results  of 
this  study  are  that  for  [Si ]  eg  1  to  2x10  cm  there  are  no  significant 
annealing-induced  changes  in  any  of  the  measured  properties  while 
for  [Si]  a  2  to  6x10  cm  there  are  major  changes  in  all  of  the 
properties.  The  detailed  nature  of  the  changes  depends  upon  the  anneal 
temperature  and  time.  The  following  conclusions  all  pertain  to  the 
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heavily  doped  material: 

(1)  For  400°C  annealing,  the  decrease  observed  in  ng  can  be 
accounted  for  quantitatively  by  the  decrease  in  the  Siga 
donor  concentration  as  measured  by  the  LVM  spectra.  None 
of  the  other  Si  defects  observed  in  the  LVM  spectra  change 
concentration  significantly.  The  mechanism  suggested  to 
account  for  the  decrease  in  Si  donors  is  pairing  with 
vacancies  to  form  the  siga-VGa  pair  defects. 

(2)  For  500°C  annealing,  the  decrease  in  ng  cannot  be  attributed 

to  changes  in  any  or  all  of  the  LVM  defects,  and  it  is  necessary 
to  proposed  to  an  "unknown  acceptor"  of  concentration  compar¬ 
able  to  that  of  the  other  defects.  Si  defect  concentrations 
change  with  prolonged  annealing  times,  however  these  changes 
can  be  largely  attributed  to  redistribution  of  the  Si  among 
the  various  defect  species. 

(3)  For  700°C  annealing, extrinsic  Si-rich  loops  are  observed  after 
the  shortest  annealing  timesfand  the  Si  in  these  loops  appears 
to  come  from  decreases  in  the  Sigg  and  Si^s  concentrations 
which  are  also  observed  after  the  shortest  annealing  times. 

Again  a  large  concentration  of  "unknown  acceptors"  is 
necessary  to  account  for  the  decrease  in  the  measured  nfi. 

If  these  acceptors  involve  Si  then  this  Si  must  come  from  a 
source  other  than  those  defects  observed  in  the  LVM  spectra. 


i.e.  the  SiGa,  SiA$, 
Si  defects. 


Si ca~Si as  pairs  or  two  other  unidentified 


(4)  For  annealing  temperatures  above  700°C  the  changes  in  all  of 


the  properties  reduce  as  increases.  The  concentration  of 
unknown  acceptors  is  essentially  zero  by  Tft  %  1000°C. 
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(5)  The  previous  measurements  of  the  critical  resolved  shear 
stress  are  consistent  with  the  results  of  the  annealing- 
induced  changes  in  the  other  physical  properties  if  one 


assumes  a  shear  stress  dependence  of  the  form  a-aQ  = 


A(N0+$[Siga])^4.  The  data  indicate  that  the  decrease 
in  [SiGa]  at  Ta  =  400°C  is  likely  due  to  VGa  generation  followed 
by  S^Ga_VQaPai'r  formation  while  at  higher  T^values  the  formation  of  thi 
defect  pair  can  account  for  only  a  small  fraction  of  the 


decrease  in  the  Siga  concentration. 
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Table  I. 


Observed  LVM  Band  Frequencies  at  Liquid  Nitrogen 
Temperature  in  Si -Doped  GaAs  Compensated  by  e- Irradiation. 


MODE  FREQUENCY  (cm-1) 

ELECTRICAL  CHARACTER 

DEFECT 

384 

Donor 

S1sa 

399 

Acceptor 

SiAs 

393 

464 

} 

SiGa'SiAs 

367 

Donor3 

Si  Related 

369 

Acceptor3 

Si  Related 

aSee  text  for  explanation  of  this  assignment. 


Table  II.  Isothermal  Annealing  Schedule  for  all  Samples. 


Note:  a.  n  values  before  Isothermal  annealing  but  after  HT+Q. 
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AH  loops  observed  are  extrinsic  (interstitial). 


FIGURE  CAPTIONS 


Figure  1 . 
Figure  2. 


Figure  3. 

Figure  4. 

Figure  5. 


Figure  6. 


Figure  7. 


Flow  diagram  showing  the  sequence  of  the  measurements. 

Isothermal  annealing  curves  for  GaAsrSi  samples  from 

ingot  KM- 16  where  [Si]  ^2.5  x  1019  cm”3.  Samples 

having  the  same  number  and  the  same  symbol  have  the 

same  initial  given  by  the  t=0  value. 

Isothermal  annealing  curves  for  GaAsrSi  samples  from 

ingot  WA-Si  where  [Si]  *  3  6x1 0^9  cm”3. 

Isothermal  annealing  curves  for  GaAs:Si  samples  from 

the  ingot  WA-Si  where  [Si]  a.3-6xl0^9  cm”3. 

Infrared  abosrption  spectrum  at  liquid  nitrogen 

temperature  for  two  identical  GaAsrSi  samples  having  [Si] 

19  -3 

a,  2.5x10  cm  .  designates  non-annealed^HT+Q  sample, 

-+-  designates  500°/1200h  annealed  sample.  Both  samples 

were  compensated  by  e~ -irradiation  at  77°K. 

400°C  isothermal  annealing  curves  giving  defect  or 

measured  carrier  concentrations  for  GaAsrSi  samples  with 
19  -3 

[Si]  2s  2.5x10  cm  .  For  convenience,  the  value  of 
a  =  10.0x10^  cm"^  is  used  as  the  constant  of  proportionality 
in  calculating  [SiGa-SiA$].  See  text  for  the  other  [defect] 
calculations. 

500°C  isothermal  annealing  curves  giving  defect  or 
measured  carrier  concentrations  for  GaAsrSi  samples  with 
[Si]i2.5xl0  cm  .  For  convenience,  the  value  of  a  * 
10.0xl016  cm’^  is  used  as  the  constant  of  proportionality 
in  calculating  CSiga-Si^$].  See  text  for  the  othe  [defect] 
calculations. 
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Figure  8. 


Figure  9. 


Figure  10. 


Figure  11 . 


Figure  12. 


Figure  13 


700°C  isothermal  annealing  curves  giving  defect  or 

measured  carrier  concentrations  for  GaAs:Si  samples  with 
19-3 

[Si]  j i  2.5x10  cm"  .  For  convenience,  the  value  of 
a  =  lO.OxlO18  cm"1  is  used  as  the  constant  of  proportionality  in 
calculating  [Siga-Si^$].See  text  for  the  other  [defect]  calculations. 

700°C  isothermal  annealing  curves  giving  defect  or  carrier 

concentrations  for  GaAs.-Si  samples  with  [Si]  2;3-6xl019crTf  3. 

For  convenience,  the  value  of  a  =  lO.OxlO16  cm”1  is  used 

as  the  constant  of  pro  portionality  in  calculating  [Siga-Si^s]. 

See  text  for  the  other  [defect]  calculations. 

ng  (optical)  calculated  from  LVM  data  vs  ng  (measured)  for 

all  HT+Q , non-anneal ed  GaAs:Si  samples.  The  straight  line 

indicates  ng  (optical)  =  ng  (measured). 

The  annealing-induced  changes  in  ng  (optical)  for  [Si] 

19 

2*  2.5x10  .  For  comparison  the  annealing  curves  for  ng 

(measured)  for  the  same  samples  and  conditions  are  also 
shown . 

ne  (optical)  vs  ng  (measured)  for  all  annealed  samples 
where  [Si]  A  1.5xl018  cm"3  and  2.5xl019  cm"3  for  samples 
from  KM-26  and  KM-16,  respectively.  The  straight  line 
indicates  ng  (optical)  =  ng  (measured). 
ng  (optical)  vs  nfi  (measured)  for  all  annealed  samples 
with  [Si]  st  3-6x1 01 9  cm'3.  The  straight  line  indicates 
ng  (optical)  *  nfi  (measured). 
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Figure  14. 


Bright  Field  micrographs  for  samples  for  Ingot  KM-16 
which  were  annealed  at  different  temperatures  after  HT+Q. 
g  =  [220]  for  all  operative  reflections. 

(a)  400°C/2000h  anneal . 

(b)  500°C/1200h  anneal.  A  low  angle  grain  boundary 
contrast  is  shown. 

(c)  600°C/306h  anneal,  dislocation  loops  are  observed. 

(d)  700°C/400h  anneal,  loops. 

(e)  800°C/400h  anneal,  loops. 

(f)  900°C/150h  anneal,  loops. 

(g)  1000°C/150h  anneal,  loops. 

(h)  1200° C/1 h  (i.e.  HT+Q),  no  loops. 

Figure  15.  The  concentration  of  the  "unknown  acceptor"  for  the 

longest  anneal  times  as  a  function  of  anneal  temperature. 

Figure  16.  Yield  stress  vs.  annealing  time  for  samples  from  ingot 
KM-16  for  various  annealing  temperatures.  •  designated 
yield  stress  aQ  of  undoped  GaAs. 

Figure  17.  Yield  stress  vs.  annealing  time  for  samples  from  ingot 
WA-Si  for  various  annealing  temperature.  •  designates 
yield  stress  aQ  of  undoped  GaAs. 


Figure  18.  o-o0  vs.  (NQ+S[SiGa])1/4  for  the  annealed  samples 

with  T^  =  400,  500  and  700°C.  The  straight  line  was 
obtained  as  described  in  the  text.  The  calculated  N 
value  is  IxlO1 5  cm"3  for  KM-16  and  3.8xl017  cm"3  for 
WA-Si . 


I I-B3-44 


KM  — 16 


II-B3-46 


1 1 -B  3-49 


KM -16  400 °C  ANNEAL 


i — i 


I I-B3-50 


I 


0.5  1.5  15  4001200  2000 

ANNEALING  TIME  (HOURS)  — - 


(  luo)  9u  jo  [psjsg] 


I I -B3-51 


ANNEALING  TIME  (HOURS) 


ANNEALING  TIME  (HOURS 


ANNEALING  TIME  (HOURS) 


(optical) 


KM -26  KM  -16 


II-B3-56 


(optical ) 


II-B3-57 


AD-A107  598 


UNCLASSIFIED 

2-2 


UNIVERSITY  OF  SOUTHERN  CALIFORNIA  LOS  AN6ELES  ELECTR— ETC  F/6  20/12 
THE  EFFECT  OF  SOLUTE-VACANCY  PAIRS  ON  THE  PHYSICAL  PROPERTIES  0— ETCCU) 
MAY  81  Hfi  SPITZER »  J  M  WHELAN •  S  M  COPLEY  AFOSR-76-2990 

AFOSR-TR-81— 0746  NL 


END 


I  82 

Dug 


o 


Or-OV 


O 


OJ 

O 


(MPa) 

cn 

o  o 

~I - T' 


CD  -nI 

O  O 

"i - r 


+  •  X 


II-B3-65 


4. 


Strengthening  of  GaAs  Single  Crystals 
by  Formation  of  Silicon-Vacancy  Pairs 


* 

V.  Ganeson  and  S.  M.  Copley 


Department  of  Materials  Science 
University  of  Southern  California 
Los  Angeles,  CA  90007 


Several  investigations  have  shown  that  annealing  can 

produce  dramatic  changes  in  the  properties  and  microstructure 

of  heavily  Si-doped  GaAs.  These  investigations  include: 

1  2  3 

measurements  of  free  carrier  density,  photoluminescence  ' 

3  4 

and  localized  vibrational  maode  (LVM)  infrared  absorption;  ' 

characterization  of  dislocation  substructure  by  transmission 

electron  microscopy ; ^ ^  and  measurement  of  critical  resolved 
6  7 

shear  stress.  '  The  results  of  these  investigations  strong¬ 
ly  suggest  that  the  annealing-induced  changes  are  related 
to  specific  changes  in  point  defect  structure.  Indeed,  a 

qualitative  correlation  has  been  found  between  changes  in 

•  I  x 

these  properties  and  changes  in  the  concentration  of  (SiGa-VGa) 
pairs.  (The  various  defect  species  are  indicated  as  follows: 

Fm  denotes  F  atom  sitting  on  M  site  and  VM  denotes  a  vacancy 
on  an  M  site.  Superscripts  denote  the  charge  state  of  the 
species.  A  dot  represents  a  positive  charge,  a  dash  repre- 


* 
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sents  a  negative  charge  and  an  x  represents  a  neutral 
species.  Square  brackets  indicate  site  fraction.  ) 

The  difficulty  in  attaining  a  quantitative  correla¬ 
tion  of  results  from  previous  investigations  arises, 
because  the  results  were  obtained  employing  different 
ingots  or  different  sections  of  the  same  ingot.  Thus 
the  specimens  had  different  Si  concentrations  and  the 
concentrations  were  not  reliably  determined.  This  research 
was  undertaken  with  the  goal  of  establishing  a  quantitative 
correlation  between  certain  physical  properties  and  the 
point  defect  structure  and  microstructure  of  Si-doped 
GaAs . 

In  this  note,  we  report  the  results  of  stress-strain 
measurements  and  correlate  these  with  preliminary  results 
of  free-carrier  density  and  LVM  infrared  absorption.  All 
measurements  were  made  on  the  same  material. 

Two  Si-doped  crystal  ingots  grown  by  the  horizontal 
Bridgeman  method  were  investigated.  These  ingots  are 
referred  to  as  KM-16  and  WA-Si.  Various  methods  were 
employed  for  the  determination  of  Si  concentrations  in 
the  ingots.  Initially,  the  atomic  absorption  method  was 
tried  and  comparisons  were  made  with  commercially  available 
Si  standard  solutions.  However,  lack  of  reproducibility 
due  to  solution  aging  effects  and  lack  of  a  dedicated 
piece  of  equipment  made  this  technique  impractical  and 
unreliable.  Methods  of  spectrochemical  analysis  and 
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electron  probe  microanalysis  were  then  employed.  The 

electron  probe  analysis  was  found  to  be  convenient  and 

reliable.  The  data  obtained  by  spectrochemical  analysis 

agreed  with  the  electron  probe  microanalysis  measurements. 

The  electron  probe  microanalysis  measurements  indicated 

19  -3 

that  Si  concentration  varied  from  3  to  6  x  10  cm  from 

the  bottom  to  the  top  of  the  WA-Si  ingot.  The  KM-16  ingot 

showed  much  better  uniformity  with  a  Si  variation  of  less 

than  10  pet.  The  average  Si  concentration  of  the  samples 

obtained  for  mechanical  testing  from  the  KM-16  ingot  was 
19  -3 

2.5  x  10  cm  .  The  average  concentration  of  samples 
obtained  from  a  specific  location  in  the  WA-Si  ingot  was 
4.7  x  1019  cm'3. 

Single  crystal  compression  specimens  with  square  cross- 
section  (5.1  mm  x  5.1  mm  x  10.2  mm)  and  with  long  axis 
perpendicular  to  the  crystal  growth  direction  were  cut 
from  each  ingot  using  a  two-circle  goniometer  and  a 
precision  cut-off  wheel.  The  crystallographic  orientations 
of  the  long  axes  (stress  axes)  of  specimens  cut  from  the 
KM-16  and  WA-Si  ingots  are  shown  in  Fig.  1.  The  Schmid 
factor  (cos  $  cos  0,  where  0  is  the  angle  between  the 
stress  axis  and  the  slip  plane  normal  and  $  is  the  angle 
between  the  stress  axis  and  the  slip  direction)  for 
specimens  from  the  two  ingots  is  clearly  the  same  even 
though  their  orientations  are  different. 
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The  compression  specimens  were  annealed  at  1200 °C 

in  quartz  ampoules  and  then  quenched  in  oil.  The  specimens 

were  aged  at  400,  500  and  700°C  for  30,  90,  900,  7.2  x  10* 

and  1.2  x  10®  minutes  (the  1.2  x  10®  minute  aging  treatment 

was  omitted  at  700°C) .  Details  of  the  annealing  and 

aging  treatments  were  previously  described.® 

Stress-strain  curves  for  all  specimens  were  obtained 

at  400°C  by  loading  in  compression  at  a  constant  force 

rate.  The  stress  rate  during  yielding  was  14.1  MPa  s 

The  yield  stress  values  reported  were  for  a  0.2  pet  offset 

strain.  Details  of  the  stress-strain  testing  procedure 

7  8 

are  reported  elsewhere.  ' 

Figures  2  and  3  show  yield  stress  versus  aging  time 
curves  for  specimens  from  the  KM- 16  and  WA-Si  ingots, 
respectively.  The  yield  stress  plotted  at  t  -  0  is  the 
value  corresponding  to  the  sample  annealed  at  1200°C  for 
one  hour  and  then  quenched  without  subsequent  aging.  With 
the  exception  of  the  Wa-Si  specimens  aged  at  700°C,  the 
yield  stress  increases  monotonically  with  increasing  aging 
time  in  agreement  with  previous  results  by  Swaminathan 
and  Copley.®  They  attributed  this  increase  in  yield  stress 

.  •  •  X 

to  an  increase  in  the  concentration  of  (Siga-Vga)  pairs, 
which  impede  the  motion  of  dislocations.  A  large  difference 
in  yield  stress  also  can  be  noted  between  the  annealed 
and  aged  specimens  from  the  KM-16  and  WA-Si  ingots.  The 
yield  stress  values  for  specimens  from  KM-16  are  lower 
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than  those,  for  specimens  from  WA-Si  at  all  temperatures 

and  for  all  periods  of  aging.  Because  the  Schmid  factors 

for  specimens  from  the  two  ingots  are  equal ,  we  conclude 

that  the  difference  results  from  a  difference  in  critical 

resolved  shear  stress.  This  suggests  that  the  (Si*  -V_  ) 

Ga  Ga 

pairs  were  present  in  all  specimens  and  that  their  concen¬ 
tration  was  greatest  in  specimens  from  the  WA-Si  ingot, 
which  had  the  greatest  Si  concentration. 

Recently,  Chen  and  Spitzer  obtained  values  of  free 

of 

carrier  density  and  concentrations/ various  point  defects 

in  specimens  taken  from  one  of  the  ingots  (KM-16)  used  for 

g 

the  present  study.  The  free  carrier  density  for  each 
specimen  was  measured  by  determining  the  frequency  of  the 
free  carrier  plasma  minimum  observed  in  the  infrared 
reflectivity.  The  concentrations  of  various  point  defects 
were  determined  from  known  LVM  absorption  bands  by  multi¬ 
plying  the  peak  LVM  absorption  less  background,  the  band 
width  at  half  maximum  and  the  absorption  cross-section. 

Figure  4  shows  the  results  for  specimens  annealed  one 
hour  at  1100°C,  quenched  and  then  aged  various  times  at 
400°C.  As  was  previously  observed  by  Kung  and  Spitzer, 
the  decrease  in  free  carrier  density  with  increasing  aging 


time  can  be  accounted  for  entirely  by  the  decrease  in  SiGa 
donors.  Also,  the  other  Si  defect  concentrations,  [Si^sJ , 
l  (SiGa  -  S^  )*],  [367  cm  ^  defect]  and  [369  cm  *  defect] 
remain  nearly  constant  during  annealing.  Chen  and  Spitzer 
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determined  by  transmission  electron  microscopy  that  no 
precipitates  or  dislocation  loops  were  present  in  the  400°C 
aged  specimens.  Based  on  these  observations  Chen  and 
Spitzer  concluded  that  the  decrease  in  [SiGa]  during  the 
400°C  aging  results  from  the  formation  of  electrically 
neutral  (SiGa  ”  VGa)x  pairs. 

If  it  is  assumed  that  the  decrease  in  [Si*  )  results 

ud 

•  *  v 

completely  from  the  formation  of  (Si_  -  V„  )  pairs, 
then  it  is  possible  to  relate  the  pair  concentration  to 
the  yield  stress,  employing  the  equation 

a  -  a  =  A (N  -  6  [Si*  J)1"  (1) 

o  o  Ga 

where  a  is  the  yield  stress,  oq  is  the  yield  stress  for 
the  annealed  and  quenched  specimens,  Nq  is  the  number  of 

•  1  X 

(SiGa  -  VGg)  pairs  for  the  quenched  specimens  and  A  and  m 

are  constants.  Figure  5  shows  that  a  good  correlation  is 

attained  with  Nq  =  1  x  10^^  cm'3  and  m  =  1/4. 

Chen  and  Spitzer  have  pointed  out  that  the  one  fourth 

power  dependence  is  consistent  with  previous  results  of 

Chin  et  al^  for  crystals  of  NaCl,  NaBr,  KCl  and  KBr  doped 

2+  2+  2+ 

with  divalent  impurity  ions  such  as  Ca  ,  Sr  and  Ba  , 

where  it  was  found  that  the  critical  resolved  shear  stress 

varied  with  the  one  half  power  of  the  divalent  ion  concen- 

2+ 

tration.  If,  m  the  case  of  NaCl  doped  with  Ca  ,  we 

assume  that  the  defect  responsible  for  hardening  was 

(Ca*  -  v'  )x,  then  the  concentration  of  the  defect  is 
Na  Na 
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determined  by  the  reaction 


CaNa  +  VNa  "  (CaNa 


V'  )X 
vNa' 


(2) 


Apply  in  the  law  of  mass  action  gives 

<<CaNa  -  VNa)X)  '  VCaNa!<VNa!  <3) 

If  tCa^al  is  large  compared  to  the  number  thermal  vacancies, 
then  [v^  ]  =  [Ca^a]  to  mainta^n  charge  neutrality.  Thus 

<(CaNa  -  Via>Xl  -  VCaNa)2  <4’ 

so  (o  -  oQ)  =  (Caf’ a  «  [  (Caf‘a  -  V '  )  *J  ^  as  we  have  found 
in  the  case  of  the  (Si*  -  V*  )x  pairs. 

Lid  Lad 
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Figures 


Figure  1:  Stress  axis  orientation  of  samples  cut  from  the 
KM- 16  and  WA-Si  ingots  determined  by  the  Laue 
back  reflection  method. 

Figure  2:  Yield  stress  versus  annealing  time  for  samples 

from  ingot  KM-16  for  various  annealing  temperatures 
Figure  3:  Yield  stress  versus  annealing  time  for  samples 

from  ingot  WA-Si  for  various  annealing  temperatures 
Figure  4:  Isothermal  annealing  curves  at  400°C  giving  defect 
or  measured  carrier  concentrations  for  GaAs:Si 
samples  (after  Chen  and  Spitzer) . 

w 

Figure  5:  A  plot  of  a  -  a  versus  (N  -  6  [Si  ])  for 

O  O 

samples  annealed  at  400°C. 
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Stress  axis  orientation  of  samples  cut  from  the 
KM-16  and  WA-Si  ingots  determined  by  the  Laue 
back  reflection  method. 
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5.  Measurements  of  [Si] 

Knowledge  of  the  absolute  Si  concentrations  coupled  with  Hall  effect  and 

optical  data  would  be  helpful  in  assigning  absolute  concentrations  of  Si-Si 

pairs  and  possibly  the  more  complex  Si  defects  in  heavily  Si  doped  GaAs. 

To  be  useful  the  Si  concentrations  should  be  precise  to  within  10%.  This  is 

presently  beyond  the  capability  of  secondary  ion  mass  spectroscopy  (SIMS)  even 

when  used  with  ion  implanted  standards.  Atomic  absorption  spectroscopy  may  be 

18  -3 

feasible  for  concentrations  in  the  mid  10  cm  range  and  higher  -  ones  which 
contain  more  than  simple  Si^a  and  Si^s  defects.  Ideally  one  would  like  an 
analytical  method  capable  of  the  10%  precision  for  GaAs  samples  containing 

17  _3 

5x10  Si  atoms  cm  in  which  the  primary  defect  is  SiGa.  This  would  allow 

the  establishment  of  the  correspondence  between  electrical,  optical  and 

the  independent  analytical  methods  for  establishing  the  total  Si. 

We  proposed  to  develop  such  a  procedure  which  offered  a  10%  precision 

17  -3 

at  Si  concentration  of  1 0  cm  .It  has  not  been  successfully  demonstrated 
but  has  shown  enough  encouragement  to  cause  us  to  continue  its  study  on  a 
no  cost  basis.  The  approach  was  to  dissolve  the  GaAs:Si  samples  in  a  Sn  solution 
containing  a  known  concentration  of  0  at  600°C .  Upon  dissolution  the  Si 
reacts  with  the  dissolved  0  to  form  S.C^s).  The  amount  of  Si  is  estimated 
from  the  charges  in  the  0  concentration  in  the  Sn  solvent.  Thermodynamically 
the  approach  is  sound.  Oxygen  concentrations  were  to  be  determined  by  measurinq 
the  0  fugacities  by  making  the  Sn  solvent  one  electrode  of  a  calcie  stabilized 
zirconia  which  also  serves  as  the  solution  container.  A  pt  reference 
electrode  in  a  controlled  atmosphere  completed  the  cell  .  The  emf  of  the  cell 
is  proportional  to  the  log  of  the  ratio  of  the  oxygen  fugacities  at  each 
electrode.  The  0  concentrations  is  then  calculated  using  Henry's  l.aw.  Thus 
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the  Henry  law  constant  was  to  have  been  measured. 

An  apparatus  was  built  for  allowing  the  insertion  of  the  GaAs  samples 
into  the  Sn  electrode  after  removing  surface  oxides.  A  major  complication 
developed  in  that  the  emf's  of  the  cells  were  not  as  stable  as  those  in 
which  Ga  was  the  solvent  and  upon  which  the  approach  was  originally  based. 

A  simplified  cell  has  been  constructed  with  provision  for  controlling  the  0 
fugicity  in  the  gas  phase  above  the  Sn.  The  object  is  to  establish  the  stability 
limits  of  the  cell.  This  phase  should  be  completed  by  the  end  of  this  summer. 

An  extra  incentive  is  that  the  above  method  should  be  applicable  to  precisely 
measuring  the  oxygen  concentrations  in  GaAs.  Lack  of  such  information  has  lead 
to  much  speculation  about  the  role  of  0  in  semi  insulating  GaAs  substrates. 
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